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PREFACE 


The  conclusions  given  in  this  report  are  based  upon  the  "Engineering 
Manual  (PMlOO-1)"  method  for  culaculation  of  "protection  factors".  Since 
an  error  analysis  is  not  presently  available,  the  conclusions  should  be 
regarded  as  tentative,  pending  the  development  of  such  an  analysis.  In 
addition,  a  redistribution  of  fallout  and/or  changes  in  the  7-ray  spectrum 
emitted  by  the  fallout  may  introduce  further  uncertainties  into  these  con¬ 


clusions  . 
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This  is  Volume  II  o£  four  separately  bound  volumes  that  report  the  research 
completed  in  fulfillment  of  Office  of  Civil  Defense  Work  Unit  No.  3233B,  "Radio¬ 
logical  Recovery  Requirements,  Scructures,  and  Operations  Research."  This  describes 
six  supporting  studies  all  previously  reported  to  the  Office  of  Civil  Defense  in 
research  memoranda.  Volume  I  describes  the  general  aspects  of  the  investigations 
and  presents  the  conclusions  and  recommendations.  The  abstract  for  each  of  the 
volumes  is  presented  on  the  following  pages. 
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ABSTRACT  FOR  VOLUME  I 


This  study  examines  the  apiplication  of  decontamination  strategies  to 
extensive  urban  areas.  Urban  areas  of  various  sizes  (from  a  few  acres  to  an 
interconnected  system  involving  hundreds  of  acres)  are  examined  with  regard  to 
decontaminating  vital  sections  and  their  connecting  links.  Feasible  creation  of 
decontaminated  "islands"  or  marshalling  areas  is  determined.  The  nature  and 
scope  of  command  and  control  system  elements  required  for  conducting  effective 
decontamination  in  practical  situations  are  defined  together  with  the  preattack 
and  postattack  data  required  by  such  a  system.  Several  models  of  detector 
response  to  gamma  radiation  developed  during  the  course  of  the  project  are  briefly 
discussed. 
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ABSTRACT  FOR  VOLUME  II 


Volume  II  contains  six  supporting  studies.  Described  are  a  number  of  models 
for  determining  the  cost  and  effectiveness  of  decontaminating  imralclpal  areas. 

One  study  examines  the  system  components  of  a  command  and  control  system  for 
municipal  decontamination.  These  studies  cover  the  following  subjects: 

1 .  A  Feasibility  Study  of  the  Application  of  Analog  Computers  to  the  Analysis 
of  Decontamination:  This  study  describes  the  research  completed  on  the  applica¬ 
tion  of  analog  computers  to  the  analysis  of  decontamination.  A  simplified  analog 
model  of  the  effect  of  a  single  decontamination  effort  on  a  single  plane  of 
contamination  is  explained  and  described.  Sample  analog  records  of  output  from  a 
prototype  of  this  model  are  presented.  The  design  of  a  more  elaborate  model 
requiring  substantially  more  analog  equipment  is  described,  and  applications  are 
indicated . 

2 .  A  Circular  Model  for  Approximating  Gamma  Ray  Intensity  at  a  Single 
Detector  Location:  An  approximate  procedure  for  determining  the  gamma  ray 
intensity  at  a  point,  due  to  fallout  radiation,  is  investigated.  A  reasonably 
accurate  model  using  circular  annulus  geometry  was  developed  wlilch  Includes  the 
effects  of  gamma  ray  attenuation,  build-up,  backscatter,  and  skyshlne.  This 
model  was  based  on  NBS  Monograph  42. 

A  single  annulus  model  Is  described  In  order  to  show  clearly  the  basic 
premises  upon  which  the  full  "circular"  model  Is  based.  The  development  of  the 
multiple  annulus  model  Is  treated  in  a  manner  to  show  how  each  of  the  parameters 
involved  is  handled.  An  example  demonstrating  the  use  of  the  model  Is  Included. 

^ .  A  Square-Grid  Model  for  Approximating  Gamma  Ray  Intensity  at  a  Sl..Kle 
Detector  Location:  This  study  describes  a  si:nple  and  practical  procedure  for 
determining  the  approximate  fallout  ganma  ray  Intensity  at  a  point  as  a  function 
of  the  geometry  of  the  contributing  planes  and  related  shielding.  The  procedure 
described  employs  a  "square-grid"  technique  for  modeling  the  contaminated  planes. 
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This  model  is  shown  to  be  reasonably  accurate  and  easily  adapted  to  practical 
situations  with  the  help  of  scaled  overlays  for  city  maps  to  be  used  in  conjunc¬ 
tion  with  the  tables  and  graphs  developed  in  the  research.  An  example  analysis 
is  carried  out  to  illustrate  the  method. 

4.  A  Point  Source  Model  for  Approximating  Gamma  Ray  Intensity  at  a  Single 
Detector  location  (The  Equivalent  Planes  Method) :  A  method  for  analyzing  the 
dose  contributions  from  plane  areas  of  fallout  is  developed;  each  contaminated 
plane  is  treated  as  a  weighted  point  source  of  radiation  intensity.  Only  the 
area,  the  location  of  the  center,  the  eccentricity  of  the  approximating  rectangle, 
and  the  intervening  shielding  need  be  considered.  Limits  are  developed  for  the 
eccentricity  of  an  area,  centered  over  the  detector,  which  can  be  represented  by 
a  square  of  the  same  Area  with  an  error  of  less  than  ten  percent  in  dose  rate 
contribution.  Limits  are  also  developed  for  the  area  of  off-center  contaminated 
areas  such  that  the  product  of  area  tima<s  the  dose  rate  contribution  per  unit 
area  is  within  ten  percent  of  the  true  contribution  from  the  area. 

A  sample  analysis  is  given. 

A  step  by  step  procedure  based  on  this  model,  called  "The  Equivalent  Planes 
hethod,"  is  presented.  A  small  booklet  which  could  be  used  as  an  "Equivalent 
Planes  Method"  workbook  (together  with  some  work  sheets)  is  included  in  an  envelope 
at  the  back  of  this  volume. 

5.  A  FORTRAN  Program  for  Decontamination  Analysis:  This  study  describes  a 
debugged  and  tested  FORTRAN  computer  program  to  compute  the  ef fectlvene.ss  parameters 
used  to  analyze  municipal  decontamination.  The  program  was  written  in  FORTRAN 
64  to  he  used  on  large  scale  computers  such  as  the  CDC  3600, 

b .  The  Nature  and  Scope  of  Command  and  Control  System  Eletmenta  Required  ior 
Cof.ductlnR  Effective  Decontamination  in  Municipalities}  This  study  serves  to 
determine  the  nature  and  scope  of  command  and  control  system  elements  whicM  are 
required  co  effect  practical  municipal  decontamination.  The  preattack  and 
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postattack  data  requlrementa  ^or  decontamination  are  specified,  ard  the  essential 
components  of  an  Information  system  for  decontamination  are  Identified  and  related. 
The  influence  of  direct  weapons  efiects  on  the  decontamination  system  Is  examined. 


ABSTRACT  FOR  VOLUME  III 


Volume  III  contains  the  cost  and  effectiveness  data  for  decontamination 
analyses  of  sixteen  sites  and  facilities  in  San  Jose,  California.  Costs  are 
measured  in  team-hours  of  effort.  Decontamination  effectiveness  is  measured  in 
terms  of  fractions  of  dose-rate  remaining  at  specified  detector  Iccc-tions  and 
fractions  of  dose  remaining  for  persons  who  perform  functions  requiring  specified 
daily  activity  patterns  at  the  sites  and  facilities  chosen. 


ABSTRACT  FOR  VOLUME  IV 


Volume  IV  contains  the  cost  and  effectiveness  data  for  decontamination 
analyses  of  twelve  sites  and  facilities  in  Detroit,  Michigan.  Costs  are  measured 
in  team-hours  of  effort.  Decontmination  effectiveness  is  measured  in  terms  of 
fractions  of  dose-rate  remaining  at  specified  detector  locations  and  fractions  of 
dose  remaining  for  persons  who  perform  functions  requiring  specified  daily  activity 
patterns  at  the  sites  and  facilities  chosen. 
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.Radiological  Recovery  Requlretceats ,  Structures,  and  Operations  Re«»garch 
Volume  ~I.  Development  of  Analytical.  Computer,  and  Systems  Models  it.  Support 

of  Decontamination  Analysis 

I .  INTRODUCTION 

This  volume  reports  six  studies  performed  under  Office  of  Civil  Defense 
Subtask  3233B,  Radiological  Recovery  Requirements,  Structures,  and 
Operations  Research.  It  is  addressed  to  technical  personnel  concerned  with  the 
planning  of  postattack  recovery  operations  as  sunmarized  and  desexibed  in  Volume 
I  of  this  report.  All  of  the  studies  presented  in  Volume  II  are  concerned  with 
the  development  of  tools  which  can  be  used  to  examine  the  effectiveness  and 
costs  of  decontamination  when  applied  to  accelerating  recovery  of  an  activity  in 
a  postattack  environment. 

When  this  work  began,  it  seemed  that  large  complexes  comprised  of  several 
buildings  could  not  be  analyzed  quickly  or  efficiently  without  modifying  the 
analysis  techniques  that  were  previously  developed  for  single-facility  shielding 
analysis.  Thus,  to  meet  what  then  appeared  to  be  a  requirement  for  efficient 
decontamination  analysis,  a  number  of  models  were  developed  for  approximating 
garma  ray  intensity  at  a  point  due  to  complex  contaminated  plane  configurations. 
Each  of  these  models  is  discussed  in  a  separate  appendix  (Appendices  A  through 
D)  in  this  volume. 

Paralleling  the  development  of  the  analytical  models  and  the  analog  computer 
model,  two  computer  programs  written  in  FORTRAN  for  the  GDC  3690  were  completed 
and  debugged  to  perform  most  of  the  computation  required  to  analyze  decontamination 
operations  as  applied  to  several  sites  and  facilities  selected  from  San  Jose  and 
Detroit,  The  first  of  these  programs  (Reference  1)  was  developed  and  debugged 
under  another  contract.  This  program  computes  the  plane- by- plane  contributions  to 
intensity  at  a  specified  detector  location.  The  second  of  these  programs  computes 
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•  the  effectiveness  parameters  of  the  individual  decontamination  analyses.  This 

program  is  described  in  Appendix  E. 

Lastly,  the  nature  and  scope  of  the  connand  and  control  system  elements 
required  to  affect  practical  municipal  decontamination  are  defined.  The  pre¬ 
attack  and  postattack  data  requirements  are  identified  and  related  to  the  system 
as  a  whole.  The  influence  of  direct  weapons  effects  on  the  decontamination  system 
are  examined.  This  work  is  described  in  Appendix  F.  The  six  supporting  studies 
included  in  this  volume  are: 

’ .  Appendix  A:  A  Feasibility  Study  of  the  Application  of  Analog  Computers 
to  the  Analysis  of  Decontamination 

2.  Appendix  B:  A  Circular  Model  for  Approximating  Gamma  Ray  Intensity  at 
a  Single  Detector  Location 

3.  Appendix  C:  A  Square-Grid  Model  for  Approximating  Gamma  Ray  Intensity 
at  a  Single  Detector  Location 

4.  Appendix  D:  A  Point-Source  Model  and  the  Equivalent  Planes  Method  for 
Approximating  Gamma  Ray  Intensity  at  a  Single  Detector  Location 

5 .  Appendix  E :  A  FORTRAN  Program  for  Decontamination  Analysis 

6.  Appendix  T:  The  Nature  and  Scope  of  Command  and  Control  System  Elements 
Required  for  Conducting  Effective  Decontamination  in  Municipalities 
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Appendix  A 


A  Feasibility  Study  of  the  Application  of  Analog 
Computers  to  the  Analysis  of  Decontamination 


Note;  The  material  in  this  Appendix  was  originally  submitted  to  USNRDL 
as  Research  Memorandum  RM-OU-214-2*. 


*  J.  T.  Ryan.  A  Feasibility  Study  of  the  Application  of  Analog  Computers 
to  the  Analysis  of  Decontamination.  RM-OU-214-2.  Durham,  North 
Carolina:  Research  Triangle  Institute,  Operations  Research  and 
Economics  Division,  1  April  1965 
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Appendix  A 


A  Feasibility  Study  of  the  Application  of 
Analog  Computers  to  the  Analysis  of  Decontamination 

I .  INTRODUCTION 

To  the  author's  knowledge,  only  digital  computer  programs  are  presently  used 
to  calculate  the  protection  factor  associated  with  specific  detector  locations. 
These  programs,  with  some  modifications,  can  be  used  to  calculate  the  individual 
intensity  contributions  from  each  of  the  contributing  contaminated  planes.  Thus, 
these  programs  can  be  effectively  used  to  estimate  the  value  of  decontaminating 
one  or  more  of  these  planes  of  contamination.  These  programs  are  very  detailed. 
The  inputs  are  an  itemized  description  of  all  of  the  planes  of  contamination  and 
of  the  intervening  shielding. 

The  over-all  time  required  to  prepare  such  inputs  is  usually  very  long,  and 
thus  reduces  the  number  of  facilities  or  activities  which  can  be  considered  in 
the  analysis  of  decontamination  for  a  large  municipal  area.  This  deficiency  led 
to  the  investigation  of  an  analog  model  that  uses  continuous  variables  for  the 
environmental  parameters  as  well  as  tho  parameters  associated  with  the  decontami¬ 
nation  operations.  Analog  models  are  characteristically  smaller  in  scope  than 
digital  models,  less  minute  in  detail,  have  a  shorter  running  tinve,  and  can  be 
more  easily  changed. 

This  study  is  a  result  of  a  very  brief  effort  to  determine  the  applicability 
of  analog  coinpnters  to  the  analysis  of  decontamination.  This  appendix  presents 
a  detailed,  non-technical  description  of  the  analog  computer  study.  The  initial 
pages  are  devoted  to  a  description  and  discussion  of  analog  coraputerr.  The  sub¬ 
sequent  pages  describe  both  the  limited  decontamination  simulation  fitted  to  the 
available  hardware  at  the  Research  Triangle  Institute  and  a  more  elaborate  model 
which  requires  substantially  more  equipment.  Sample  runs  are  Included.  The  final 
section  of  the  appendix  presents  conclusions  and  recoraniendatlons  for  future  work. 
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Section  VI  io  a  technical  description  of  the  mathematical  equations  and  hard 


ware  components  (both  existing  and  proposed)  of  the  decontamination  models. 
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II.  DISCUSSION  OF  ANALOG  COMPUTERS* 

A.  General 

The  history  of  computers  is,  in  fact,  two  histories.  Most  authors  (References 
A-1  and  A-2  for  example)  recognize  that  there  are  very  basic  factors  which  dis¬ 
tinguish  between  the  family  of  digital  computers  and  the  family  of  analog  com¬ 
puters.  Usually  the  abacus  and  the  slide  rule  are  used  to  exemplify  a  simple 
computer  from  each  of  these  families  respectively.  For  the  purpose  of  this  report, 
however,  it  is  not  necessary  to  elaborate  on  the  philosophy  or  historical  divergence 
of  the  two  families  of  computers.  It  is  only  necessary  to  discuss  the  one  primary 
factor  which  distinguishes  between  analog  and  digital  computers--the  way  in  which 
data  is  handled  in  each  machine. 

In  the  digital  computer,  data  is  comprised  of  discrete  numbers  (represented 
by  digits,  magnetic  senses,  bit  positlous — or  in  the  case  of  the  abacua--bead 
positions,  etc.)  whereas  in  the  analog  computer  data  is  represented  as  a  continuous 
variable  (represented  by  wheel  positions,  electric  voltages--or  in  the  case  of  the 
slide  rule--slide  bar  positions,  etc.). 

The  general-purpose  electronic  analog  computer  is  the  type  of  computer 
considered  in  this  report.  In  order  to  evaluate  the  features  of  this  type  of 
analog  computer,  it  is  helpful  also  to  consider  electronic  digital  computers. 

The  following  several  paragraphs  extracted  from  Reference  A-3  very  clearly  present 
the  pertinent  differences  between  analog  and  digital  models  es  they  affect  the 
analysis  and  conclusions  presented  In  this  report: 


Readers  already  fatolliar  wltn  analog  computers  should  omit  this  section.  Readers 
who  would  like  to  learn  more  about  the  use  of  analog  computers  are  referred  to 
Reference  A-1. 

For  a  detailed  (technical)  comparison  of  digital  and  analog  codfriters .  see 
Reference  A-2,  Chapter  2. 
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'‘Digital  computers  are  composed  of  a  large  number  of  multiple  inter¬ 
connected  bistable  devices.  These  on-off  components  along  with  trans¬ 
ferring  devices  form  the  basis  of  the  digital  computer  capability — 
they  enable  formation  and  storage  of  numbers,  simple  arithmetic  opera¬ 
tions,  and  two-valued  logic  operations.  Any  function  which  can  be 
expressed  or  approximated  as  a  sequence  of  these  operations  can  be  per¬ 
formed  on  the  digital  computer. 

"Analog  computers  are  not  as  simply  composed  but  consist  of  many 
components  with  different  functions  or  operating  modes.  Their  capa¬ 
bility  depends  on  their  having  a  large  and  varied  collection  of  com¬ 
ponents  to  draw  on.  Electronic  circuits  provide  this  required  variety, 
including  sumjusrs,  integrators,  multipliers,  function  generators, 
switches,  comparators,  and  many  simple  circuits  V7hich  may  be  built 
quickly  from  resistors  and  capacitors  to  provide  special  functions. 

"Mechanizi'.g  an  analog  computer  model  consists  in  connecting  these 
various  components  together  in  the  same  way  as  the  equations  are  formed 
cj  as  the  physical  system  is  envisioned.  The  model  is  built  in  parallel; 
for  each  operation  In  the  system  there  is  a  component  in  the  model  that 
perfu  ms  the  same  operation.  If  there  are  ten  multiplications  in  the 
system,  then  there  must  be  ten  multipliers  in  the  model,  which  will  all 
be  operating  at  once.  This  is  just  the  opposite  of  digital  computers, 
which  are  organized  ir  series.  That  is,  every  operation  in  the  system 
is  perforioed  by  one  component.  If  there  are  ten  multiplications  in  the 
system,  then  the  digital  model  will  perform  one  multiplication  at  a 
time  in  the  programmed  sequence. 

"One  consequence  of  parallel  organization  is  that  increases  in 
system  size  or  complexity  bring  a  corresponding  increase  in  the  number 
of  components  required  in  the  model.  With  series  organization,  an 
increase  in  system  size  or  complexity  results  in  a  longer  running  time 
for  the  model. 

"A  second  consequence  of  the  parallel  organization  is  that  if  the 
model  is  interrupted  at  any  time  during  the  run,  all  variables  will  be 
found  to  be  at  the  value  corresponding  to  the  one  the  system  would 
have  at  that  time.  That  is  to  say,  there  Is  an  isomorphism  between 
time  in  the  system  and  time  in  the  model.  This  is  not  the  case  with 
series  organization  as  in  digital  computer  models,  where  the  model  can 
be  interrupted  only  at  certain  specific  times  to  be  meaningful. 

"Computing  time  on  au  analog  computer  model  is  usually  very  short. 
While  the  basic  switching  time  of  the  digital  computer  may  be  one 
thousand  to  ten  thousand  times  faster  than  the  response  tine  of  an 
aq^alog  computer  component,  the  number  of  operations  required  in  the 
dil^htal  computer  for  any  complicated  function  is  so  large  that,  in 
fact,  the  total  operation  may  not  be  faster  than  the  analog  computer 
component.  Thus,  speed  In  performing  the  total  operation  coupled 
with  parallel  organization  allows  the  time  variable  on  the  analog  to 
be  scaled  as  a  fraction  or  multiple  of  real  time.  Running  times  on 
the  analog  are  usually  scaled  to  range  between  10  seconds  and  2 
minutes.  For  example,  an  analog  model  of  a  valve  which  opens  In  2 
seconds  might  be  scaled  to  have  1  second  of  real  time  equal  10  seconds 
of  computer  time,  thereby  slowing  do'm  the  operation  for  each  observa¬ 
tion.  A  battle  model  in  which  the  time  of  flight  to  target  Is  5  hours 
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might  be  scaled  to  have  1  hour  of  real  time  equal  1  second  on  the 
computer,  vastly  speeding  up  the  time  required  for  a  run, 

"Analog  computer  models  require  a  considerably  more  contracted 
treatment  of  large-scale  problems  but  provide  a  more  tractable  model. 
That  is,  large  volmties  of  data,  readily  handled  by  digital  computers, 
are  impossible  on  the  analog  unless  they  can  be  treated  in  an  aggre¬ 
gated  manner;  but  changes  can  be  made  easily,  while  operating,  by 
simply  resetting  a  dial. 

"There  is  considerable  difference  between  the  form  of  the  output 
on  the  two  computers.  The  usual  output  of  a  digital  computer  is  a 
tabulation  of  the  results.  Analog  o<itput  is  usually  in  a  graphic 
form,  whether  as  a  chart  showing  the  variation  of  a  number  of  para¬ 
meters  with  time  or  as  a  plot  of  one  variable  against  another.  This 
graphical  output  is  adapted  very  v;ell  to  problems  where  a  physical 
understanding  of  the  interactions  is  desired,  since  the  effect  of 
variations  can  be  immediately  seen. 

"In  summary  then,  an  analog  computer  model  with  its  treatment  of 
data  and  variables  quite  aggregated  would  be  expected  tc  be  considerably 
smaller  in  scope  than  a  digital  computer  model.  Operational  changes 
could  be  simply  made  and  input  data  easily  varied.  Output  would  be 
graphical  and  running  time  short,  so  that  results  of  input  changes 
would  be  immediately  available  and  easily  interpreted.  This  is  in 
contrast  with  present  digital  models  which  are  detailed,  longer  run¬ 
ning,  less  readily  changed,  and  considerably  broader  in  scope." 


B.  Analog  Computers  and  Decontamination  Analyses 

Some  comments  augmenting  and  tying  the  above  remarks  to  the  problem  of 
analyzing  decontamination  operations  are  appropriate  here. 

The  "parallel  organization"  feature  described  above  is  particularly  useful  in 
analyzing  the  effects  of  decontaminating  a  number  of  independent  contaminated 
surfaces.  Instead  of  having  to  compute  the  effects  of  the  separate  operation 
serially,  one  can  build  the  model  so  as  to  show,  as  a  function  of  (scaled)  real 
time,  the  effect  of  the  sum  of  these  decontamination  operations,  where  they  could 
be  performed  simultaneously,  or  sequenced  in  any  way  whatsoever. 

Furthermore,  different  differential  equations  can  be  used  to  govern  the 
efficiency  of  the  individual  decontamination  operations  where  the  parameters 
describing  these  differential  equations  can  be  varied  manually  at  run  time. 

This,  of  couise,  would  be  impossible  to  do  with  a  digital  computer.  The  ability 
to  easily  solve  very  complex  differential  equations  (not  easily  soluable  by  paper 
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and  pencil  methods)  is  one  o£  the  chief  capabilities  cf  analog  computers. 

The  "graphical  output"  feature  can,  of  course,  be  performed  on  a  digital 
computer  with  the  appropriate  peripheral  equipment,  but  the  ability  to  respond  at 
contputer  running  time  to  the  information  being  graphed  is  impossible  with  normal 
digital  equipment.  Thus,  analog  coooputers  can  be  a  useful  tool  for  training 
personnel  responsible  for  large-scale  municipal  decontamination.  This  will  be 
discussed  in  more  detail  later  in  this  appendix. 


ill.  DESCraPTION  OF  II.-nXED  DEC(»frAMl^ATIO^■  SIMULATION 


A .  Introduction 

The  limited  decontamination  analog  model  was  developed  to  explore  the  possi¬ 
bilities  arid  determine  the  application  of  analog  computers  to  decontamination 
analysis.  As  such,  this  model  does  not  provide  realistic  answers  to  decontamina¬ 
tion  problems  and  only  roughly  simulates  any  real  situation.  Very  gross  treatment 
of  the  details  was  used  in  order  to  simulate  the  main  features  of  a  single  decon¬ 
tamination  operation  on  a  small  sized  computer. 

B.  The  Hardware  Used 

The  limited  analog  model  described  in  this  section  was  fitted  to  a  basic 
installation  of  the  model  3400  Desk  Top  Donr.er  Analog  Computer  (Figure  A-1) .  The 
basic  model  3400  computer  contains  ten  amplifiers,  usable  as  summers  cr  integrators. 
The  Research  Triangle  Institute  installation  also  includes  a  Donner  Model  3430 
Problem  Board  and  a  Donner  Model  3073  Poteatiosieter  Strip.  The  problem  board 
(shown  in  Figure  A-2)  is  used  to  interconnect  coiiponents  for  the  solution  of 
particular  problems  and  the  potentiometers  are  used  to  adjust  constant  and 
parametric  coefficients  in  the  equations  which  describe  particular  problems.  The 
output  from  the  analog  model  is  graphed  using  an  Offner  Dynograph  Amplifier- 
Recorder  Model  542.  The  recorder  is  also  shown  on  che  right  in  Figure  A-I.  A 
more  detailed  description  of  the  actual  compments  used  n  t;  is  model  is  con¬ 
tained  in  section  VI. 

C.  The  Limited  Decontamination  Model 

In  the  limited  decontamination  analog,  each  contaminated  plane  which  con¬ 
tributes  to  the  intensity  at  a  detector  location  is  approximated  by  a  point  source 
at  the  centroid  of  the  contaminated  plane. 

Figure  A-3  is  a  schematic  -.bowing  the  relationships  between  a  single  detector 
and  four  contaminated  planes. 
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Figure  A-1 

Model  3400  Desk  Top  Dotmer  Analog 
Computer  Installation 


Figvire  A-2 

A  Close-up  of  the  Problem  Board 
Set  up  to  Analyze  a  Simple 
Decontamination  Operation 


Contamitnatcd 
Plane  #  1 


Figure  A-3 

Schematic  Diagram  of  Simplified  Decontamination  Operations 
with  Four  Contaminated  Planes. 
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The  contributions,  C^,  to  radiation  intensity  from  ground  level  sources  of 
contamination  are  com{>uted  using  Equation  A-1. 


"i  "  T  ^ 


(A-l) 


where  *  the  area  in  square  feet  of  the  i  contaminated  plane, 

r^  =  the  distance  in  feet  from  the  centroid  of  the  i*"^  contaminated  plane 
to  the  defector,  and 

=  the  shielding  factor  associated  with  the  i^^  contaminated  plane. 

Hera  represents  the  fraction  of  the  total  intensity  (received  at  the  detector) 

whi;h  comes  from  the  plane.  The  shielding  factor  is  a  dimensionless  number 

which  attenuates  the  contribution  according  to  the  shielding  between  the  detector 

and  the  i^^  plane.  (Note  that  0  <  <  1) . 

The  fraction  of  the  total  intensity  received  at  the  detector  which  comes  from 

the  roof  of  the  building  in  which  the  detector  is  centrally  located  is  found  by 
*■ 

Equation  A-2 


=  n  In  (1  +  —  ) 


(A-2) 


where  r  ~  the  radius  of  a  circular  roof,  with  an  area  equal  to  the  area  of  the 

actual  roof  over  the  detector, 

h  =  the  height  of  the  roof  above  the  detector,  and 

S  =  the  shielding  factor  associated  with  the  roof, 
r 

Equation  A-2  can  also  be  used  to  determine  the  contribution  to  the  intensity  at  an 
unshielded  detector  (out-of-doors)  from  the  plane  directly  below  the  detector. 

Here  =  1  and  h  is  a  height  corresponding  to  the  ground  roughness  factor. 


See  page  743  of  Reference  A-7  for  a  derivation  of  Equation  A-2. 

"kif 

See  Reference  A-5  for  an  explanation  of  how  the  ground  roughness  factor  affects 
the  "effective"  height  of  the  detector. 
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In  th«  simplified  model  both  — n-  S  and  IT  In  (1  +  — y  )  S  are  treated  as  single 

r  ^  h  ^ 

constant  inputs.  Any  variations  in  these  individual  parameters  can  be  simulated 
by  appropriately  changing  the  capacitors  which  represent  these  aggregate  variables 
before  running  the  model. 

** 

Decontamination  is  governed  (for  each  individual  plane)  by  Equation  A-3. 


*  * 

Fi  =  Fi  +  (1  -  F^)  e 


(A.3) 


th 


where  F^^  =  the  fraction  of  fallout  remaining  on  the  i"  contaminated  plane  after 
it  has  been  decontaminated. 

t  h 

F.  =  th(5  limiting  fraction  of  fallout  remaining  on  the  i  contaminated 
plane  after  infinite  decontamination. 

K  =  a  constant  associated  with  a  given  method  and  the  physical  nature  of 
i 

the  i^^  contaiiiinated  plane. 

E.  -  the  effort  (usually  measured  in  man-hours)  which  is  applied  to  the 
i*"^  contaminated  plane. 

Hnre.  both  1  -  F  and  K  E  are  treated  as  single  variable  inputs.  Both  of 

j  ^  ^  - 

these  variables  can  be  varied  by  changing  settings  of  the  potentiometer  which 
is  assigned  to  the  particular  variable  pirameter.  This  can  be  done  while  the 
model  is  being  run. 

Table  A-I  is  a  sumr.Tary  of  the  equations  used  in  the  limited  model  of  decon¬ 
tamination. 

While  the  model  is  being  run,  the  individual  contributions  are  summed  and 
the  radiation  intensity"  at  the  detector  is  output  as  a  function  of  time.  The 
total  dose  than  an  individual  would  receive  at  the  detector  location  is  also  out¬ 
put  as  a  function  of  time.  The  dose  is  measured  over  the  time  interval  from  t^ 
to  t^.  This  is  simply  given  by 


** 


See  Reference  A-6  for  a  derivation  of  Equation  A-3. 
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H(t)  denotes  the  value  (represented  by  a  voltage  on  the  analog  computer)  associated 
with  the  sum  of  the  intensity  contributions  from  the  individual  planes  as  a 
function  of  time.  For  a  system  with  a  number  of  instantaneous  decontamination 
operations,  H(t)  will  be  a  monotonically  decreasing  step  function. 


A-12 


TABLE  A-T. 


System  Equations  fee  Simplified  Model 


1.  Intensity  contribution  from  the  1  ‘  ground  level  contaral-nated  plane: 

Af 

C.  =  — r  3.  where 
1  2  1 

"i 

til 

A,  *  the  areas  in  square  feet  of  i  contaminated  plane. 

r^  “  distance  in  feet  from  centroid  of  fallout  on  i*^  contaminated 
plane  to  the  detector  location. 

th 

S^  =  shielding  factor  associated  i  contaminated  plane. 


2.  Intensity  contribution  from  roof: 

2 

C  =II  ln(l+-^)S  where 


r  =  the  radius  of  a  circular  roof  of  equal  area, 
h  *  height  of  roof  (above  detector) . 

=  shielding  factor  associated  with  roof. 

3.  Decontamination  efficiency: 

F  =  F*  +  (1  -  F*)  e'^^i^i  where 
i  i  i 

th 

F^  =  fraction  of  fallout  remaining  after  decontaminating  i  plane 

•k 

F.  =  fraction  of  fallout  which  cannot  be  removed. 

1 

K.  =  constant  associated  with  i*"^  surface  and  the  method  used  to 
decontaminate  it. 

t  h 

E^  =  the  amount  of  effort  applied  to  decontaminate  the  i  contami 
nated  plane.  (Usually  measured  in  man-hours  of  effort). 


For  the  simplified  model  the  following  parts  of  the  above  equations  were  set  to 
single  constant  parameters; 


i 

1. 

"l 


2.  =  n  In  (1  +  ^^) 

2  J  ^ 


The  following  parts  of  the  above  equation  were  set  to  single  variable  parameters 


i  * 

1.  V^=l-F^ 


2. 
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IV. 


SAMPLE  RUNS  WITH  PROTOTYPE  MODEL 


Several  test  runs  vere  made.  Seven  ?  these  runs  are  included  here  to 

illustrate  the  nature  and  format  of  the  answer.s  obtainable.  The  situation 

simulated  in  all  seven  runs  was  the  decontamination  of  a  single  contaminated 

-1  2 

plane.  For  simplicity,  ao  fallout  build-up  function  was  used.  The  t  ’  decay 

curve  was  approximated  as  the  sum  of  two  exponential  decay  curves  (i.e., 

“1.2  “K  t  “K  t 

t  ‘  ^e  1  -i-e  2).  A  fuller  discussion  of  this  approximation  is  contained 
in  Chapter  VI. 

In  all  of  these  sairple  runs,  no  attempt  was  made  to  quantitatively  scale 
dose-rate  and  total  dose.  The  time  variable  was  also  left  arbitrary.  Thus,  two 
inches  of  horizontal  distance  on  the  graph  outputs  might  correspond  to  two  weeks, 
two  hours,  etc.  The  fiCtual  graph  outputs  from  these  seven  runs  are  presented 
as  Figures  A-4  through  A-10.  Time  runs  from  right  to  left  on  these  graphs.  For 
Figures  A-4  through  A-9,  decontamination  is  assumed  to  take  n  small  finite  time. 

In  Figure  A-10  decontaiiination  is  instantaneous. 

Figure  A-4  shows  the  effect  of  deconcatainating  a  single  plan  of  contamination 
where  90%  of  the  fallout  material  has  been  removed.  Figures  A-5  and  A-6  show 
the  effect  of  decontaminating  the  saiue  plans  of  contamination  where  the  decontami¬ 
nation  efforts  expended  are  only  75%  and  50%  of  the  effort  expended  in  the  opera¬ 
tion  simulated  by  the  run  which  produced  Figure  4.  (Sell]  assuming  the  percent 
of  removable  fallout  equals  90%). 

Figure  A-7  shows  the  effect  of  decontaminating  the  same  plane  where  only 
50%  of  the  fallout  can  be  removed  and  a  "full"  effort  (the  same  effort  is  assumed 
for  Figure  A-4)  is  expended.  Figure  A-8  shows  the  effect  of  decontaminating 
this  plane  where  only  50%  of  "full"  effort  is  expended. 

For  all  of  the  Figures  A-4  through  A-8,  both  total  dose  ..nd  dosa-rate  are 
shown  where  total  do‘?e  begi  is  to  accutnulate  at  time  tj  and  stops  accurjj  lat  Ing  at 
time  t^-  Note  that  the  vertical  scale  for  total  dose  in  Figures  A-7  and  A-8  is 
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exactly  half  of  the  scale  used  for  total  dose  in  Figures  A-4,  A-5,  and  A-6, 

In  Figure  A-9,  total  dose  begins  to  accumulate  at  t  where  no  build-up 
function  is  assumed.  The  total  dose  continues  to  accui.iulate  throughout  the  graph 
on  F)  ’ure  A'9.  Here,  the  horizontal  scale  for  total  dof-j  was  again  cut  to  half  of 
thft  used  in  Figure  A-8,  The  same  fraction  of  removable  fallout  and  effort 
expended  are  used  in  Figures  A-8,  A-9,  and  A-10. 

Figure  A-10  shows  how  the  graphs  would  look  if  decontamination  took  place 
instantaneously.  Otherwise  Figure  A-10  is  the  same  as  Figure  A-9. 
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Analog  Record  of  Dose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  I 


Legend 

t  *  instantaneous  time  of  arrival  of  fallout  (no  buildup), 
a 

t,  •  time  when  decontamination  begins. 

G 

(ti*  t2)  *  interval  of  time  for  which  total  nose  is  calculated. 


Thus, 


t  <  t 

1 


<  t  < 


t 


2 


t  >  t 


2 
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Figure  A-5 

Analog  Record  of  Dose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  il 


A- 18 


Figure  A-6 

Analog  Record  of  Pose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  III 


Legend 

t  = 
a 


time  of  arrival  of  fallout  (no  buildup) . 
time  when  decontamination  begins . 

t^)  =  interval  time  for  which  total  dose  is  calculated. 


Thus , 


Dose 


t  <  t 

1 


<  t  <  t2 


t  <  t 


2 
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Figure  A- 7 

Analog  Record  of  Dose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  IV 


d 

t  =  time  of  arrival  of  fallout  (no  buildup) . 
t^  =  time  when  decontamination  begins. 

(ti,  ^2)=  interval  of  time  for  which  total  dose  is  calculated. 


Analog  Record  of  Dose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  V 


Legend 


t  =  time  of  arrival  of  fallout  (no  buildup) . 
a 

t,  =  time  when  decontamination  begins, 
d 

(ti,  t^)  =  interval  of  time  for  which  total  dose  is  calculated, 


Thus , 


Dose  = 


dx 


dx 


t  <  t. 


t  >  t. 
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Figure  A-9 

Analog  Record  of  Dose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  VI 


Legend 

t  =  time  of  arrival  of  fallout  (no  buildup) . 
a 

t,  =  time  when  decontamination  begins, 
d 

Total  dose  is  calculated  from  t^  to  any  time  t  . 

Thus, 


0 


t  <  t 
—  a 


Dose  = 


t 

H(x)x 

a 
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t  <  t  <  08 
a 
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Figure  A-10 

Analog  Record  of  Dose  Rate  and  Total  Dose  as  a  Function  of  Time 

for  Case  VII 


Legend 

t^  =  time  of  arrival  of  fallout  (no  buildup). 

tj  =  time  when  decontamination  begins. 

Total  dcse  is  calculated  from  t  to  any  time  t  . 

a 


Thus , 


t  <  t 

—  a 


Dose  = 


/. 
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H(x)x  dx 


c  <  t  <  “ 
a 
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V.  AN  ElABORATE  MODEL 


Although  quite  adequate  to  demonstrate  feasibility,  the  limited  analog 
model  described  in  Section  III  is  not  suitable  for  complex  decontamination  prob¬ 
lems.  A  more  elaborate  model  with  less  aggregation  of  variables  i.tnd  a  capacity 
for  handling  a  larger  number  of  contaminated  planes  is  necessary.  The  basic 
system  equations,  however,  remain  the  same  for  the  larger  analog  model  described 
in  this  section. 

As  in  the  model  already  implemented  and  described  in  Section  III,  decontami¬ 
nation  is  measured  by  a  "graph"  of  dose-rate  (radiation  intensity)  at  some  detector 
location  of  interest.  Total  dose  is  measured  for  selected  time  intervals  at  the 
same  point.  Complex  buildup  functions  can  be  used  to  simulate  the  initial  stages 
of  fallout  arrival.  The  number  of  contaminated  planes  which  can  be  considered  for 
a  given  detector,  as  well  as  the  number  of  detector  locations,  is  in  theory, 
unlimited.  A  practical  limit  on  the  numbers  of  analog  components  would  limit  the 
number  of  contaminated  planes  to  20  and  the  number  of  detector  locations  to  three. 
By  appropriately  using  relays  and  switching  circuits,  the  total  dose  received  by 
a  person  changing  environments  instantaneously  (such  as  moving  from  indoorr  to 

outdoors)  could  be  simulated. 
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All  functions  (log,  t  ‘  ,  etc.)  can  be  quite  accurately  simulated  by  using 
suitably  adapted  function  generators.  This  will  improve  the  precision  of  the 
simulation. 

Furthermore,  a  quick-response  capability  could  be  included  to  serve  as  an 
analysis  and  training  tool.  This  is  accomplished  by  using  ranual  over-rides 
which  permit  the  operator  to  apply  decontamination  efforts  on  the  basis  ol  graphs 
which  portray  previous  decontamination  efforts  during  the  same  run. 

A  technical  description  and  flow  diagram  of  the  more  elaborate  model  is 


included  in  Chapter  VI. 
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VI.  A  TECHNICAL  DESCRIPTION  OF  THE  ANALOG  DECONTAMINATION  MODELS 


A.  Calculations  to  be  Performed 

The  calculations  to  be  performed  by  the  computer  are  expressed  in  Equations 


A-5  and  A-6; 


n  Aj  S . 

R;t)  =  2  1  -  u  (t  -  t.  ) 


-<t  -  t  ) 

c  e  - 

"i  ^i 


+  Jt  I  In 


D(t)  =  /  R(t)dt  =  /  v(t;  t^.t^)  R(t)dt 


t-l-‘ 

1  - 


(A-5) 


(A-6) 


where  the  symbols  are  defined  as  follows: 

R(t)  =  dose  rate  in  roentgens  per  hour, 

D(tJ  =  dose  (normalized  units)  in  roentgens, 

n  =  number  of  contaminated  planes  considered, 
th 

A^  =  area  of  i  contaminated  plane  in  square  feet, 

l^h 

r^  =  distance  from  detector  to  centroid  of  i  contaminated  plane  in  feet, 

th 

=  shielding  factor  for  i  location  (C  <  <  1) , 

th 

K.  “  constant  for  i  contaminated  plane  and  a  particular  decontamination 

^  method, 

E^  =  effort  applied  to  i*"^  coitaminated  plane  in  man  hours/1000  sq  ft, 

“K  E 

C„  =  effort  coefficient  =  1  -  e  i  i, 

£ , 

1 

th 

C  =  fraction  of  removable  mass  for  i  contaminated  plane, 
m. 

1 

T.  =  decontamination  time  constant, 


effective  radius  of  surrounding  contanrlnated  plane  (roof  or  detector 
surface) , 

effective  height  of  surrounding  contaminated  plane  (roof  or  detector 
surface) , 

shielding  factor  for  surrounding  contaminated  plane  (roof  or  detector 
surface) , 


A->/ 


i 


“  Integration  limits  for  dose  D(t) . 

The  function  u(t-t^  )  ic'  a  unit  step  function  which  has  the  values: 

i 

u(t-t  )  ■=  0  0  <  t  <  t 

^i  “  ^1 

»  1  t  <  t  <  CO 

^i  “ 

The  functioi  v(t, ;t,,t«)  has  the  values: 

1  JL 

vCt.t^.t^)  “0  0  <  t  <  t^ 

=  1  tl  ^  ^  ^2 

-2  tj  <  t  <  00 

A  functional  diagram  showing  the  calcjlations  to  be  performed  is  shown  in  Figure 
A-11.  Only  one  of  the  n  identical  channels  (each  channel  corresponding  to  a 
single  contaminated  plane)  is  shown.  For  n  areas,  the  calculations  represented 
within  the  dotted  outline  will  be  repeated  ti  times. 

B .  Range  of  Constants,  Parameters,  and  Variables 

The  units  of  the  various  quantities  and  the  expected  ranges  of  these  quantitit 
are  listed  below: 


SYMBOL 

UWITS 

RANGE 

R(t) 

r/hr 

- 

t 

hrs 

2A-336  (Iv  OO) 

n 

dimensionless 

1-10 

A 

ft' 

100-I00,u00 

I 

ft 

10-500 

S 

dimensionless 

0-1 

CL 

L 

dimensionless 

0-1 

K 

2 

ft  /man-hour 

.1-5 

A- .^3 


I'thiT  Char.no Is 


Figure  A-11 

Functional  Diagram  of  Calculacions 


SYMBOL 

UNITS 

RANGE 

E 

man  hrs/1000  ft^ 

.01-1 

C 

dimensionless 

0-1 

m 

-1 

T 

hr 

10-0.1 

t. 

hr 

10-1000 

^i 

P 

ft 

10-200 

h 

ft 

1-100 

S 

dimensionless 

0-1 

s 

^2 

hr 

10-1000 

hr 

10-1000 

D(t) 

roentgens 

- 

For  computational  purposes  the  dose  rate  R(t)  and  dose  D(t)  are  measured  in 
normalized  units.  To  convert  to  roentgens  per  hour  or  roentgens,  one  point  on  the 
dose  rate  curve  muse  be  hand  calculated,  and  the  proper  scale  factor  determined. 

C.  Computer  Scaling 

1.  Choice  of  Time  Scale 

The  time  range  of  major  interest  is  from  Ik  to  336  hours.  For  convenience, 

the  computer  t'  ’  scale,  t  ,  may  be  chosen  so  that  one  hour  of  real  time, 

c 

t^,  corresponds  to  ,5  second  of  computer  time.  Time  zero  on  the  computer 
will  correspond  to  time  24  hours  in  real  time.  Thus: 

=  ~  (c^  -  24)  seconds,  (A- 6) 

where  t^  ==  computer  tinu-  in  seconds,  and 

t  =  real  time  ^n  hours, 
r 

2 ,  Dose  Rate  Scale 

At  t^  =  24  hrr,  the  dose  rate  R(t)  will  be  asstr  led  the  valui  100  in 
arbitrary  units.  The  computer  scale  factor  for  R(t)  will  be  taken  as  1 
volt/unit.  To  convert  from  arbitrary  units  to  roentgens  per  hour,  the  scale 
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factor  is  determined  by  hand  calculation  of  the  dose  rate  at  t  -  2k  hrs. 

r 

The  dose  D(t)  will  be  scaled  at  .01  volt  per  arbitrary  unit  to 
prevent  the  integrator  from  overloading  with  long  integration  times. 


3.  Coefficient  Scaling 

For  adjusting  values  of  the  various  parameters  that  remain  constant 
during  a  particular  computer  run,  the  computer  potentiometers  will  be  used. 

In  order  to  simplify  the  equipment  requirements,  it  is  convenient  to  grcmp 
certain  parameters  into  dimensionless  ratios  where  possible.  These  parameter 
groups  will  t!  be  set  into  the  computer  as  potentiometer  settings.  Addi¬ 
tional  coefficients  are  defined  as  follows: 


C  =  rt  In 
s 


range  .01  -  10 


range  0  -  20 


For  potentiometers  with  dial  divisions  from  0  to  100,  the  following  scale 
factors  are  assigned  to  the  parameters  that  are  to  be  set  on  potentiometer 
dials: 


Coefficient 

Range 

Sc^ie  Factoi 

(units  per  dial 

C 

"i 

.01-10 

.1 

C 

s 

0-20 

.2 

^i 

0-1 

.01 

C 

m 

i 

0-1 

.01 

^i 

0-1 

.01 

S 

s 

0-1 

.01 
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It  should  be  noted  that  with  additional  equipment  complexity,  each  value  sucn 
as  A^,  etc.  may  be  entered  separately  as  a  dial  setting,  and  the  necessary 
multiplications  and  divisions  accomplished  by  the  computer. 

4.  Dynamic  Range  Requirements 

In  some  cases,  it  would  be  desirable  to  simulate  a  time  period  longer 
than  15  days.  In  Increasing  the  time  period  of  a  single  computer  run,  diffi¬ 
culties  may  be  encountered  because  of  tue  small  signal  levels  that  will  be 
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present.  For  example,  the  function  t  *  decreases  by  a  factor  of  approxi¬ 
mately  100  between  24  hours  and  1100  hours.  If  the  computer  output  were 
initially  set  at  100  volts,  the  output  at  1100  hours  (550  seconds  computer 
time)  would  be  approximately  1  volt  if  no  decontamination  had  taken  place. 

This  level  is  not  too  small;  however,  the  signal  level  from  each  of  the  i 
channels  must  be  considerably  less  than  1  volt  (with  10  channels  giving  equal 
contributions,  the  signal  level  would  be  .1  volt).  On  certain  computers, 
the  noise  level  prevents  accurate  operation  at  these  low  voltage  levels. 

In  order  to  avoid  the  inaccuracies  inherent  in  using  low  signal  levels, 
the  best  approach  would  appear  to  be  to  break  the  problem  into  time  periods 
of  inter'est.  If  the  period  of  interest  were  from  10  to  20  days,  the  computer 
time  would  be  scaled  such  that  the  problem  starts  at  t^  =  240  hrs  instead  of 
24  as  discussed  previously. 

D.  Computer  Diagram 

A  computer  diagram  to  perform  the  necessary  calculations  is  shown  in  Figure 

A-13  where  the  symbols  are  defined  in  Figure  A-12.  A  function  generator  generates 
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a  voltage  proportional  to  t  '  .  This  voltage  is  adjusted  to  a  value  analogous  to 

the  contribution  of  each  of  the  contaminated  planes  by  the  potentiometers  C  ,  S  , 

C  ,  and  S  .  The  voltages  analogous  to  the  radiation  received  from  each  area  are 
s  s 

summed  in  a  summing  amplifier,  and  the  resultant  voltage  is  recorded  on  a  strip 
chart  -.ecorder. 
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Cymbol 


Name 


Operation 


Coefficient 


e  =  K  e,  (0  <  K  <  1) 
o  1  “ 


High  Gain  Amplifier  e^  »  k 


k  >  10^ 


e^  =  energized 
e^  »  o  unenergized 


Function  Generator  *^0  *  ^  '1^ 


Figure  A-12 

Table  of  Flow  Chart  Symbols  and  Operations 


dt  +  E 

o 
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Computer  Diagram 


To  simulate  decontamination  of  a  given  area,  the  relay  Ry  is  closed,  A 

*.  f 

*  U  /  'T 

time  lag  of  the  form  e  i  is  introduced  by  operational  amplifier  #7,  and  the 
original  voltage  is  reduced  by  the  proper  amount  by  subtraction  in  amplifier  #1. 

The  amount  of  voltage  subtracted  depends  upon  the  effort  coefficient  and  the  per¬ 
cent  of  removable  material  as  set  on  the  potentiometers  C^,  and  . 

To  calculate  the  dose  from  the  dose  rate,  an  Integrator  is  energized  by  relay 
Ry  B.  This  relay  remains  closed  for  a  time  period  t2  -  as  set  on  the  sequencing 
function  generator. 

Note  that  the  computer  diagram  shows  only  one  of  n  identical  setups  that 
must  be  patched  into  the  computer. 

E ,  Function  Generators 

-1,2 

To  simulate  a  function  varying  with  time  as  t  ,  several  methods  are  avail¬ 
able.  A  conventional  diode  function  generator  with  a  linear  time  function  input 
may  be  used.  Another  method  is  to  use  an  x-y  plotter  as  a  function  generator, 

supplying  a  linear  time  function  as  a  drive  for  the  x  axis,  and  reading  out  a 
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voltage  proportional  to  t  ‘  on  the  y  axis  by  means  of  a  curve  following  probe. 
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The  simplest  technique,  however,  appears  to  be  that  of  approxistating  the  t 

function  by  a  sum  of  exponentially  decaying  functions.  For  example,  with  three 

exponentials,  the  t  curve  may  be  fitted  by  the  sum  at  six  points  on  the  curve. 

The  exponential  functions  are  gevierated  very  easily  on  the  analog  computer  by 

setting  up  an  initial  charge  on  a  capacitor,  and  discharging  the  capacitor  slowly 
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through  a  resistor.  Figure  A- 14  shows  the  approximation  of  the  t  '  function 
in  the  range  24-350  hrs  by  the  sum  of  two  exponentials  as  follows: 

[.86  .-0313(^-24)  ^  ^-.00375(6-24) 

45.3  I 

The  computer  diagram  for  accomplishing  this  approximation  would  appear  as  shown 
in  Figure  A-15  (.5  sec  computer  time  ■  1  hr  real  time). 
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1.0 


.2 


.1 

.9 

.8 


.04 


.01 


.001 


Time  Hours 


Figure  A-14 

Approximation  of  the  Function  of  the  Form  t 


by  Two  Exponentials 
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Quantity 

11 


Item 


Relays 

1  Sequence  Function  Generator 

(11  Channel) 

1  Dual  Channel  Recorder 

A  40  amplifier  analog  computer  should  be  sufficient  for  the  10  area  simulation. 

The  sequence  function  generator  may  be  construe, ed  In  the  laboratory  If  a  suitable 
type  cannot  be  purchased.  To  buy  this  equipment  would  cost  about  $40,000.  Total 
rental  costs  (including  set-up  costs)  would  be  under  $1,000  for  a  one  year  period. 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Conclusions 

The  basic  conclusion  of  this  study  is  that  the  analog  computer  provides  a 
potentially  useful  and  quite  feasible  tool  for  simulating  and  analyzing  dccontami" 
nation  operations. 

The  development  and  implementation  of  this  first  analog  model  of  decontamina¬ 
tion  has  pinpointed  the  areas  of  difficulty  and  the  direction  in  which  further  work 

appears  most  promising.  The  most  difficult  problem  encountered  was  accurately 

-1.2 

approximating  the  t  ‘  decay  rate.  This  problem  is  considered  in  detail  in 

Chapter  VI;  since  a  high-accuracy  analytic  function  generator  would  be  very  costly 

(about  $1000),  some  form  of  approximation  involving  the  sum  of  exponential  functions 
-K  t 

of  the  form  K^e  2  is  a  practical  compromise  solution. 

Even  though  a  number  of  the  problems  associated  with  the  simple  modv-^l  using 
only  one  plane  of  contamination  were  a  result  of  complex  scaling  necessary  to  keep 
the  range  of  variables  within  machine  limits,  future  models,  to  be  useful,  must 
include  even  more  detail  and  so  must  use  even  more  equipment.  Thus,  scaling  will 
beco;ne  an  Increasingly  difficult  problem. 

B.  Future  Work 

Any  future  work  must  be  vndertaken  in  the  light:  of  the  major  difficulty  found 
in  using  the  limited  a-.i&log  model--the  large  amount  of  equipment  required  for  an 
adequate  simulation.  Nonetheless  it  is  recommended  that  two  major  directions  be 
explored ; 

The  development  of  a  larger  more  detailed  simulation,  but  still  with 
simplified  aggregation  of  variables  where  possible,  end 
2.  The  development  of  a  quick  response  capability  to  serve  as  a  ready 
analysis  and  training  tool. 
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Appendix  B 


A  Clrciilai;  Model  for  Annrov limiting 

Intensity  at  a  Detector  Location 


Note:  The  material  In  this  Appendix  wa;  orglnally  submitted  to  USNBDL 
as  Research  Memorandum  RM-OU-214-3*. 


. . — 

W.  R.  Oavl»,  M.  K.  Moss,  and  A.  J.  Benjamin.  A.  Circular  Model  foe  Aporoxl- 
Bating,  Gaapi  Bay  Intaniitv  at_a  GiVfa  kuttor  Locarion.  RH-OU>2L4-3. 
lAirhas,  North  Carolina;  Research  Triangle  Institute,  Operations  Re.scarch 
and  Economics  Olvlslon,  30  April  1%5. 
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Intercept  value  used  in  the  analytical  approximation  o£ 
L  (X,u))  functions. 

Vertical  height  of  the  detector  above  the  source  plane. 


Total  detector  response  at  a  standard  reference  position. 

Gamma  dose  rate  received  at  a  detector  from  the  annular 
region. 

Gamma  dose  rate  received  at  a  detector  from  the  annular 
sector  of  the  annulus. 


D/D  *  Dose  reduction  factor, 

o 

h,  .  =  Height  of  contaminated  annulus  above  the  ground  or  reference 

^  plane. 


L  (X) 

L  (X.q)  - 

\  (X.cu)  “ 


Total  detector  response  for  infinite  plane  isotropic  source. 

Generalized  expression  for  the  geometry  factor. 

Geometry  factor  describing  detector  response  due  to  circular 
plane  sources  of  fallout  radiation. 

Fractional  area  of  the  k^^  annular  region  covered  by  the 
contaminated  annular  sector. 

Total  number  of  sectors  in  the  k^*^  annulus. 


S 


0) 

X 


Slope  value  used  In  the  analytical  approximation  of  L  (X,a}) 
functions. 

Solid  angle  fraction. 

Effective  mass  thickness. 

Height  of  annular  region  above  the  detector. 


A  Circular  Model  for  Approximating  Gamma  Rav 
Intensity  at  a  Detector  Location 


1.  INTRODUCTION 

The  purpose  of  this  appendix  Is  to  show  the  Initial  results  of  recent  studies 
directed  toward  the  development  of  rather  simple  and  practical  procedures  for 
determining  the  gacxna  ray  Intensity  at  a  point  due  to  nuclear  weapon  Induced 
fallout.  The  consensus  of  those  Involved  In  this  research  effort  was  that  If 
such  an  approach  proved  feasible,  It  would  result  in  a  considerable  reduction 
In  both  labor  and  computer  costs  for  those  cases  where  apprc'.'*’nate  results  are 
considered  adequate.  In  addition,  the  development  of  such  an  analytical  tool 
would  be  of  considerable  value  In  simplifying  those  tasks  associated  with  the 
ensuing  analysis  of  decontamination  effectiveness  In  municipal  areas.  This 
research  is  a  first  step  in  the  development  of  several  alternative  approaches  to 
a  model  for  field  analysis  of  fallout  decontamination  effectiveness. 

In  an  effort  to  assure  a  reasonable  degree  of  accuracy  in  the  analytical  model 
presented  in  the  following  discussions,  such  pertinent  factors  as  gamma  ray  attenu¬ 
ation,  build-up,  backscatter,  and  skyshine,  are  incorporated.  Another  basic 
premise  closely  adhered  to  in  connection  with  the  development  of  this  model  Is 
that  it  be  built  around  the  familiar  results  contained  in  NBS  Monograph  42 
(Reference  B-l). 

The  general  notatlonal  scheme  presented  here,  together  with  selected  elementary 
geometric  considerations  and  the  information  contained  in  Reference  B-l,  provides 
the  basis  for  what  is  called  a  "circular  model".  This  model  is  utillaed  to  furnish 
values  of  the  so-called  "reduction  factor,"  D/D^,  doe  to  various  planes  of  con¬ 
tamination  at  different  heights  with  or  without  the  inclusion  of  barrier  shielding. 


In  particular,  these  planes  of  contamination  are  divided  into  annular  sectors  wit., 
the  detector  located  in  each  instance  at  the  origin  of  the  polar  coordinate  system. 

The  approach  taken  in  describing  this  circular  model  is  that  of  developing 
initially  the  detector  response  for  the  single  contaminated  annulus  case  (Section 
II)  following  with  increasingly  more  complex  annular  geometries  (Section  III)  and 
terminating  with  a  brief  discussion  on  barrier  sheilding  considerations  (Section  IV). 

Each  of  the  annular  geometric  configurations  utilized  in  the  analytical 
development  of  the  model  required  the  extensive  use  of  the  L  (X)  and  the  L(X,w) 
functions  developed  in  Reference  B-I.  It  was  found  during,  the  course  of  the 
investigation  tha^  the  L  (X,u))  functions  could  be  approximated  over  their  range 
of  greatest  utility  by  simple  analytical  expressions  (Section  VI).  In  turn,  such 
expressions  can  be  quickly  and  easily  evaluated  with  the  aid  of  a  desk  calculator 
and  standard  mathematical  tables,  or  efficiently  scored  in  a  computer  memory  without 
the  necessity  of  lengthy  subroutines. 

The  procedure  for  treating  attenuation  through  barriers  makes  effective  use 
of  the  L  (X)  function  and  the  (X,u))  geometry  factor,  through  judicious  choices 
and  combinations.  It  should  be  noted  that;  considerable  refinement  of  these 
procedures  and  associated  calculations  may  be  necessary  for  situations  near  the 
detector  location  at  the  origin  of  the  coordinate  system. 
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II.  SINGLE  ANNULUS  MODEL 


The  single  annulus  is  chosen  initially  since  it  will  permit  the  development 
of  the  basic  premises  from  which  more  complex  geometries  can  later  be  described. 


A.  Gomoletd  U.iiforml|Y  Contaminated  Annulus 

In  order  to  enpr'-sb  the  detector  response  for  the  case  of  a  single  uniformly 
contaminated  annulus  (see  Figure  B'l),  it  is  necessary  to  first  define  the  solid 
angle  fractions  and  subtended  at  the  detector  by  the  contaminated 

annulus. 

1.  Solid  Angle  Expressions 

From  Figure  B-1  we  define  first  the  solid  angle  fraction 
subtended  by  angle  as 


coj^a  *  1  (Se®  Section  V  for  (?''rivation  of  solid 

angle  expressions' 


where 


cos 


V(ka)^  + 


Similarly, 


“(k+U.  •  1  “(k+1). 


where 


cos  a 
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2.  Detector  Response 


The  detector  response  for  the  single  annulus  can  be  computed  by 


appropriate  differencing  of  the  geometry  factor  L^(X,u))  (see  Figure  B-2 
and  Section  VI)  for  the  values  of  and  differencing  can  be 

Illustrated  In  terms  of  the  detector  response  fov  the  uniformly  contaminated 
annulus  In  Figure  B-l  as  follows: 


^  -  L  (X) 

o 


[ 


I<c  (^•“(kH)a^  *  ^ 


(B-5) 


where  Is  the  detector  response  at  a  standard  reference  position*,  and 
Is  the  response  due  to  the  annulus.  The  function  L  (X)  Is  graphed 
In  Reference  B*l. 

It  should  be  noted  that  all  dlstancis  d  from  the  detector  to  the 
contaminated  plane  of  Interest  may  be  expressed  In  terms  of  effective  mass 
thickness  (X)  for  the  medium  between  the  source  and  detector.  For  the  case 
of  air  (at  20°C,,  76  cm  Hg)  as  this  medium,  the  following  relationship  may  be 
written  (see  p.  19,  Reference  B-l): 


d/X  -  13.3  Ft.  of  alr/psf. 


(B-6) 


llie  Interchangeability  of  d  and  X  as  the  penetration  variable  In  these 
analytic'’*!  expressions  Is  permlssnble  with  the  use  of  the  conversion  factor 
indicated  In  Equation  B-b. 

Single  Annular  Sector  Case 

The  a’utular  sector  shown  in  Figure  B>i  represents  a  fraction 

•  •  e 


2* 
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ti  defined  here  as  the  detector  response  at  an  unprotected  position  in  an  infinite 
hoaofaneoue  eir  ecdlun  loceced  three  feet  i^ove  e  hypothetical  Infinite  p lance  source. 
Thie  source  is  ceneidered  to  be  <f  the  same  character  as  tha  radiation  fallout 
unif  >mly  dietributed  on  the  tonteminetion  planc(s}  being  inveetigeted. 


C«oiiii»try  Factor  (L^  (X,w})  D«picting  D«t«ctor  B«spon«a 


Du»  «.o  Ci'cvir.r  Pi*nc  C 


Radiation  Source 


o£  the  total  annulus.  Consequently,  the  detector  response  due  to  the  annular 
sector  In  the  annulus  is  simply 


D 

o 


1.  (X) 


"■c  <*-‘"(k+l)a>-*'c 


<B-e) 


C.  Multiple  Annular  Sectors 

For  the  case  of  N  annular  sectors  comprising  all  the  contaminated  areas 
ch 

of  the  k  annulus,  the  detector  response  is  given  as  follows: 


j-N 
L  C. 


4  S 

M 

d"  D 

0  0 


j-N 


L  (X) 


\  (*»“(k+l)a^  ■  ^c 


(B-9) 


where  (see  Section  li'B)  represents  that  fraction  of  the  annulus  covered  by  the 
uniformly  contaminated  annular  sector. 
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III.  MULTIPLE  ANNUUR  REGION  MODEL 


The  following  discussion  outlines  the  expansion  of  the  previously  developed 
concepts  foi  a  single  uniformly  contaminated  annulus  to  that  of  the  more  frequently 
encountered  multiple  annular  region  case. 

The  notation  utilized  In  developing  mathematical  expressions  for  this  model 
will  relate  to  the  physical  description  shown  In  Figure  B~3.  The  use  of  the 
subscript,  j  ,  which  refers  to  the  (here)  nonexistent  annular  sectors,  Is  not 
required,  but  is  retained  for  consistency  In  notation. 

A .  A  Set  of  Uniformly  Contaminate  ■  Annular  Regions 
1.  Detector  Response 

First  it  is  necessary  to  determine  the  detector  response  due  to  the 
annular  region.  This  quantity  Is  given  by  the  expression 


L  (X) 


(B'lO) 


t 


k 


it  Che  height  of  Che  k^^ 


annular  region  above  Che  detector. 
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Figure  B-3 

Hultlpl*  Annular  Kcgions  tottflguraf  Shoving 
Unifomly  Contaalnatad  B«*gion3  of  Arbitrary  Haight 
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2 .  Contaminated  Regions  of  Arbitrary  Height 

With  the  annular  region  detector  response  defined,  we  next  explore 
a  condition  which  has  considerable  practical  usage  when  attempting  to  deter 
mine  detector  response  in  urban  areas  where  radioactive  fallout  can  be 
considered  uniformly  distributed.  It  Is  desired  to  compute  the  detector 
response  to  N  annular  uniformly  contaminated  regions  of  arbitrary  height. 

The  detector  response  from  a  set  of  such  annular  regions  Is  given  by 


Here  it  is  assumed  that  each  annular  region  is  not  effectively  "shadow" 
shielded  by  another.  A  more  generalized  analysis  concerning  mutual  shielding 
and  barrier  attenuation  will  be  given  in  Section  IV. 


As  a  means  of  illustrating  the  summing  of  the  detector  responses 
associated  with  equivalent  and  uniformly  contaminated  annular  planes  of 
different  h^^^  values,  the  following  example  is  included. 


From  the  information  provided  In  Figure  an  expression  for  Che  total 

detector  response  can  be  written  as  follows: 


11) 
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Figure  B-4 

Detector  Reaponee  Due  to  Uniforaly  Conteaineted 
Plcnce  vith  Different  h..  Veluee 
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where  D  «  D 
°1  °2 


B.  Equivalent  and  Uniformly  Contaminated  Annular  Sectors  for  a  Set  of  Annular 
Regions. 

The  detector  response  which  is  attributed  to  the  j  annular  sector  of  the 
annular  region  may  be  defined  by 
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M.  L  (X) 


'■( 


d-h.  d-h,  . 
kj  _ 111 


V 


(k+1)^  a^  + 


V 


(ka)^  + 


J  .  (B-13) 


where  is  again  the  fractional  area  of  the  k^^  annular  region  covered  by  the 
contaminated  annular  sector. 

It  is  then  apparent  that  the  detector  response  of  the  k^^  annulus  is 

found  by  summing  the  expressions  as  follows  where  we  let  be  the  total 

o  th 

number  of  sectors  in  the  k  annulus: 


^  : 

D  - 

o 


D  ■^"**k  r 


‘'kj'^^Ck+Da^  ■  \ 


w-\j-  “k.)] 


(B-K) 


B-IJ 


By  knowing  the  detector  response  for  the  annulus,  it  is  now  possible  to 
express  the  detector  response  for  any  number  of  such  annular  regionj  by 


MjL  (X) 


h  "^(k+Da^'  (‘^■‘^hk'  “ka^ 


(B- 
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IV.  BARRIER  ATTENUATION 


» 


One  o£  the  most  common  situations  encountered  In  the  analysis  o£  contaminated 
regions  located  In  core  areas  of  large  municipalities  Is  the  frequent  occurrence 
of  physical  barriers  between  the  detector  and  a  contaminated  plane  of  Interest. 

As  an  example  of  a  rather  simplified  approach  to  the  barrier  problem,  we  will 
consider  here  the  case  of  two  annular  sectors  of  separate  annular  regions.  The 
heights  of  the  respective  annular  sectors  are  assumed  to  be  such  that  one  Is 
completely  shielded  by  the  other.  This  particular  case  Is  Illustrated  In  Figures 
B-5  and  B-6. 

The  gamma  radiation  dose  at  the  detector  from  a  shielded  contaminated 
annula..  sector  Is  essentially  due  to;  (1)  direct  gamma  radiation  passing  through 
the  barrier;  and  (2)  scattered  gamma  radiation  transported  by  various  Interaction 
processes  In  air  to  the  detector  location.  Here  for  simplicity  the  actual  radiation 
dose  from  a  shielded  contaminated  annular  sector  will  be  approximated  by 
(see  Figures  B-6  and  B-7) 


L  (X) 


L 

c 


(X.  u^.) 


■ 


(B-16) 


where  X  is  the  effective  distance  through  the  barrier.  This  approximation  will 

not  be  found  to  properly  treat  the  scattered  radiation  contribution  to  the  detector. 

The  functions  L  (X),  L  (d,a;),  and  L  (X,u))  limit  the  detector  response  to  a 

c  c 

circular  area  of  '’nifom  gamaia  radiation  intensity  whose  center  is  on  the  perpen¬ 
dicular  from  the  detector  to  the  pi  'ne  and  separated  by  either  a  thickness  of  air  d 
or  barrier  thickness  X  (see  i^igure  8-7).  A  numerical  example  of  these  considera¬ 
tions  is  given  in  Section  VII. 


8-15 


Figure  B-5 


Two  Annular  Sectors  of  Separate  Annular  Regions 
One  Sector  Shielded  From  the  Detector 


Figure  B-6 

Detector  Response  From  a  Shielded  Sector 
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Gcoactry  Asvociatcd  With  (X,u)  «nd  (D»u}  FacCO's 
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V.  der:vation  of  solid  angle  fraction  expression 

The  purpose  of  introducing  the  concept  of  solid  angle  in  this  memorandum  is 
due  primarily  to  the  fact  that  the  chosen  direction  for  viewing  the  geometrical 
configurations  associated  with  the  circular  model  is  not  limited  to  that  of  a 
plane.  As  a  result  the  more  generalized  notion  of  angle,  i.e.,  solid  angle,  is 
utilized  in  a  somewhat  modified  manner. 

One  of  the  imist  common  means  of  envisioning  a  solid  angle  is  to  draw  a  cone 

frou  a  point  as  the  apex.  Then  the  ratio  of  the  area  Intercepted  by  that  cone  on 

any  sphere  centered  at  the  apex  to  the  square  of  the  radius  of  the  sphere  is  a 

measure  of  the  solid  angle  Q  subtended  by  the  cone.  When  measured  in  this 

manner  the  unit  of  solid  angle  is  called  the  steradian.  Since  the  total  area  of 

2 

a  spherical  surface  is  4i(R  ,  the  solid  angle  subtended  by  any  surface  which 
completely  surrounds  a  point  cannot  exceed  Arc  steradians.  It  is  then  aisc 
apparent  that  the  maximum  solid  angle  subtended  on  one  side  of  a  plane  surface 
is  2n  steradians.  On  the  basis  of  this  premise,  the  unit  of  angular  measurement 
00  adopted  in  this  memorandum  was  chosen  to  be  2x  steradians.  This  unit  is 
consistent  with  that  found  in  Reference  B-1,  where  a;  is  designated  as  the  "solid 
angle  fraction"  in  an  effort  to  avoid  confusion  in  angular  terminology. 

By  using  the  paramatars  indicated  in  Figure  B’B,  the  derivation  of  the  basic 
solid  angle  fraction  expression  found  in  Section  ll'-A*l  of  this  appendix  can  be 
obtained . 

With  the  geosMtry  defined  as  illustrated  above,  an  exprashion  for  the  area 
intarcepced  on  tht  sphtrical  surface  can  be  written  as 
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Figure  B>8 

Solid  Angle  Geoaetry 
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2 

k  ^  t  2jf  (R  sin  a)  R  da  =  -2nR  cos  a 


/ 


2jtR  (1  -  cos  Oj^) 


(B-17) 


From  the  basic  definition  of  a  solid  angle,  the  following  relationship 
is  obtained. 


“k.  ■  7  ■  “ka> 


(B-18) 


Since  our  unit  of  solid  angle  is  2n  steradians,  Equation  B-18  must 
be  rewritten  in  the  following  manner: 


2n 


1  -  cos 


“ka  *  “k. 


(B-19) 


where  is  the  so-called  "solid  angle  fraction." 
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VI.  APPROXIMATE  ANALYTICAL  FIT  FOR  L^  (X,a>)  FUNCTION 


It  was  discovered  during  the  course  of  developing  the  circular  model  that  an 
analytical  fit  over  a  rather  extensive  range  of  values  for  the  L^  (X,oj)  curve 
found  on  page  41  of  Reference  B-l  is  possible  by  employing  an  expression  of  the 
general  form: 

D/D^  -  L  OL,ta)  “  cu®e^  .  (B-20) 


The  parameter  S  is  obtained  from  the  previously  referenced  curves  by 
applying  the  following  familiar  analytical  relationship. 


In  L  (X,a),)  -  In  L  (X,u),) 

g  a  I  I  iifc  I  ...  —I.— 1 1 1, .  ■ 

In  -  In 


(B-2L) 


The  parameter  A  is  then  calculated  by  solving  Equation  8-20  for  A  as 
follows : 

A  a  . L  ,CX»a).  . In  oi 

^  X 


(B-22) 


The  best  analytical  fit  for  the  L^  (X,«)  curve  was  achieved  by  obtain¬ 
ing  the  L^  (X,cu)  values  given  on  page  41  of  Reference  B-1  corresponding  to 
and  values  of  0.05  and  0.2  respectively  and  solving  for  the  intercept 
at  cu  ■  0.2  for  each  X  value  selected.  The  numerical  results  of  these 
c^ilculations  are  presented  in  Tabic  B-I,  and  the  resultant  curve  is  available 
as  Figure  B-2  in  this  appendix. 

It  is  apparent  when  comparing  the  L  (X,a>}  curves  developed  here  with  those 

appearing  in  Reference  fi-1  that  the  accuracy  of  the  L  (X.cjd)  values  obtained  rrom 

Equation  B-20  deteriorate  rapidly  due  to  the  lack  of  convergence  as  This 

is  particularly  true  for  those  values  of  cu  >  0.3  and  cu  >  0.3  associated  with  the 

(X,  )  function.  It  should  also  noted  that  there  are  no  rcstricvlons  associated 

with  this  model  that  would  prohibit  Che  exclusive  use  of  exact  values  (X,  )  as 
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TABLE  B-I 


Calculated  Values  of  S  and  A  for  Selected 
Effective  Mass  Thickness  Values 
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VII .  EXAMPLE 


In  an  effort  to  demonstrate  the  type  of  situations  best  handled  by  the 
circular  model,  and  at  the  same  time  effectively  summarize  the  analytical 
methods  presented  in  this  appendix,  the  following  example  is  included. 

The  problem  is  that  of  determining  the  c.etector  response  (reduction  factor) 
dt^e  to  the  contaminated  annular  sectors  (A,  C)  in  Figure  B-9.  For  the 
purpose  of  illustration,  an  value  of  0.2  has  been  arbitrarily  seL;:!Cted  for 
each  of  the  three  planes  of  contamination  to  be  invescigated. 


Contaminated  Sector  A: 

First  we  define  the  solid  angle  fraction  subtended  by  angle  0!^^  by 


=  1  -  cos 


(B-23) 


where 


cos 


+  (d)^ 


^(50)^  +  (100)^ 


0.893  .  (B-24) 


Therefore,  =  0.107. 


We  can  now  convert  the  distance  in  air  from  the  source  plane  into  an 
effective  mass  thickness  value  by  employing  the  conversion  factor  provided 
in  Equation  B-6  as  follows: 


X 


JLQO.  f-t-..  _ 

13.3  Ft.  of  air/psf 


a  7.53  psf  . 


(B-25) 


The  detector  response  for  this  case  is  given  by  the  expression 


Exaaple  Illustration 
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By  substituting  the  values  of  M.,  X,  and  oj.  into  the  previous  expression 

J 

we  obtain 


D^/D^  -  0.2  L  (7.53) 


(7.53,  0.107) 


(E-27) 


Obtaining  the  (X,u))  value  from  Figure  B-2,  the  detector  response  is 
found  to  be 


D./D  =  0.2  (0.375)  (0.054)  *  0.00405  .  (B-28) 

A  O 


Contaminated  Sector  B: 


The  expression  for  the  detector  response  from  the  annular  sector 
of  the  contaminated  annulus  is  derived  in  Section  III  and  appears  below 
in  only  a  slightly  modified  form. 


Db/D^  •  Mj  I.  a) 


Lc  I  d-h^,  1  -  ^ 


2  2 

r  +  (d-hj^)^ 
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(k)^  +  (d-h^/ 
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In  this  case  the  effective  mass  thickness  for  air  is 


d-h. 


13.3  Ft.  of 
air  psf 


.26  psf 


(B-30) 


With  the  X  value  and  the  Info.'matlon  provided  In  Figure  B-9,  we  can 
proceed  to  find  the  detector  ret|>onse  In  the  following  manner; 
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0.2  L  (2.26) 


Db/D 


o 


-  0.2  L  (2.26) 


L  (2.26,  0,554)  -  L  (2.26,  0.487) 
c  c 


(B-31) 


Here  the  L  (X)  and  (X,u))  values  are  obtained  from  page  33  and  page  41  respectively 
of  Reference  B"!,  since  the  od  values  in  this  case  approach  the  limits  of  accuracy 
for  the  approximation  method  described  in  Chapter  VI.  Finally, 


D„/D  -  0.2  (0.57)  (0.045)  =  0.00513  .  (B-32) 

o  o 


Contaminated  Sector  C; 

In  finding  the  detector  response  for  Sector  C,  a  barrier  is  encountered 
between  the  detector  and  the  contaminated  plane  being  investigated.  From  the 
discussion  given  in  Section  IV,  we  can  express  the  detector  response  for  this 
barrier  case  as  foiiovs: 

D^/D^  •  L  (X)  ^  (X,  u)^2>  '  \ 

The  values  for  and  found  by 

I  -  cos  -  1  -  _ iOQ  •  0.142.  (B-34) 

^(60)*  +  (100)^ 

and 

'-‘ii2  “  ^  ^  _ —  —  “  °  •  (8-35) 

^(80>^  ♦  (100)^ 
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*  The  X  value  chosen  for  the  barrier  is  75.3,  which  is  equivalent  to  0.523  feet 

of  concrete.  By  substitution  of  these  values  into  the  expression,  we 

obtain 

D^/D^  -  0.2  L  (75.3)  ^  (75.3,  0.22)  -  (75.3,  0.142)j  .  (B-36) 

By  using  the  values  for  (X,u))  in  Figure  B-2,  the  detector  response  for  Plane 
C  is: 

Dg/D^  “  (0.2)  (0.038)  0.453-0. 285  j  -  0.00128.  (B-37) 

Despite  certain  limitations  inherent  in  the  circular  model  presented  here, 
the  methods  described  can  produce  simple  and  practical  approximate  solutions  to 
many  types  of  problems  requiring  kiuwledge  of  gamna  ray  intensity,  due  to 
weapon- Induced  fallout,  at  a  single  detector  location. 
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VIII.  SUMMARY  AND  CONCLUSIONS 


From  the  information  provided  in  the  previous  sections,  It  Is  apparent  that 
the  circular  model  Is  somewhat  restricted  In  application.  The  principal  limitation 
Is  that  circular  or  annular  planes  of  contamination  are  a  poor  approxltnation  to 
most  real  situations.  However,  It  may  be  possible  in  many  cases  to  reduce  non- 
circular  geometries  of  Interest  to  circular  regions  of  equal  area.  Of  course, 
such  an  innovation  mu.st  be  approached  with  caution  since  for  certain  configurations 
the  error  Introduced  in  the  detector  response  due  to  radical  changes  in  the 
original  geometry  may  be  quite  significant.  (Of  course,  the  "Engineering  Manual" 
azimuthal  sector  approach  makes  these  approximations.)  It  is  noted,  however,  that 
this  appendix  represents  but  one  of  a  series  of  different  approaches  to  the  over¬ 
all  problem  of  quick  and  efficient  approximations  to  the  gamma  ray  intensities  due 
to  fallout  radiation.  Other  approaches  will  not  have  the  particular  limitations 
noted  here. 
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12.  s 
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14.  X 
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Width  of  unit  annulus. 

Intercept  value  of  approximate  L(X,u;)  functions. 

V 

Dose  reduction  factor  (D  /D  is  the  dose  rediiction  factor  for 
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The  coordinates  of  the  center  of  the  square  in  the  i^^ 
annulus. 


A  SQ' are-Grld  Model  for  Approximating  Gamma 
Rav  Intensity  at  a  Single  Detector  uocation 


I.  INTRODUCTION 

A.  Purpose 

This  appendix  describes  the  development  o£  a  simple  and  practical  procedure 
for  estimating  fallout  gamma  ray  intensity  at  a  detector  as  a  function  of  the 
geometry  of  the  contributing  planes  of  contamination  and  related  shielding.  The 
need  for  such  a  procedure  arises  from  requirements  for  analyses  of  large  municipal 
areas  with  many  contaminated  planes  and  detector  locations  scattered  throughout  the 
area.  Conceptually,  the  "circular-model"  procedure  described  in  Appendix  B  is 
simple  to  use,  but  accuracy  suffers  when  the  contaminated  planes  and  barriers  do 
not  have  approximately  polar  symmetry  about  the  detector  location.  A  square-grid 
model  can  accurately  approximate  real  situations  more  easily  sir.::e  most  planes  of 
contamination  (and  most  buildings)  are  rectangular.  Further,  the  square-grid 
model  can  be  applied  easily  to  an  arbitrary  detector  location  with  the  use  of 
scaled  map  overlays. 

B.  Some  Analysis  Considerations 

A  problem  encountered  during  the  development  of  the  square-grid  model  was 
the  increased  difficulty  in  computing  the  solid  angles  subtended  at  the  detector 
by  the  square  areas  Adjunct  to  this  problem  is  the  subsequent  determination  of 
the  radiation  at  a  given  detector  location  using  Spencer's  L- functions  (Reference 
C-l).  A  simple  sector-type  weighting  as  was  used  in  the  circular  model  is  not 
possible  with  the  square-grid  model. 

The  solid  angles  subtended  at  a  given  detector  location  by  a  square  area  can 
be  computed  rigorously.  For  example,  the  .<!olid  angle  subtended  by  any  square 


C-l 


can  be  computed  (Reference  C-1)  by  the  equation 


(C-l) 


where. 


»i  “ 


^‘i 


Z 


+1  if  corner  1  la  in  quadrant  1  or  3 
-1  if  corner  i  ia  in  quadrant  2  or  4  ’ 

coordinatea  of  corner  1, 

IVil  /  1x^1  . 

1*1  /  1*^1 

perpendicular  distance  from  the  detector  to  the  plane 
containing  the  plane  of  contamination  under  analysis, 
and  directly  above  the  origin,  and 

solid  angle  of  a  rectangle  with  height  z  and  corners  at 

(Xj^.y^).  (*i»"yi)‘ 


A  full  discussion  of  Equation  C~1  and  the  variables  involved  is  given  in  Reference 
C-1.  Unfortunately,  when  computing  the  solid  angles  for  the  situations  considered 
in  this  appendix,  it  was  founa  that  the  T(e,q)  curves  of  Reference  C«1  were  not 
sufficient  for  computing  the  solid  angle  cu.  This  is  because  differences  must  be 
taken  between  numbers  of  approximately  the  saiae  magnitude.  The  r  (e,q)  functions 
thus  were  computed  from  the  analytic  expression 


(C-2) 


C.  The  Model 

The  general  program  which  wan  followed  for  the  development  of  a  model  lor 
approximating  fallout  gamnia  ray  intensities  at  a  detector  is  now  described.  The 
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model  is  primarily  aimed  at  the  approximate  eatimation,  in  the  field,  of  decontamina¬ 
tion  effectiveness.  High  speed  computer  programs  are  not  practical  in  these 
circumstances.  Of  course,  building  and  area  characteristics  must  be  knovm 
approximately.  A  general  discussion  of  the  model  and  the  information  required 
for  its  application  is  presented  next, 

I*  General  Method 

The  over-all  approach  to  the  square-grid  model  for  studying  the 
effects  of  decontamination  of  fallout  radiation  is  based  on  the  following 
criteria; 

(a)  The  techniques  must  enable  local  civil  defense  personnel  in 
the  field  to  calculate  the  effects  of  large-scale  decontamination  of 
municipal  areas; 

(b)  The  approximate  techniques  and  procedures  must  not  require 
computers;  l.e.,  graphical  and  tabular  forms  should  be  applicable. 

As  noted  in  Section  11,  the  square-grid  model  is  used  to  estimate 

the  reduction  factor  0/D  for  the  various  olanes  of  falli'ut  contamination  - 

o 

each  at  differing  heights.  The  planes  of  contamination  are  divided  into 
square-grid  areas,  and  the  detector  may  in  principle  be  located  at  any 
arbitrary  position,  at  height  z.  The  size  of  the  square-grid  sections  was 
chosen  as  SO  feet,  an  area  which  is  large  enough  that  subsequent  computations 
are  not  unwieldy,  but  which  is  also  sufficiently  small  to  yield  an  acceptable 
approximation  of  major  contaminated  planes  in  almost  any  area.  Consequently, 
any  plane  of  contamination  in  this  model  is  approximated  by  50  foot  square 
grids.  An  irotroplc  detector  is  assumed  to  be  located  at  a  height  z  above 
the  center  square  grid.  It  is  at  this  point  detactor  that  the  reduction 
factor  0/D^  of  any  given  contaminated  plane  area  is  computed. 

Since  there  will  in  general  be  several  planes  of  contaminacion  at 
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varying  heights,  contributing  to  the  fallout  ganuna  radiation  at  a  given 
detector  location,  the  analysis  will  functionally  depend  to  a  large  extent 
on  the  variable  z,  the  detector  height.  By  varying  the  detector  height 
above  the  ground  plane,  one  can  generate  functions  which  can  be  used  for 
the  analysis  of  the  contributions  of  planes  of  contamination  at  various 
heights. 

In  this  model,  the  z  dependence  of  the  detector  response  is  treated 
separately  from  the  other  functional  dependences.  Tnis  is  accomplished  by 
separately  considering  the  detector  response  from  planes  of  contamination  at 
each  vertical  height.  The  responses  due  to  planes  of  coi.taminatlon  at  all  heights 
are  then  superimposed  to  obtain  the  total  detector  response  or  reduction  factor 
D/0^  at  a  given  point. 

The  model  will  thus  work  as  follows:  A  given  plane  of  contamination  at  a 
particular  height  is  approximated  by  50  foot  square  grid  areas.  The  detector 
is  always  vertically  centered  on  the  center  of  the  center  square  in  the 
square-grid  network.  The  squares  are  numbered  as  discussed  in  Section  II. 

When  the  detector  response  due  to  a  given  square  is  desired,  the  solid  angle 
subtended  by  Che  square  at  the  detector  and  the  effective  distance  from  the 
square  area  to  the  detector  must  be  computed. 

Once  the  solid  angles  are  determined,  the  detector  response  due  Co  this 
particular  square  may  be  determined  by  use  of  appropriate  curves.  Every  square* 
grid  area  for  the  given  plane  of  contamination  may  be  so  treated  and  the  solid 
angles,  weighting  functions,  and  detector  responses  (reduction  factors)  of 
each  of  these  squares  may  be  tabulated.  This  procedure  may  be  then  applied  to 
all  other  planes  at  another  height,  and  so  on.  For  thc>  square*grid  overlay,  one 

then  has  a  tabulation  of  the  unique  (in  the  absence  of  barriers)  detector  response 

e 

for  each  individual  square*grid  area. 

y~" . . . 

If  barriers  are  present,  one  refers  to  specially  contructed  tables  and  curves 
Chat  will  be  discussed  in  paragraph  3,  Section  11. 
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When  the  square-grld  overlay  Is  appropriately  scaled  to  detailed  municipal 
maps  (for  example,  Sanborn  Maps),  then  the  properly  scaled  overlays  may  be 
placed  over  the  maps  and  the  contribution  to  the  detector  due  to  given  areas 
of  contamination  may  be  quickly  assessed.  This  is  accomplished  by  noting  the 
square-grid  nximbers  of  the  areas  covering  the  contaminated  areas  of  interest, 
determining  the  reduction  factors  fot  each  of  the  Individual  square 

areas  from  the  tabulated  values,  and  summing  the  resultant  individual 
contributions. 

D.  .Sfifltgat,^ 

Section  II  of  this  appendix  presents  the  explicit  formulation  of  the  square- 
grid  model  including  the  necessary  background  material  and  mathematical  expressions 
used.  Section  IV  illustrates  the  use  of  the  model  with  a  concrete  example.  An 
overlay  grid  is  constructed  and  applied  to  a  hypothetical  municipal  area. 


II.  EXPLICn’  FOEimATlON  OF  THE  SQUARE-GRID  MODE!. 


A.  Approach 

The  notation  used  throughout  this  research  memorandum  is  similar  to  that  used 
in  Reference  C-1.  The  planes  of  contamination  are  divided  into  SC-foot  square 
sections  (determined  by  a  grid  overlay)  with  the  detector  centered  vertically  over 
any  desired  grid  (chosen  in  each  instance  to  be  the  origin  of  a  rectangular 
coordinate  system). 

The  detector  response  due  to  the  individual  square  lattice  grids  is  then 
computed  using  the  L(X,a))  functions,*  provided  the  effective  attenuation  distance 
X  and  the  solid  angle  oj  are  known.  It  should  be  noted  here,  however,  that  the 
determination  of  the  detector  responses  requires  special  precautions  since  Che 
L(X,aj)  curves  of  Spencer  are  all  computed  for  a  contaminated  plane  which  is 
symmetric  about  the  detector  location.  This  is  not  the  case  with  square-grids 
in  Che  square  lattice  model.  Consider  the  situation  depicted  in  Figure  C-1.  If 
one  desires  to  know  the  contribution  due  to  the  square-grid  in  the  upper  right 
hand  side  of  the  Figure,  it  is  first  necessary  to  compute  the  L(X,u))  function  for 
the  solid  angle  subtended  by  the  entire  square  annular  region  enclosed  by  the 
broken  lines  in  Figure  C-1.  Then  an  appropriate  weighting  must  be  applied  to 
estimate  the  detector  response  from  the  individual  square-grid.  This  weighting  is 
not  as  simple  to  compute  as  in  the  case  of  the  circular  model  wnore  a  simple  propor¬ 
tion  was  sufficient.  The  weighting  scheme  devised  fir  the  square-grid  model  is 

^it  was  found  that  these  functions  could  be  approximated  over  their  range  of 
greatest  utility  by  simple  analytical  sxpresaions  (Saction  111).  Such 
expressions  can  be  evaluated  with  the  aid  of  a  desk  calculator  and  standard 
nathomatical  tables. 
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Figure  C-1 

Concentric  Squaret  Configurations 
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discussed  in  Section  tx-b. 


1.  Assumptions 

The  basic  assumptions  underlying  the  square-grid  ^del  are: 

(a)  Each  50  foot  square  is  contaminated  uniformly; 

(b)  The  gamma  radiation  dose  received  by  the  detector  from  a  shielded 
contaminated  square-grid  is  due  to  both  direct  gazmna  radiation 
passing  through  the  barrier  and  air-scattered  gamma  radiation; 

(c)  The  detector  response  from  disjoint  sources  is  equal  to  the  sum 
of  the  contributions  from  the  individual  soui'es;  and 

(d)  It  is  verified  that  the  orientation  of  a  square  '•^rid  with  respect 
to  the  detector  location  need  not  be  considered.  The  distance 
from  the  center  of  an  individual  30  foot  square  to  the  detector 
and  the  distance  between  lattice  points  provide  the  geometric 
characteristics  reqdred  in  the  computation  of  detector  response. 
This  property  is  of  considerable  practical  importance  to  the 
construction  of  simple  graphs  and  tables  for  determining  the 
solid  angle  subtended  at  the  detector  by  any  lattice  square  as  a 
function  of  distance. 


2.  Geometry 

The  square  lattice  or  square-grid  model  is  constructed  In  a  manner 
similar  to  Che  circular  model  of  Appendix  B.  The  two  basic  differences 
between  the  circular  model  and  the  square-grid  model  lie  in: 

(a)  The  procedures  for  computing  the  solid  angles;  and 

(b)  Ttie  techniques  for  creating  the  oositlon-deaendent  parameters 


used  to  compute  the  individual  contributions  of  square-grids.  The 
lattice  uses  a  unit-square  annulus  (Figure  C-1);  azimuthal  sector 
contributions  are  used  in  the  circular  model. 

Consider  the  labeling  scheme  for  the  unit  squares  of  the  lattice  Indicated 

in  Figure  C-2.  In  this  Figure,  the  first  number  of  a  pair  (1)  in  a  unit  square 

identifies  the  unit-square  annulus.  The  second  niuaber  (j)  specifies  the 

position  of  the  unit  square  in  the  square  annulus  by  indexing  counter  clockwise. 

In  addition,  there  is  a  third  number  associated  with  each  unit  square,  h. ,, 

^  J 

which  specifies  the  height  of  a  unit  square  above  the  base  reference  plane. 

If  we  let  a  be  the  width  of  a  unit  annulus,  then  the  solid  angle 
fractions  subtended  by  the  n^^  square  annular  region  is  given  by 

a  oi  -  ^  -  T  (Vi- 

The  solid  angle  fraction  subtended  by  a  unit  square  is  given  by 


’3  ?•  ■  "ij) 

-  -14  2-  J'ij  •  ‘  •  '’lj>  J 


where 


T(x,y,z  -  h) 


I  ^  -1 

T  tan 


ihLixl _ 

\(z  -  h)  If  +  y^  ■♦•  (*  -  j 


(C-4) 


(C-5) 


**1 


(C-6) 


C-IO 


3,10 

3,9 

3,8 

D 

3,5 

B 

3,11 

D 

2,6 

2,5 

B 

2,3 

3,3 

■ 

3,12 

2,8 

B 

B 

B 

2,2 

3.2 

3,13 

2.9 

B 

B 

B 

B 

B 

3,14 

2,10 

B 

B 

1,8 

2,16 

3, 2  A 

3,15 

2,11 

2,12 

1 

2,13 

2,14 

2,15 

3,23 

3,16 

3,17 

3,18 

3,19 

3,20 

3,21 

3,22 

Figure  C-2 

Square-Grid  Labeling  Scheme 
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and  the  coordinates  (X.  , 

^  J 

1 ,  j  square , 


and  y^j) 


are  the  coordinates  of  the  center  of  the 

/■* 


The  above  equations  are  essentially  those  presented  in  Spencer's 
Monograph  42  (Reference  C**l)  on  pages  68  and  69. 

Differences  in  the  L^(X»w)  functions  tabulated  in  Spencer's  Monograph 

are  appropriate  to  give  the  reduction  factor  (^)  ior  the  square  annular 

o 

regions  described  above.  Of  course,  the  proper  solid  angles  must  be  substi¬ 
tuted  in  their  arguments. 


In  this  section,  barrier  shielding  of  contributing  source  squares  is  not 
treated.  This  is  treated  separately  In  the  next  section  of  this  report. 
Barriers  require  the  use  of  the  L  (X,u))  function  in  addition  to  the  L  (X,cu) 

A  C 

function.  These,  functions  may  be  approximated  by  simple  functions  of  the  form 


L  =  AtD°'exp  (BX) 


(C-7) 


where  A  is  the  Intercept  value  and  B  is  the  slope  of  the  curve  (See 
Section  III). 

In  order  to  show  specifically  the  scheme  for  computing  the  detector 
response,  Figure  C-3  is  used. 

In  terms  of  this  notation  and  the  discussion  of  this  section,  it  follows 
D 

chat  the  ^  of  this  typical  square  uniform  annulus  region  is  given  by 
*■0 

Since  there  are  8n  unit  squares  in  the  square  annulus  (n  ^  1),  the 
contribution  of  a  typical  unit  square  of  the  square  annular  region  may  be 
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Figure  C-3 
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written  in  the  form 


D 

o 


u 


nU 

8nu) 

n 


L  (X  ,  CD  ) 
c'  n,j,  n' 


L  (X 
c  ' 


h-l.jj^^n-pj 


L(X) 
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where 


CL- 


8n 


8n 


(C-10) 


03. 


The  factor  “  8r^  ^'n  *  “^n"  %-l  weighting  which 

corrects  for  the  effective  position  of  the  different  unit  squares  of  a  given 
square  annular  region. 


It  is  easily  seen  that  the  contribution  from  each  square  of  a  given 
square  annulus  n  is  weighted  according  to  the  solid  angle  fraction  it 
subtends.  The  total  geometry  reduction  factor  for  the  whole  plane  is  thus 


given  by 


=  2  s 

D  n  j  8nw 

o  n 


or 


r.  03  , 

iL  =  _  ^  -Ill 

D  2  ?  Act) 

o  n  j  n 


L  (X  .  03  )  -  L  (X 
c'  n,j,  n'  c'  ' 


L  (X  ,  CO  )  -  L  (X 
c'  n,j,  n'  c'^  ■ 


h-l,j,Vl^j 

n-1,  j/Vl^  j 


L(X)  (C-11) 


L(X)  .  (C-12) 


In  order  to  illustrate  the  manner  in  which  the  square-grid  model  is 
applied,  an  overlay  was  constructed  of  50  foot  square-grid  areas  (Joe  Figure 
C-2).  In  order  to  us  the  overlay  for  decontamination  analysis,  the  reduction 
factor,  D/Djj.  for  each  of  the  50  foot  squares  is  required.  Tlie  overlay  scaled 
to  apprc3priato  local  maps  could  then  be  quickly  and  effectively  utilized  in 
accessing  the  effects  of  large  scale  decontamination  in  a  particular  area  of 
ii.terest.  In  order  to  obtain  the  D/D^  values,  the  solid  angle  of  cacti 
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Individual  50  foot  square  must  first  be  determined.  An  appropriate  weight  is 
then  obtained  by  computing  the  total  solid  angle  of  the  annular  region  in 
which  a  particular  50  foot  square  is  located;  the  solid  angle  fraction  of  the 
individual  squares  may  then  be  determined. 

The  individual  calculations  of  the  solid  angles  are  quite^  length y .  In 
keeping  with  the  spirit  of  the  model,  an  approximate  but  highly  accurate 
procedure  for  computing  the  solid  angles  as  a  function  of  horizontal  distance 
is  derived  for  various  detector  heights.  Computations  of  the  solid  angles 
of  the  50  foot  square-grids  were  made  for  appropriately  spaced  on-axis  center 
distances.  (The  on-axls  case  has  desirable  symmetry  properties  which  simplify 
the  calculations.)  To  Illustrate  this  construction,  »  defector  height  of 
100  feet  was  used  (z  =  100  feet).  Calculations  of  the  solid  angles  subtended 
at  the  detector  by  the  50  foot  squares  were  made  for  each  of  the  successive 
on-axis  square-grids.  Calculations  for  several  off-axis  lattice  squares  were 
performed  (e.g.,  x  *  y  =  150  feet).  It  was  found  that  the  solid  angles  of  the 
50  foot  squares  depended  (to  a  high  degree  of  accuracy)  only  on  the  absolute 
distance  of  the  square  lattice  from  the  detector.  With  z  given,  the 
dependence  is  a  function  of  horizontal  distance.  Consequently,  the  curve  of 
Figure  C-4  was  drawn  showing  the  approximate  solid  angle  of  any  square  lattice 
for  three  detector  heights.  The  solid  angles  read  from  Figure  C-4  are  given 
in  Table  C-I  for  z  =  100  feet.  The  solid  angle  fractions  used  lor 

weighting  tne  L(X,a')  functions  are  also  given.  These  solid  angle  fractions 
are  computed  by  use  of  other  calculations  suamarlzed  in  Table  C-II.  ".'he 
reduction  factors  0/0^  (no  barriers)  for  a  particular  square  annular  region, 
n,  at  a  horizontal  distance  (X)  from  the  plane  of  contamination  is  then  given 
by 
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Solid  Axttle  Fraction 
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V'l  2 
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(Horizontal  Dietance  from  Detector) 


Figure  €*4 

Solid  Angle  Fraction  of  30  Foot  Square*  as  a  Function  of 
Horizontal  Distance  fros  Detector 
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TABLE  C-I 


Solid  Angles  and  Solid  Angle  Fractions  for 
the  50  Foot  Squares,  (z  «  100  Feet) 


HHRI 

HffiHI 

HH 

IMMU 

^^L(7.6)  /a.g(7.6,a)^) 

0.0 

0 

— 

0 

wm 

1.00 

.00632 

1,1 

50 

0 

msm 

HQI 

.00620 

1,2 

50 

50 

. 00478 

2,1 

100 

0 

.01416 

^5 

.00399 

2,2 

100 

.0122 

mSm 

.00346 

2,3 

100 

■■ 

.0079 

Bsl 

.00223 

3,1 

— 

0 

.00689 

.0547 

3,2 

WSm 

.00620 

.0493 

3,3 

.0046 

.0366 

3,4 

mSm 

mSM 

.0032 

.0254 

.00104 

4,1 

200 

0 

.0036 

.0421 

.00148 

4,2 

200 

50 

.0034 

.0397 

.00139 

4,3 

200 

100 

.0029 

.0339 

.00118 

4,4 

200 

150 

.0022 

.0257 

. 00090 

4.5 

200 

200 

.0016 

.0187 

.00066 

5,1 

250 

0 

.0022 

.0366 

.00107 

5,2 

250 

50 

.0021 

.0350 

.00103 

5,3 

250 

100 

.0018 

.030 

.00087 

5,4 

250 

150 

.0014 

.  0233 

.00068 

5,5 

250 

200 

.0010 

.0167 

.00049 

5,6 

250 

250 

.00079 

.010 

.00029 

6,1 

300 

0 

.00127 

.0288 

.00051 

6,2 

300 

50 

.00120 

.0272 

.00048 

6,3 

300 

100 

.00105 

.0238 

. 00042 

6.^ 

300 

150 

.0008 

.0182 

.000338 

6,5 

300 

200 

. 00063 

.0143 

.000255 

6.6 

300 

2  50 

. 00045 

.0102 

.000180 

6,7 

300 

300 

. 00040 

.0091 

.000157 
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TABLE  C-II 

Solid  Angles  Subtended  at  the  Detector  (0.0.100') 
bv  the  Square  Annular  Regions  Centered  at  the  Origin 


louare  Annular  Reeion 

“n  . 

n  n-1 

n  -  0 

0.0375 

0.0375 

1 

0.2344 

0.1970 

2 

0.4i:5 

0.1830 

3 

0.5436 

0.1261 

4 

0.6292 

0.0856 

5 

0,6892 

0.06C1 

6 

0,7333 

0,0441 

These  reduction  factors  for  various  annular  regions  are  given  for  z  =  100'  in 
Table  C-III.  With  this  information,  the  contributi'-  (and  hence  the  effect  of 
de;;ontamination)  of  each  of  the  square-grid  areas  shown  in  the  overlay  diagram 
(Figure  C-2)  can  be  tabulated.  Consjquently ,  the  overlay  can  then  be  placed  on 
a  scaled  city  map  and  the  direct  contribution  of  any  area  can  be  «.?siiy 
determined  from  the  tabulated  contributions.  A  table  of  contributions  can  be 
prepared  for  any  desired  detector  height,  z.  For  large-scale  considerations, 
the  50  foot  square  areas  are  small  enough  to  be  linearly  divided  fo..-  partial 
contributions  with  little  error. 


TABLE  C-  n 


Solid  Anxics  Subteaded  at 

'  fu  Ijetcctor  (0.0. 100') 

iouare  Annular  fteaion  n 

/D^  -  L(7.6)a..^(7.b,w^; 

n  *  b 

U . OOb  3 

1 

0 . 0^44 

3 

0 ,  (>409 

*4 

').o35l 

5 

U  ,  U  .’''''' 

b 

0.0l7t> 
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Barriers 


An  approximate  treatment  of  barrier  attenuation  and  scattering  is  no^; 
considered.  Since  the  model  is  constrained  to  stay  within  the  framework  of 
the  results  of  Spencer  (Reference  C-l),  a  formulation  is  developed  in  terms  of 
Spencer's  L  (X,w)  and  L  (X,w)  functions  (Reference  C-l).  Here,  the  detector 

C  o 

is  assumed  to  be  located  in  an  open  area.  Howeve'r,  this  assumption  <^occ  not 
restrict  the  application  of  the  model  to  regions  in  which  the  detector  is  in 
an  open  area.  If  the  detector  is  in  a  structure,  it  will  be  necessary  to 
analyze  the  shielding  by  that  structure.  Procedures  boch  approximate  and 
more  detailed  are  available  for  this  analysis  (Reference  C-2). 

The  purpose  of  this  section  will  not  be  to  enter  into  the  treatment  of 
barrier  problems  in  general.  Rather,  the  purpose,  is  to  characterize  rhe 
treatment  of  barriers  in  a  manner  that  is  consistent  witli  the  approximate 
character  of  the  over-all  model  and  which  is  simple  enough  to  be  practical 
for  field  applications.  The  formal  treatment  of  barriers  in  the  spirit  of 
the  above  discussion  Is  first  considered. 

Consider  the  contribution  of  the  unit  square  of  the  n^^  square 

annular  region  to  the  detector  located  at  the  point  (x=0,  y=0,  z)  in  the 
absence  of  any  barriers; 


D  ,/D 
nj  o 


Jli 


L  (X  , 
c  air , 


.)  - 


L  (X  , 
c'  air 


■W>] 


L(X) 


(C-13) 


wher’  X  ,  is  the  effective  mass  thickness  of  air  between  the  detector  and 
air 

the  source  plane.  If  barriers  arc  present  which  intersect  the  solid  angle 

cone  subtcndid  at  the  detector  by  the  squsrc,  the  direct  and  the 

scattered  components  of  the  radiation  arc  affected.  With  intervening 

barriers  present  in  the  n,  n-l,  ....  square  annular  regions,  we  have, 

in  principle,  a  total  mass  thickness  X  defined  wir,"i  the  respect  to  the 
r  r  t  n.g 
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square  and  the  given  detector  location,  i.e., 


,.X  =  va  .X  ,  +  ra  ,  +  +  ?:  a  ,X  ,  + 

t  n,g  "  n,n,l  n,i  "  n,n-l,i  n-l,i  ...  "  n,n-m,i  n-m,i 


i  j  ^alr 


Here  the  a  .  are  weightings  which  ^ft€  most  appropriat^ilr^aei'ined  by  the 

n ,  n-  Gi ,  l 

fraction  or  the  solid  angle  (defined  with  respect  to^';^j)  int^^^sected  by  the 
barrier  at  the  square.  The  last  term  of  the  abpvd'^pression  treats 

approximately  the  unobstructed  fraction  of  the  solid  angle  which  is  filled  with 
air  It  si.juld  be  noted  chat  it.  practice  there  would  probably  never  be  more 
thar.  three  terms,  and,  in  general,  fewer,  in  these  sums  over  the  mass  thickness 
parameters  cf  the  barriers  located  in  any  square  annular  region. 


.th 


.fter  determining  the  mass  factor,  X  , ,  of  the  barriers  shielding  the 

^  » j 


n,j  square,  an  expression  fc.  the  relative  contribution, 

derived.  With  a  barrier  of  effective  mass  -hickness,  ^  symmetrically 

suriounding  the  detector  between  the  detector  and  the  n*"^  annular  region. 


0  /D  =  L(  X  .) 
n'  o  't  n.j'' 


c  r  n,J  n 


L  (,.X  ,,.c  ,) 

c 't  n, j  n-1'  I 


(C-14) 


The  fractional  part  (c  ,/A )  (D  /D  )  is  not  the  desired  contri- 
but ion  because  this  fraction: 

(a)  includes  i-f' j  *  scattering  contribution  not  actually  present; 
and 

(b)  does  not  include  the  scattering  arou.  the  barriers  of  ri.>.dir-tion 

criglnatlng  In  ,  , . 

n.  J 


Also,  differences  in  tho  L  (X  ,  )  fu'ictions  can^iot  be  used  in  their 

a'  n,J,  n' 
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usual  form  to  compute  this  contribution  because  they  Include  a  large  component 

of  scattering  into  cu  ,  from  the  infinite  plane  beyond  co  .  and  no  contribution 

n ,  J  n ,  j 

from  the  shielded  area  which  scatters  and  arrives  at  the  detector  outside  of 


However,  by  defining  a  function  in  terms  of  the  L^(X,co)  in  a  manner 
analogous  to  that  in  which  L^(X,u))  is  tefined  in  terms  of  L(X),  we  obtain  a 
function  approximating  the  contribution  from  a  given  area. 


et 


12^ 

L(X) 


This  function  eliminates  the  contribution  of  scattering  into  co  .  fror. 

n,j 

the  Infinite  portion  of  the  plane  not  Included  in  cu  ..  Differences  of  the 

n,  j 

L  ,  (i.e.,  AL  )  Include  a  scattering  contribution  into  the  solid  angle  Aa. 
ac  ac  n 

from  that  portion  of  the  source  plane  corresponding  to  the  area  defined  by 

the  solid  angle  a  , .  The  AL  functions  also  include  the  direct  contribution 
^  n-1  ac 

from  the  plane  into  Aco  .  However,  the  AL  does  not  include  a  contribution 

n  ac 

to  the  detector  by  radiation  originating  in  od  -  od  ,  which  is  scattered  out 

n  n“  1 

V-  -  “^n-l  subsequently  scattered  back  into  the  detector. 

These  are  approximately  compensating  contributi.^'’s  and  we  will  use  this 
function  to  compute  the  approximate  contribution  from  barrier  obstructed 
source  squares. Clearly,  for  situations  where  most  of  the  radiation  from 
any  given  source  square  is  unobstructed  by  barriers,  this  approrimatlon  will 
be  adequate.  When  most  of  the  source  planes  are  obstructed  by  very  large 
barriers,  then  (3.2)  would  be  a  better  approximation.  Tlius,  we  cake  for  the 
approximate  contribution  of  the  square  source  area  with  effective  mass 

bfere  exact  expressions  could  be  developed  for  this  contribution;  however,  it 
was  felt  rhat  this  effort  would  not  be  Justified  until  after  the  potentiali¬ 
ties  of  a  square-grid  point  aourc*^  njodel  sd  been  explored. 
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thickness,  the  expression 


D_  ,  o)  , 


(C-15) 


We  will  show  below  that  the  values  of  (C-15)  may  be  presented 

graphically  for  a  large  range  of  X.  This  represents  a  conservative  approxi¬ 

mation  when  most  of  the  source  squares  are  not  obstructed  by  barriers. 

The  detector  response  due  to  a  given  contaminated  30  foot  square  area  is 
given  by  the  function 

f(a)L(X)AL^^ 

where  ^  is  defined  by  the  relation 
ac  ^ 

with 

"■.c  -  Tw  -  L  (•;;-)  L.(j^,a.) 

Hic  quantity  f(a))  is  the  Appropriate  weighting  for  a  given  square  grid  area 
and  is  given  by  aVZlci.  The  function  L(X)^L  is  highly  dependent  on  the 

barrier  mass  thickness  X.  However,  for  a  given  annular  square, 

f C- )L(X)(3l.^^  is  dependent  on  the  solid  angle  of  the  Individual  square  only 
through  the  quantity  fvi.).  The  quantity  f(a)  varies  with  the  position  of  the 
50  foot  lattice  in  a  given  annular  square  1.  Consequently,  the  quantity 
L(X)AL^^  may  be  computid  as  a  function  of  X  for  any  given  square  annulus. 
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Then,  through  the  f(cD)  function,  the  reduction  factor  of  each  of  the  individual 

50  foot  squares  penetrating  a  barrier  of  any  thickness  S  may  be  easily 

obtained.  The  function  L(X)i^L  as  a  function  of  X  was  computed  for  rings 

&c 

0  through  6,  and  the  results  tabulated  in  Table  C-IV.  Figure  C-5  also 

presents  these  data  graphically.  It  is  noted  that  the  variation  of 

L(X)AL  as  a  function  of  X  is  for  all  practical  purposes  a  single, 
etc 

exponential.  (The  variations  from  the  straight  line  on  the  semi- logarithmic 
plot  are  well  within  the  accuracy  with  which  Spencer's  curves  L(X)  and  L  (X,a)) 
can  be  read. 

Table  C-IV  gives  the  values  of  f(co)  -  cjd/^cjd  for  each  of  the  50  foot  square- 
grids  through  ring  six.  Although  these  weighting  functions  for  any  given 
annular  sector,  i,  vary  from  one  square-grid  to  another,  an  average  value 
of  f(co)  for  each  ring  can  be  obtained  through  the  function  AcoSi,  where  8i  is 
the  number  of  50  foot  square-grids  in  ring  number  i.  The  accuracy  of  using 
the  average  solid  angle  for  every  square-grid  in  an  annular  ring  is  certainly 
well  within  the  general  accurcy  obtainable  for  the  reduction  factor  in  any 
problem  in  which  complex  barriers  are  present.  Thus,  we  assign  this  average 
value  to  each  50  foot  square-grid  in  a  particular  riug.  With  this  assignment, 
the  reduction  factor  of  any  square-grid  through  a  barrier  of  arbitrary  thickness 
may  be  computed  and  presented  in  graphical  form.  Table  C-V  presents  the 
necessary  data  from  which  Figure  C-6  is  drawn.  Figure  C-6  thus  represents 
the  reduction  factor  D/D^  the  individual  square-grid  areas  as  a  function 
of  the  barrier  mass  thickness  X. 

Thus,  the  ease  of  readily  superimposing  the  values  of  D/D^  for  the 
various  square-grid  areas  for  the  case  of  air  with  no  barrier  Is  now  -arried 
completely  over  to  the  case  with  barrier.  One  needs  th.ee  ingredients:  the 
barrier  mass  thickness  X;  Che  square-grid  cverlay;  and  the  graphical  repre¬ 
sentation  of  Che  functions  f for  the  various  rings,  for  example, 
Figures  C-5  and  C-b. 
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TABLE  C-IV 


L(X)AL  as  a  Function  of  Barrier  Mass  Thickness  X 

_ 0C _  -  ■_  ■ 


L(X)AL^^  (X.co) 

_ aC _ 


X 

DSf 

Ring 

0 

Ring 

1 

Ring 

2 

Ring 

3  .  _ 

Ring 

4 

Ring 

^--5_  1 

Ring 

6 

24 

.00019 

.0054 

.0214 

.02788 

.02712 

.01740 

.00689 

72 

. 00014 

.0054 

.0109 

.00802 

.00399 

.00197 

.0012C 

144 

.00007 

.00164 

.00205 

. 000864 

.00031 

.00014 

216 

.00002 

.00037 

.000371 

.000105 

.00004 

.00003 

.00002 

288 

.00001 

.00009 

.OOOOul 

.00001 

<io“^ 

B 

TABLE  C-V 

Reduction  Factors  for  the  50  Foot  Squares  as 
a  Function  of  Barrier  Mass  Thickness  X 


D  ./D  -  -~L(X)  L  (X.  .  ) 

n  j  o  8n  ac 


DSf 

Ring 

0 

Ring 

1 

Ring 

2 

Ring 

Ring 

3 

Ring 

6 

24 

* 

7.  12(-6 

r~ 

1.33(-4 

2.44(*4 

1.47 (-4 

7.25(-5 

2.62  (-5 

5. 31  (-6 

72 

5.25(-h 

4.21(-5 

2.96(-6 

1.  l0(-6 

144 

2.52(-6 

4. 04  ('5 

2. 34 (-5 

4. 54 (-6 

8. 30 (-7 

2. 10(-7 

1.  19(-7 

216 

7.50(-b 

9.  l0(-6 

4.23(-6 

5.53(~7 

1.07(-7 

4.50(-b 

1.83(-8 

288 

3.75(-6 

2. 22  (-6 

6.9S(-7 

8. 62  (-8 

m 

-.o-' 

- - — - - 

•  lo'^ 

*7  . 1-  (-6  *  7.  I-  X 


L(X)/5I.^^  For  Ring  i 


* 


Figure  C-5 

Variation  ot  L(X)AL  vs.  X 

AC 

For  Ring  i  •  3 
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Figure  C-6 

Reduction  Factor  D/D^  for  the 
Individual  Square  Grid  Areas 
as  a  Function  of  the  Barrier 
Mass  Thickness  X  for  Ring 
i  ”  3. 
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The  expression  (C-14)  which  gives  the  total  effective  mass  thickness  from 

the  detector  through  the  intervening  barriers  to  the  n,j^^  square-grid  area 

is  in  ,:,eneral  quite  difficult  to  evaluate.  However,  one  can  argue  on  physical 

grounds,  that,  for  the  regions  of  interest  in  the  particular  applications  of 

this  model,  the  computation  of  the  effective  mass  thickness  can  be  greatly 

simplified.  Indeed  it  is  possible  to  show  theoretically  that  the  expression 

(C-14)  is  approximately  related  in  a  rather  simple  manner  to  the  effective 

th 

line  distance  from  the  detector  to  the  n,j  ‘  square  area  for  large  distances 

and  barrier  detail  in  keeping  with  the  approximate  character  of  this  model. 

Then,  the  total  effective  mass  thickness  can  be  practically  estimated 

t  n,j 

as  follows: 

(a)  Compute  the  line  distance  from  the  detector  to  the 
center  of  the  n,j^^  square  radiation  source  area; 

(b)  From  considerations  of  the  barriers  (i.e.,  wall,  floor, 
roof,  etc.,  data)  collected  from  earlier  building 
surveys,  compute  the  effective  distance  through  which 
the  radiation  travels  in  the  barriers;  i.e.,  simply 
the  thickness  of  all  the  walls  through  which  the 
radiation  travels  (note:  the  obliqueness  cancels 

when  the  area  is  nrojected  to  the  vertical). 

(c)  Add  the  effective  thickness  for  each  type  of  barrier 
material  separately.  Then  use  the  formula  for  X 

in  conjunction  with  the  material's  mass  thickness 
table  (see  Reference  C-I,  p.  15)  for  each  material, 
then  sum  X's  for  all  the  materials  to  obtain 
the  total  barrier  mass  thickness. 
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III.  APPROXIMATE  ANALYTICAL  FIT  FOR  L  (X.ui)  AND  L  (X,a,)  FUNCTIONS 

’’  AC 


It  was  found  during  the  course  of  developing  the  circular  model  (Appendix  B) 

that  an  analytical  fit  over  a  rather  extensive  range  of  values  for  both  the 

L  (X  ,u))  and  L  (X,(jo)  curves  found  on  pages  41  and  42,  respectively,  of  Reference 
s  c 

C-1  is  possible  by  employing  an  expression  of  the  general  form 

a  AX 

D/D^  =  L(X,a;)  =  a.®  e^  ,  (C-15) 

where  s  is  the  slope  and  A  the  Intercept  for  selected  values  of 
X  and  uj-  The  slope  in  turn  is  obtained  by 

s  =  In  L(X,a)2)  -  In  L(X,uu^)  ,  (C'16) 

while  the  intercept  is  calculated  by  solving  Equation  C-15  for  A  as 
follows : 

A  .  In  LiX.g;)  -  s  In  g. 

X 

The  best  analytical  fit  for  the  L^CX.sl)  curves  was  achieved  by  obtaining 
the  L^(X,  )  values  given  on  page  41  of  Reference  C-2  corresponding  to  and 
U2  values  of  0.05  and  0.2  respectively  and  solving  for  the  Intercept  at  ,  ■  0.2 
tor  each  X  value  selected.  The  same  procedure  wts  followed  for  the  L^(X,..) 
curves  depicted  on  page  42  of  Reference  C-7.  The  slope  in  this  case  was 
calculated  for  ,  ■  0.01  and  ■  0.1  while  the  intercept  was  obtained  at 

..  "  0.1.  The  numerical  results  of  these  calculstions  are  presented  in  Table 
C-VI,  -nd  the  resultant  curves  are  shown  In  rlgures  C-7  and  C-8. 

It  is  apparent  when  comparing  th*.  L(X,  . )  curvet  developed  here  with  those 
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Figure  C*S 

CeoMtry  FdcCor  (L^(X,^))  For  Detector  Reeponec  to  Redietion 

Received  by  a  tetecCot  Froai  a  Limited  Cone  of  Directions 
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appearing  in  Reference  C-2  that  the  accuracy  of  the  L(X,od)  values  Deteriorate 
rapidly  due  to  the  lack  of  convergence  as  o)  0.  This  is  particularly  true 
for  those  values  of  U3  >  0.3  and  o)  >0.5,  associated  with  the  L  (X,u))  and 
L^(X  ,a0  functions  respectively.  Despite  the  accuracy  limitations  inherent 
with  this  analytical  approximation  of  L(X,cu)  functions,  there  is  no  evidence 
that  the  approach  would  not  produce  the  desired  results  when  appropriately 
used  iu  the  square-grid  model.  It  should  be  noted  that  there  are  no 
restrictions  associated  with  this  model  that  would  prohibit  the  exclusive  use 
of  tNe  exact  values  for  L  (X.oj)  and  L  (X,a))  presented  in  Reference  C-1. 

A  C 


TABLE  VI 

Calculated  Values  c  s  and  A  for 
Selected  Effective  Mass  Thickness  Values 


L^(X,cd) 

L^(X,cu) 

X 

s 

A 

X 

5 

A 

1 

1.125 

-0.725 

2.26 

1.115 

-0.396 

5 

1.055 

-0.0950 

7.53 

1.137 

-0.0531 

24 

1,055 

-0.0952 

22.6 

1.042 

-0.002.55 

72 

0,975 

+0.00445 

75.3 

0.916 

40.00726 

144 

0.957 

+0.004355 

216 

0.954 

+0.00361 

288 

0.930 

+0.00302 
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IV.  AN  FX'MPLE 


P 


Consider  the  situation  indicated  in  Figure  C-9  below.  We  wish  to  explicitly 

demonstra  ?  the  method  of  calculation  of  the  D  ./D  for  several  of  the  contri- 

n,  j  o 

buting  squares  of  the  grid.  Starting  with  the  application  of  the  overlay, 
discussed  in  Section  II-B  (see  Figure  C-3),  to  the  quadrant  of  the  idealized 
scaled  map  shown  in  Figure  C-S  for  all  of  the  contributions  j/D^  of  the  squares 
can  be  obtained  directly  from  the  appropriate  overlay  as  discussed  earlier.  Of 
course,  the  values  will  only  be  correct  when  the  radiating  source  square  contri¬ 
buting  to  the  detector  is  not  appreciably  obstructed  by  barriers.  In  particular, 
for  the  contribution  of  the  2,3  square  (for  a  detector  located  100'  above  the 

origin)  the  value  of  the  function  D*  -/D  =  (0.0366)L(7.6''[L  (7.6,0.418)  - 

4 ,  J  o  c 

(7 .6,0.234)]  =  0.00218  is  immediately  read  off  of  the  overlay. 


We  now  assume  that  an  equivalent  thickness  of  2 '  of  concrete  for  the  building 
blocking  the  contribution  of  square  (3,7)  has  been  determined.  Also,  we  assume 
that  an  equivalent  thickness  of  1'  of  concrete  for  the  saoxi  structure  blocking 
the  contribution  of  its  roof,  square  2,3  has  been  determined.  Then  the  effective 
mass  thicknesses  of  these  barriers  are  given  by  (see  Reference  C-1,  p.  15) 


=  2(z/A)pA  «  2(0.5)(144)(2')  »  288, 

X2  5  =  2 (0.5) (144) (1')  -  144 

The  approximate  0  ,(X)/D  for  the  source  squares  obstructed  by  barriers  cannot 

n,j  o 

be  read  from  the  overlay  because  the  li^^(X,u))  fuactiont  leading  to  this  contri¬ 
bution  are  functions  cf  the  various  effective  barrier  thicknesses  .X  present. 

An  important  but  eaey-to-apply  approach  has  been  discussed  for  obtaining  the 

terms  D  ,(X)/D  in  Section  II,  which  employs  a  table  and/or  graph  giving  the 
n,  j  0 
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(The  detector  location  is  at  (0,0,100')) 


Figure  C~9 

Idealised  Scaled  Map  (to  be  used  In  conjunction 
with  Overlay  -  see  Figure  C-3  end  Table  C-I) 


»  values  of  the  obstructed  contribution  of  each  source  square  D  /D  as  a  function 

n,j  0 

of  the  effective  mass  thickness  X  (see  Table  C-V  and  Figure  C~3).  Using  the 

graph  or  Table  C-V  to  evaluate  D.  ^(X  =  238)/D  for  the  example  of  this  section, 

*  j,  /  o 

we  find  «  8.617  x  10  In  terms  of  the  discussion  of  Section  I,  the 

formal  expressions  representing  the  approximate  barrier  obstructed  contributions 
of  the  source  squares  3,7  and  2,3  (Detector  location  (0,0,100*))  are  given  by 


3,7  o 


L(288)  (288 ,0)3)  -  1^^^88,032) 


and 


Do  c/D, 


2,5' “o  0)2 


-U), 


[ 

[ 


U144)  L^^(144,.v  -I<^^(144,o3p 


Using  the  information  on  the  overlay  relating  to  these  solid  angles  we  can 
write 


D3,7/Do-  (0.055) 


I 


L(288)  L  (288,0.544)  -  L  (288,0.418) 

*  ftC  AC 


and 


°2,5^°o  “ 


[ 


1(144)  1  1  (144,0.418)  -  1  (144,0.234) 

AC  SC 


Then  using  Spencer's  tables  (Reference  C-1,  pp.  37,42)  we  find: 


or 


-  (0.055)  I  (0.00024) (0.067) 

-  (0.077)  I  (0.0063)(0.349)| 
o 

D,  ,/D  •  8.86  X  10  ^ 

/  o 
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and 

D.  JD  •=  1.54  X  10‘^ 

o 

0£  course,  these  results  could  also  have  been  obtained  from  a  table  of  the 

function  L  as  a  function  of  (o  with  X  as  a  parameter. 

ac 


C-36 


V.  SUMMARY  AND  CONCLUSIONS 

It  has  been  shown  that  a  square*grld  overlay  technique  may  be  used  to  analyze 
the  radiation  dose  contribution  to  a  detector  from  individual  planes  of 
contamination.  Once  these  contributions  are  known,  the  effectiveness  of 
decontaminating  any  plane  may  easily  be  assessed.  Data  required  for  appli¬ 
cation  of  tne  square-grid  model  are: 

(a)  A  grid  overlay  with  each  square  labeled  in  a  coordinate 
scheme  (a  50  foot  grid  size  is  proposed  with  each  square 
labeled  (l,j)  where  i  is  the  square  aitnulus  containing 
the  square  (i,j)  and  j  is  the  position  of  the  square 
(i,j)  in  the  annulus.  The  detector  is  vertically  above 
the  center  of  square  (0,0)). 

(b)  A  map  of  the  area  under  analysis  which  is  of  the  same 
scale  as  the  grid  overlay. 

(c)  Data  on  the  effective  mass  thicknesses  of  structures  in 
the  area. 

(d)  A  set  of  tables  giving  the  contribution  from  each  square 
in  the  grid  to  the  detector  located  a  given  distance 
above  the  plane  of  the  grid.  The  tables  are  constructed 
as  a  function  of  the  detector  height  above  the  plane. 

(e)  A  set  of  curves  giving  the  effects  of  barrier  shielding  (as 
a  function  of  X)  on  the  detector  response  for  each  square 
in  the  grid. 

Ihe  required  tables  and  curves  are  derived  for  a  detector  height  of  100  feet 
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and  a  square-grid  length  c.f  3C  feet.  It  Is  shown  for  this  case  that  the  grid  may  be 
approximated  by  a  function  of  the  dificuice  of  the  center  of  the  square  from  the 
center  of  the  (0,0)  square.  Changing  the  orientation  of  the  square  with  respect 
to  the  detector  produces  only  a  very  small  perturbation  in  the  value  of  the  solid 
angle. 

It  ia  concluded  that  when  detailed  maps  of  a  standard  scale  are  available 
for  an  area  (such  as  Sunbom  maps)  a  square-grid  overlay  technique  can  offer  a 
practical  tool  for  analysis  of  decontamination  effectiveness.  It  is  recommended 
that  data  be  derived  for  a  useful  range  of  detector  heights  (say  3',  10',  30',  as 
well  as  the  present  100')  and  that  the  method  be  applied  to  some  real  situations. 

Further  work  in  barrier  effects  is  required.  A  more  refined  approach  to  the 
analysis  of  barrier  effects  should  be  sought. 
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V 

e 

"i 

L(X) 


U) 

X. 


th 

"  area  In  square  feet  of  the  1^^  contaminated  area 

■  fraction  of  the  vertical  barrier  voided  by  apertures 

■  percent  of  the  vertical  barrier  voided  by  apertures 

>  an  area  of  contamination  such  that  the  line  from  the 
detector  to  the  center  point  of  the  area  is  perpendicu¬ 
lar  to  the  plane  of  the  area 

«  horizontal  barrier  shielding  factor  for  centered  areas 

>  horizontal  barrier  shielding  factor  for  off-center  area 
»  vertical  barrier  shielding  factor  (no  horizontal 

barvier  between  source  and  vertical  barrier) 

■  vertical  barrier  shielding  factor  with  horizontal 
barr  's  between  the  source  and  vertical  barrier 

“  doee  con -.ribution  per  unit  area 

«  doee  cont’^ibution  from  centered  area 

»  eccentricity  ratio,  length/width 

«=  height  of  the  i  area  above  reference  plane 

•»  total  dose  from  an  infinite  plane  source  of  fallout 

■  horizontal  distance  from  the  detector  to  the  center  of 
the  1  area 

<■  solid  angle  fraction 

*  horizontal-shield  equivalent  mass  thickness  in  psf 
<■  vertical-shield  equivalent  ma^s  thickness  in  pcf 
height  of  detector  above  reference  plane 
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A  Point “Source  Model  for  Approximating 
Gamma  Ray  Intensity  at  a  Single  Detector  Location 

I.  INTRODUCTION 


A.  Purpose 

This  appendix  presents  a  procedure  for  estimating  the  relative  radiation  dose 
contributions  at  a  detector  from  nultiple  planes  of  contamination.  The  two 
objectives  in  developing  this  model  are: 

1.  to  provide  accurate  approximations  to  the  relative  dose  contributions 
by  many  different  planes  of  fallout  contamination,  and 

2.  to  reduce  to  a  minimum  the  charts  and  tables,  maps,  equipment,  and  time 
required  to  calculate  the  contributions  to  a  given  detector  location. 

B.  Background 

Two  other  models  with  similar  objectives  have  been  developed  under  this 
contract.  Conceptually,  the  "circular  model"  (Appendix  B)  is  simple,  and  it 
is  easy  to  use,  but  accurate  approximations  are  difficult  to  make  when  the  con¬ 
taminated  planes  under  analysi.3  are  not  segments  of  circular  annuli.  The  "square 
grid  model"  (Appendix  C)  can  be  used  to  approximate  accurately  the  more  usual 
geometric  arrangement  of  rectangular  planes,  but  application  requires  a  suitable 
scale  map  for  use  with  overlays  of  the  grid.  This  point-source  model  presents 
a  scheme  for  approximating  contaminated  planes  by  the  appropriately  located 
point-sources  of  radiation.  The  point-source  model  is  formulated  such  that 
accuracy  limits  are  known.  Scale  maps  are  not  required,  but  only  certain  dimen¬ 
sions;  only  a  small  number  of  concise  charts  are  required  for  its  application. 

C.  The  Model 

In  this  model  the  detector  and  source  point  geometry  are  presented  in  cylin- 
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drlcal  geometry.  Since  the  detector  ie  isotropic,  the  geometry  reduces  to  two 
distances.  These  are  [h-sj,  the  vertical  separation  of  the  detector  from  the 
source  plane,  and  r  ,  the  horizontal  separation  of  the  detector  axis  from  the 
point  representing  a  given  contaminated  area.  Here  h  is  the  height  of  the  source 
plane,  and  z  is  the  height  of  the  detector,  both  above  a  reference  plane.  Figure  D-1 
illustrates  this  geometry. 


Figure  D-1 
Geometry  of  Model 

The  shape  of  the  contaminated  areas  to  be  represented  by  a  point  is  of  con¬ 
cern.  It  is  shown  that  for  the  area  centered  vertically  below  (or  above)  the 
detector,  a  circle  or  square  of  a  given  area  subtends  essentially  the  tame  solid 
angle  and  thus  yields  the  same  relative  dose.  Off-center  areas  are  all  assured  to 
be  squares.  Off-center  areas  which  are  not  square  should  be  approximated  as  closely 
as  possible  with  squares,  unless  the  maximum  dimension  meets  the  criterion  of 
Figure  D-11. 
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It  has  been  shovm  (Appendix  C)  that  the  orientation  of  a  square  area  within 
a  source  plane  does  not  appreciably  affect  the  solid  angle  subtended  at  a  detector. 
Thus,  to  a  good  approximation  the  relative  contribution  by  a  contamlnatec  area  to 
a  detector  is  a  function  only  of  the  vertical  height  |h-z|  separating  source  and 
detector,  the  horizontal  distance  r  from  the  center  of  the  area  to  the  vertical 
detector  axis,  the  size  of  the  area,  and  the  intervening  shielding. 

In  this  paper  the  relative  dose  rate,  or  reduction  factor  delivered  to  the 
detector  from  a  given  area  is  developed.  This  relative  dose  rate  for  area  i  is 
given  by  the  ratio  where  is  the  dose  rate  contributed  to  the  detector  by 

plane  i  and  is  the  total  relative  dose  rate  three  feet  above  an  infinite, 
uniform  plane  of  contamination.  The  reduction  factor  RF  is  then  given  by 

RF  =  1  D  /D 

all  planes  ^  ° 

The  relative  contribution  from  all  areas  beyond  those  analyzed  individually  may 
be  obtained  by  subtracting  the  "centered  area"  relative  contribution.  Figure  D-3, 
from  the  infinite  plane  contribution  L(lh-z|)  given  in  the  right  margin  of 
Figure  D-3. 

A  major  consideration  in  the  development  of  this  model  was  the  establishment 
of  the  maximum  area  of  a  square  which  could  be  accurately  represented  by  a  point- 
source  at  the  center  of  the  square.  A  criterion  was  devised  and  maximum  areas 
were  calculated  as  a  function  of  horizontal  separation  for  several  vertical 
separation  distances  (|h-z|  ■  3',  10',  20',  50',  and  100').  The  criterion 
establishes  that  for  a  square  area  leas  than  or  equal  to  the  maxinum  area,  the 
relative  contribution  calculated  on  an  area  basis  from  Spencer's  curves  (Reference 
D-1)  are  within  107.  of  the  contribution  calculated  from  the  equation 

C(A^.lh^-zI,r^'  .  C^(|h^.z|,r^)A^ 
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v;here  (|h^-z|,r^)  is  the  contribution  per  unit  area  as  a  function  of  |h^-z|  and 
r^  and  where  is  the  contaminated  area  in  square  feet. 

The  factors  for  the  shielding  effects  of  barriers  between  the  source  and  the 
detector  are  derived  from  Spencer's  curves  and  other  established  sources  (Ref¬ 
erences  D-3,  D-^ ,  f-j.  D-5) .  They  are  presented  in  the  standard  form  of  multiplicitive 
adjustments  of  the  unshielded  dose  from  a  source. 
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11. 


DEVELOPfENT  OF  THE  POINI -SOURCE  MODEL 


The  analysis  of  the  effects  of  decontamination  of  multiple  contaminated 
planes  surrounding  a  detector  usually  presents  unique  requirements.  This  model 
seeks  to  provide  a  concise,  rapid  method  of  calculating  accurately  the  fraction  of 
the  total  fallout  dose  rate  at  a  detector  contributed  by  each  contaminated  area. 

In  decontamination  analyses,  the  detector  positions  of  interest  are  usually  unpro¬ 
tected  or  are  in  low  protection  factor  structures.  This  is  because  decontamination 
is  most  effective  in  hastening  the  resumption  of  operations  in  poorly  protected 
areas  (Reference  D-2) .  Shielding  is  treated  in  an  approximate  manner  consistent 
with  its  relative  importance  in  the  model  and  with  the  requirement  of  conciseness. 
Techniques  of  the  "Equivalent  Building  Method"  (Reference  D-5)  and  the  "Engineering 
Manual"  (Reference  D-3)  are  combined. 

A .  Areas  Centered  Opposite  Detector 

The  area  immediately  below  or  above  a  detector  is  always  of  interest  in 
determining  radiation  doses.  This  area  may  be  a  major  contributor  to  the  dose  in 
the  absence  of  significant  shielding.  In  any  case,  the  difference  between  a 
centered  area  dose  and  the  infinite  plane  dose  (assuming  uniform,  unshielded 
contributions)  is  required  as  the  contribution  from  beyond  the  areas  which  are 
analyzed  specifically. 

1 .  Centered ‘Circles  and  Squares 

In  this  model  we  wish  to  approximate  the  dose  from  centered  areas  by  a 
function  of  only  area  and  detector  height.  If  we  arc  to  eliminate  shape  as 
a  parameter,  the  most  obvious  concern  Is  whether  or  not  centered  circles  and 
squares  of  equal  areas  contribute  equally  to  the  detector.  Assuming  that 
the  contribution  from  a  .entered  circle  or  square  is  a  function  only  of 
detector  heigi.t  and  o*^  the  solid  angle  .subtended  by  r!\e  area,  the  desired 
apprcximatlon  of  an  equivalence  criteria  between  squares  and  circle.s  is  tliac 


equal  solid  angles  are  subtended  by  equal  areas  at  equal  relative  detector 
heights.  Comparisons  of  solid  angle  fractions  subtended  are  given  in  Figure 
D-2.  It  is  seen  from  this  figure  that  for  the  range  of  areas  and  detector 
heights  used  in  this  model,  th'?  solid  angle  fraction  subtended  by  centered 
circles  and  squares  of  equal  area  are  equal  within  the  limits  of  accuracy  which 
can  be  used  in  reading  dose  rate  contributions  from  graphs.  Thus,  the  dose 
rate  contribution  from  centered  circular  areas  is  taken  as  being  equal  to  the 
dose  rate  contribution  from  centered  squares  of  equal  areas,  as  is  the 
standard  practice  (References  D-1  and  D-3) . 

In  Figure  D-3  the  relative  dose  rate  contributions  C,  from  a  contaminated 

d 

circular  or  square  area  centered  above  or  below  the  detector  is  presented  as 
a  function  of  area  for  several  detector  heights  jh^^-zl.  The.se  values  were 
derived  in  exactly  the  same  manner  as  were  Spencer's  (X,(jo)  functions  where 
X  is  the  equivalent  mass  thickness  between  the  source  and  detector  and  ui 


is  the  solid  angle  fraction  subtended  by  the  source  (Reference  D-1) .  That 


is,  letting  L(|h^-z|)  be  the  total  dose  rate  received  at  a  detector  located 
|h^-zj  feet  above  an  infinite  plane  of  fallout  contamination  of  such  intensity 
that  L(3')  *  1,  then 


where 


L((h^-zi)  -  L 


(D-2) 


A. 

I 


Thus,  Cj(jh^-zj,  k^)  corresponds  to  Spencer's  L(X)L^(X,cj)  function.  The 
L( I h^-z|)  values  given  in  the  margin  of  Figure  D-3  are  the  relative  dose  rates 
from  an  infinite  area  of  fallout  contamination,  A  <■  «,  at  a  detector  at 
height  |hj^-zj  . 
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Figure  D>2 

Coaperison  of  Solid  Angles  for  Centered  Areas 
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2.  Centered  Rectangles 

Areas  encountered  in  practical  situations  are  most  often  rectangular. 

We  shall  thus  investigate  the  maximum  eccentricity  of  a  centered  rectangular  area 
such  that  the  dose  rate  contribution  is  within  a  giver  margin  of  the  dose 
rate  from  a  centered  square  (or  circle)  of  the  same  arec. 

A  centered  rectangular  area  (i.e.,  L/w  ^  1)  will  always  crjntrlbute  less 
dose  to  the  detector  than  a  centered  square  of  the  same  area.  Thus  by  approx¬ 
imating  a  centered  rectangle  by  a  centeied  square  of  the  same  area,  one  will 
overestimate  the  dose  contribution  from  the  area  and  consequently  will  over¬ 
estimate  the  effectiveness  of  decontaminating  the  area.  This  means,  of  course, 
that  the  dose  rate  overestimates  should  be  kept  to  a  small  percent  of  the  doss 
rate  contribution  from  the  area.  This  overestimate  of  the  relative  dose  rate 
of  decontaminating  the  centered  area  will  mean  that  tho  reduction  of  decon¬ 
taminating  other  areas  will  be  underestimated.  We  shall  establish  bounds  on 
the  error  introduced  sufficiently  stringent  to  insure  acceptable  accuracy 
in  the  application  of  the  modal. 

Considering  the  precision  of  other  factors  in  the  total  decentamination 
model,  an  error  limit  of  ten  percent  seemed  acceptable.  The  criterion  estab¬ 
lished  below  assures  that  the  error  in  unshielded  dose  contribution  introduced 
by  this  approximation  is  less  than  ten  percent. 

Consider  the  two  geometries  of  Figure  D-4. 


figure  n-4 

Approxlaetely  Equivalent  Areas  •  Circular  Segaents 
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We  cee  that 


AA  ■ 

AJ  -  M/£ 

.fld  r'  -  ^  (f  +  *  (f )  ^ 

where 

A  ■  area  of  square  "center”  of  rectangle 
A  +  M  ®  urea  of  rectangle 
i  •  width  of  rectangle 
I  +  M  “  length  of  rectangle 
Thus,  the  eccentricity  e  of  the  rectangle  is 


e  =  i-LM 

i 


(D-3) 


We  have  shown  (Figure  D-2)  that  the  dose  from  centered  squares  and  circles  of 
equal  areas  are  identical  within  the  accuracy  of  the  graphical  data.  There¬ 
fore,  let  denote  the  area  A  in  a  centered  circular  or  square  configu¬ 
ration,  and  let  A^  denote  the  area  A  in  a  centered  rectangular  configuration. 

Thus,  the  dose  rate  from  the  area  on  the  left  in  ligure  D-4  is  (  [A  + 

2 

Letting  A^  •  n  r'  ,  the  dose  rate  from  the  area  on  the  right  is 

CjCIA  +  M  1^,)  -  C,(A^)  +  (CjCAp  -  C^(A^)1  .  (D.4) 

c 

Since  the  areas  are  equal,  we  make  the  approximation 


Cj  ([A  .  ([A  .  (D-5) 

Now  consider  the  way  in  which  the  off-center  M  is  appioxlmated  in  Equation 
0-5.  Figure  D-S  depicts  this  approximation,  and  it  is  clear  that  the  circular 
segment  approximatinn  of  Equation  D-5  uses  areas  farther  frexn  the  detector  and 
thus  gives  a  low  approximation  to  the  dose  rate  from  (A  +■  M)„. 
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Figure  D-5 

Approximately  Equivalent  Areas  -  Circular  Segment 


Figure  D-6  depicts  the  way  in  which  areas  are  equated  In  the  approximation 
Cj  ([A  ([A  +M]^).  (0-6) 


Figure  D-6 

Approxin«tely  Equivalent  Areas  -  Centered  Square 
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In  this  case  it  is  clear  that  the  centered  square  approximation,  Equation  n-6, 

uses  areas  closer  to  the  detector  and  thus  gives  a  liigh  approximation  to  the 

dose  rate  from  (A  +  M)-  •  Therefore, 

R 

([A  +M]^)  >  ([A  +My  >  ([A 

and 

([A  +M]p  -  (fA  +My  <  ([A  +M]p)  “  ([A  +  Ml^^)  • 

Cj  ([A  +  M]^)  Cj  ([A  +  M]^) 

The  dose  from  the  rectangular  area  ([A  +AA]g^)  is  then  within  ten  percent 

of  the  estimate  C,  ([A  +ziA]  )  if 
d  c 

Cj  ([A^-M]^)  -C,  ([A-^MlJ  <  0.1  .  <D-1) 

cJTTT+T^Tp 

Satisfaction  of  Equation  D-7  is  thus  taken  as  the  criterion  for  the  maximum 
eccentricity  of  a  centered  rectangle  of  area  A  +  M  to  be  approximated  by 
a  centered  square  of  area  A  by  using  the  dose  curves  in  Figure  0  3. 

The  locus  of  maximum  allowable  eccentricity  e  from  Equation  D-3,  as 
a  function  of  total  area,  A  +  is  derived  as  described  belcw  for  the 
detector  heights  of  3,  10,  20,  50,  and  100  feet.  To  obtain  a  point  on  the 
locus  for  a  given  height,  the  following  steps  are  taken: 

(a)  an  area  A  is  selected  for  the  center  square 

(b)  incremental  Al's  are  added  to  A  and  the  ratio  of  Equation  D-7 
is  computed  for  each  /A 

(c)  the  maximum  M  satisfying  Equation  D-7  is  found  graphically 

(d)  the  maximum  allowable  eccentricity  e  is  calculated  by 
Equation  0-3  for  the  area  A  +  M. 

The  e  values  thus  obtained  are  shown  as  functions  of  area  A  +  M 
max 

In  I'lgurc  D-7.  The  form  of  these  curves  Is  anticipated  hy  noting  that 

Lira  I  C  ,  ( f  A  +  ,  A  ]  1  “  C ,  (  (A  +  M  ]  )  i  -0 
\  d  c  d  es  I  ^ 

A  *^  « 
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Area  (ft 


Eccentricity  e  »  LengthA^ldth 


Figure  D-? 

Maxioum  Eccentricity  Ratloe  For  Centered  Rectengaler  Area* 
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since  as  the  center  square  area  Increases,  the  dose  rate  from  A  approaches 
the  dose  rate  from  the  entire  plane.  And 


Lim  I  C^([A  +  M]^)  -  ([/  AA^  ,g)  |  »  0 

[a  +  AA]“^0 


since  a  zero  area  contributes  zero  >se  te  regardless  of  the  "eccentricity*'. 


B.  Off-center  Areas 

As  implied  by  the  title  "Point-Sou  Me  ;1",  we  wish  to  establish  a  scheme 
to  represent  areas  of  contamination  by  ]  int  :  ources.  There  was  no  requirement 
to  introduce  point  sources  for  centered  areas  in  the  previous  section.  For  off- 
center  areas  we  shall  establish  criteria  by  which  off-center  areas  may  be  represented 
by  an  equivalent  point  source  at  the  center  of  tbe  area.  Actually  we  will  not  use 
point  source  data.  To  approximate  the  dose  contribution  from  the  area,  we  will 
multiply  the  area  by  the  contribution  per  unit  area  at  the  center  of  the  area. 

This  approach  requires  two  sets  of  data  for  its  application.  First,  the  dose  rate 
contribution  per  unit  area  as  a  function  of  horizontal  distance  from  the  detector 
is  required  for  each  detector  height.  Next,  a  set  of  criteria  for  the  maximum 
area  which  can  be  treated  in  this  manner  must  be  established. 

1 ,  Dose  Rate  Contribution  Per  Unit  Area 

The  dose  rate  contribution  per  unit  area  as  a  function  of  horizontal 
distance  from  the  detector  was  calculated  for  each  detector  height  using 
the  equation 


Cj  ( |h-z| ,r) 


C.(|h-z!,r  -  C .  (}h-z|.r-^) 


2xrAr 


(D-8) 


where 


h-z 


r 

( jh-z ! .r) 


vertical  aeparation  of  detector  and  source  area 
horizontal  separation  of  detector  and  center  of  source  area 
contribution  per  unit  area  at  center  ot  source 
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Dom  Eat«  ContrlbutloD  Par  Unit  Area 


(|h-z|,r)  “  contribution  from  a  centered  area  of  radius  r  as 

defined  in  Equation  D-2,  (the  change  of  argument  from 
A  to  r  is  only  for  clarity  in  Equation  D-8) 

2irrAr  =  area  between  circles  of  radius  and  . 

The  function  C^(|h-z|,r)  is  presented  in  Figure  D-8  fur  the  selected  |h-z| 
values  of  3,  10,  20,  50,  and  100  feet. 

The  differencing  in  Equation  B-8  to  obtain  these  values  must  be  accomplished 
before  the  maximum  allowable  Ar  is  established.  To  assure  the  validity  of 
these  values,  the  increments  are  chosen  well  belov;  values  which  previous  work 
(Appendix  C)  and  judgement  indicated  would  be  the  maximum  values  of  Ar. 

These  judgements  are  confirmed  by  the  maximum  values  of  Ar  established  below. 

2 .  Maximum  Area  to  be  Represented  by  a  Point 

Since  there  was  obviously  some  limit  to  the  validity  of  the  point  source 
approximation,  the  following  limiting  arguments  were  devised.  Consider  the 
convex  function  shown  in  Figure  D-9. 


Figure  D-*) 

Approximations  on  a  Convex  Function 
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Because  of  the  convexity  of  the  function,  the  tangent  to  f(p)  at  any  point 

touches  the  curve  only  at  that  point,  and  the  chord  from  (f(r)  at  to 

At 

(f(r)  at  r4y-)  intersects  the  curve  only  at  the  two  end  points.  Thus, 

■  f(«f  )  ■ 

f(p)dp  <  - 

and  ~ 

'  f(r-^)  +  f(r-»^)  ■ 

f(p)dp  -  f(r)  Ar  <  - s - 5 - ^  AT  -  f(r)  at  (D-9) 

which  means  that  the  error  introduced  in  approximating  the  integral  by  f(r)Ar 
is  less  than  or  equal  to  the  difference  between  f(r)  and  the  average  of 
f(r-|^)  and  f(r+|^)  . 

The  dose  per  unit  area  functions,  C^(|h>z|,r)  are  convex  over  most  of 

their  range.  The  log-log  scale  of  Figure  D-8  obscures  this  fact,  but  the 

general  form  of  each  of  these  curves  is  like  that  in  Figure  D-10  when 

shown  on  linear  scales.  The  value  r  ,  such  that  r  >  r  ,  implies  f(r)  is 

o  o 

convex,  is  less  than  ten  feet  for  all  of  the  values  of  jh-zj  <  100  feet.  The 
inequality  (D-9)  thus  is  used  in  establishing  the  criterion  for  maximum  area 
to  be  represented  by  a  point.  Only  the  convex  portions  of  the  ({h~z|,r) 
functions  are  used. 
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As  stated  earlier,  in  accordance  with  our  objective  of  making  a  slightly 
conservative  estimate  of  decontamination  effectiveness,  we  establish  the 
criterion 


It  has  been  shown  that  the  dose  contribution  from  an  off-center  square 

area  is  insensitive  to  the  orientation.  Thus  we  set  the  maximum  area  to 

be  represented  by  a  point  located  at  it.s  center  to  be  a  square  of  area 
2 

A  (Ar  )  ,  where  Ar  is  the  roaximum  .Ar  satisfying  the  condition  of 
max  max  max  ’  ® 

Equation  D-10,  and  where  r  is  the  horleontal  distance  to  the  center  of  the 
square  area.  For  areas  which  are  not  square,  the  criterion  should  bo  that 
the  maximum  dimension  is  less  than  or  equal  to  At  The  locus  of  maximum 
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square  area  to  be  represented  by  a  point  at  the  center  is  given  in  Figure  D-11 
for  each  of  the  detector  heights  3,  10,  20,  50,  and  100  feet.  If  the  area 
does  not  meet  the  criterion,  it  should  be  subdivided  into  areas  which  do. 


C.  Shielding 

In  the  previous  sections  we  have  considered  only  unshielded  detectors.  In 
this  section  we  shall  obtain  conventional  multiplicative  factors  to  account  for 
shielding.  The  cases  considered  are: 

(a)  horizontal  shielding  between  the  detector  and  a  centered  area 
of  contamination; 

(b)  horizontal  shielding  between  the  detector  and  an  off-center  area 
of  contamination;  and 

(c)  vertical  shielding  between  the  detector  and  an  off-center  area  of 
contamination. 

1.  Centered  Areas  and  Hcrizontal  Barriers 

The  horizontal  barrier  shielding  factor  for  centered  areas  is  given 
by  the  equation 


(X^,!h-zI.A) 


(X^,|h-z!.A)  (X^.uO 

(0,Ih-z|,A)  "  (0,w) 


(D-11) 


where 


•  shield  equivalent  mass  thickness  in  psf 
|h-z|  <•  vertical  separation  of  detector  and  source  in  feet 

A  *  area  in  square  feet 

Cl)  ■  solid  angle  fraction  subtended  by  area  A  at 
height  jh-zj 


c„  (X^. 


overhead  contribution  for  the  given  X  and  w 
(Reference  D-3)  ° 
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(9|)  1  U0|3VJ»C9S  X*)uo*I^OH 
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Figure  D‘11 

Maxlanan  Square  Area  To  Be  Represented  by  a  Point 
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Horlscntal  F.arrler  Factors  For  Cantered  Areas 


Since  ^  function  onl>  of  the  solid  angle  subtended,  separate 

graphs  ate  not  given  for  each  |h-z|  in  Figure  D-12.  Rather,  a  nomogram  is  given 
in  Figure  D-12  for  converting  areas  at  different  |h-z|  values  to  equivalent  areas 
at  ih-z|  ■  10  feet.  Essentially  this  is  the  same  approach  as  that  taken  in 
Reference  D-S. 

The  (X^,]h-zj,A)  values  of  Figure  D-12  are  used  in  the  following  manner 
to  obtain  the  shielded  dose  contribution  from  a  centered  area: 

(a)  find  Cj  (jh-zj,A)  using  Figure  D-18 

(b)  in  the  Figure  D-12  nomogram  locate  the  given  area  value  on  the  A 
scale  (on  the  left) 

(c)  draw  a  straight  line  from  this  point  A  •'through  the  given  value  of 
{h-z[  on  the  middle  scale  and  continue  to  the  A'' scale  on  the  right 
side 

(d)  from  the  intercept  of  this  line  with  the  A'’  scale,  follow  a  horizontal 

line  to  tne  curve  for  the  correct  X  (if  X  falls  between  two  curves, 

0  0 

an  interpolation  should  be  made) 

(e)  read  the  (X^,jh-z|,A)  value  at  the  bottom 

(f)  multiply  ([h-z|,A)  by  (X^,[h-z|,A)  to  obtain  (X^, jh-z | ,A) . 

This  would  be  the  contribution  from  the  centered  area  Areas  also 
shielded  by  vertical  barriers  should  be  treated  as  off-center  areas. 


2 .  Off-center  Areas  and  Horizontal  Barriers 

The  horizontal  barrier  shielding  factor  for  off-center  areas 
by  the  equation 


(X^.|h-r!,n(r^^  )V;  -  (X^.  |  h-s  |  ,  iT(r-^)‘) 

'  '  . .  y"*""  . . . . .  . . y—  ■' 

(0.|h-*j,^(r4^)^)  -  (O.lh-z:  ,«(r.^«)S 


is  given 


(l)-12) 


where  th*f  terras  are  the  same  as  those  defined  for  Equation  D-ll,  except  that  the 
are?,  has  been  given  by  a  function  of  the  horizontal  separation  of  aource  and 
detector. 
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The  values  of  (XQ,jh-z|,r)  are  given  in  Figure  D-13.  Again,  separate 
curves  are  not  given  for  each  |h-z|.  The  (X^, [h-z j ,A)  values  are  functions 
only  of  the  solid  angle  subtended  by  the  Area  A  which  is  defined  by  r  in 
Equation  D-12.  The  nomogram  thus  can  be  used  for  converting  radii  at  different 
|h-z|  values  to  equivalent  radii  at  |h-z|  «  10  feet.  The  values  of  Figure 
D-13  are  used  in  analogous  manner  to  those  of  Figure  D-12.  That  is,  to  obtain 
the  shielding  dose  contribution  from  an  area: 

(a)  find  C^(|h-zi,r)  from  Figure  D-8  and  multiply  the  value  by  the  area 
(a  square  area  satisfying  the  criteria  of  Figure  D-11) 

(b)  in  the  Figure  D-13  nomograms,  locate  the  given  horizontal  separation 
distance  r  on  the  r  scale  (on  the  left) 

(c)  draw  a  straight  line  from  this  point  r  through  the  given  value 
of  |h-z{  on  the  middle  scale  and  continue  to  the  r'  scale  on  the 
right  side 

(e)  read  the  (XQ,|h-zj,r)  value  at  the  bottom 

(f)  multiply  (|h-z|,r)  j  A  by  B^^^  (XQ,|h-z|,rX  (If  there  are  no 
vertical  barriers,  this  is  the  contribution  from  that  area.) 

3.  Vertical  Barriers 

As  a  cultiplicative  barrier  factor  for  vertical  barriers,  we  again  need 
the  ratio  of  the  contribution  with  the  barrier  present  to  the  contribution 
without  the  barrier.  First,  the  case  for  which  no  horizontal  barriers  exist 
between  the  source  and  the  vertical  barrier  is  considered.  The  contributions 
through  the  walls  are  obtained  from  values  in  "The  Lquivalcut  Building  Method" 
(Reference  D-5) .  To  obtain  these  contributions,  we  find  the  protecvlon  factors 
(PF)  with  no  roof  contribution  (infinite  root  weight)  and  we  note  that  the  PP 
is  the  reciprocal  of  the  relative  contribution  through  the  wall.  Thus  we 
obtein 
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Horizontal  Barrier  Factors  for  Off-Center  Areas 


(D-13) 


C^(X^,A*)  PF  (co.O.A*) 

(0,A*)  PF  («,X^A*) 


where 


B  (X  )  «  vertical  barrier  factor  with  no  horizontal  shield 

V  V 

between  the  source  and  tl  -t  vertical  barrier 

X^  =  equivalent  mass  thickness,  in  psf,  of  the  vertical  barrier 

*  *  2 
A  =  (2r  y 

r  -  horizontal  distance  from  detector  to  the  vertical  barrier  , 

k 

The  PF  values  in  Reference  D-5  are  given  for  areas  A  in  values  ranging  from 

100  to  100,000  square  feet.  Since  these  numbers  include  wall  scattering  as 

well  as  direct  attenuation,  it  is  conceivable  tliat  the  barrier  factor  (X^) 

is  also  a  function  of  the  distance  r  ,  However,  the  ratio  of  equation  0-13, 

* 


when  obtained  for  the  complete  range  of  areas  A  given  in  Reference  D-5,  shews 
nc  appreciable  difference  between  the  values  obtained  for  tie  different  areas. 
That  is,  the  variance  introduced  by  reading  the  graphs  of  Reference  D-5  was 
greater  than  the  differences  in  the  values  for  different  areas.  Thus,  a 
single  curve  of  B^  (X^)  is  given  in  Figure  D-14. 

If  a  horizontal  barrier  exists  between  the  source  and  the  vertical 
barrier,  the  curve  derived  from  Reference  D-3  for  this  configuration  will  be 
used.  This  curve  is  reproduced  in  Figure  D-1'  as  ' 

If  there  are  apertures  in  a  vertical  barrier,  corrections  to  the  mass 
thickness  of  the  solid  portion  should  be  made  before  using  Figure  n-14.  In 
the  configuration  with  an  intervening  horizontal  barrier,  the  mass  thickness 
of  the  solid  wall  should  be  multiplied  by  the  fraction  of  the  wall  not 


voided  by  apertures.  Thus, 
X^  «  (l-.Ap)  X 
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Figure  D-14 

Vertical  Barrier  Shielding  Factors 
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where 

Ap  *  fraction  of  wall  voided  by  apertures 

X  *  mass  thickness  of  solid  portion  of  wall 

In  the  configuration  with  no  intervening  horizontal  barriers,  correction 

curves  for  apertures  are  presented  in  Reference  D-5.  However,  a  consideration 

of  their  basis  indicates  that  the  correction  curves  of  Reference  D-5  deviate 

somewhat  from  the  form  which  they  should  possess  (see  Chapter  III)  .  Because 

of  this  deviation,  and  since  the  corrections  do  not  differ  greatly  for 
* 

different  areas  A  ,  a  single  ’’consensus"  set  of  correction  curves  will  be 
used  in  this  model.  Figure  D-15  presents  the  aperture  corrections  to  be 
applied  in  this  configuration  with  no  intervening  horizontal  barrier.  To 
use  these  correction  curves,  compute  the  precent  of  the  vertical  barrier 
occupied  by  apertures  %  Ap,  and  determine  the  mass  thickness  of  the  solid 
portion  of  the  barrier  X  .  For  these  values  of  %  Ap  and  X  ,  read  fron 
Figure  D-15  the  value  of  X^  to  be  used  in  Figure  D-14  to  obtain  (X^) . 
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Correction  to  Vertical  Barriers  for  Apertures 


III. 


EVALUATION  OF  EQUIVALENT  BUILDING  METHOD  APERTURE  CORRECTIONS 


Shelter  Design  and  Analysis  Volume  II  presents  the  Equivalent  Building 
Method  of  Fallout  Radiation  Shielding  Analysis  Design.  The  Equivalent  Building 
Method  is  based  on  the  approximation  of  reducing  a  complex  shielding  situation  to 
an  equivalent  simple,  solid-wall,  single-story  structure  problem.  For  above  the 
ground,  the  basic  structure  assumed  in  TR-20,  Volume  II  is  a  solid-wall,  single- 
stc  y,  square  building.  The  wall  height  is  13  feet  with  the  detector  3  feet  above 
the  floor.  The  aperture  sill  heights  are  3  feet,  the  detector  height.  Basic 
structure  areas  of  100,  1,000,  10,000,  and  100,000  square  feet  are  used.  For 
each  of  these  areas  a  family  of  curves  which  relate  the  protection  factor  to  wall 
mass  thickness,  with  roof  mass  thickness  as  a  parameter,  is  given.  When  a  structure 
varies  from  the  basic  case,  an  equivalent  mass  thickness  is  found  from  given  curves 
such  that  a  standard  structure  with  this  wall  weight  has  the  same  protection  factor 
as  the  structure  under  study.  In  particular,  corrections  for  apertures  are  made 
using  Figure  10  of  TR-20,  Volume  II,  reproduced  here  as  Figure  D-16. 

The  purpose  of  this  appendix  is  to  examine  the  peculiarities  exhibited  by 
these  curves  at  mass  thicknesses  less  than  100  psf.  In  figure  D-16,  Ap  is  the 

I 

fraction  of  the  wall  taken  by  apertures;  so  let  Ap  be  the  fraction  of  wall  above 

the  detector  taken  by  apertures,  and  let  all  apertures  be  above  detector  height. 

Only  the  ground  contribution  C  is  affected  by  apertures  in  the  exterior  walls; 

8 

therefore,  equivalent  wall  mass  thickness  is  found  by  equating  ground  contributions 


Cj(Ap,  A.)  -  n;) 


Since 


(O.X^) 

-Hi — 


<  0 


(Figures  1  through  4  of  References  D-3  and  D-5) 


ArEKTUkCS-AtOVEGtOUNO 


(Source:  Figure  10  of  TR-20,  Volume  II  -  Reference  (1)) 


C^(Ap  ,  xp  <  C^(0,  X^)  <s«o  x;  >  Xg  . 


(D-15) 


Using  the  method  of  the  Engineering  Manual  (TR-20,  Volume  I),  we  have 

Cg(4.  •*  <l-*p)  G.  li-s„<Xg)l 

+  id-A^)  G^(u^)  +  G^(ug)]  S^(X^)E  I 


+  Ap  C  (oj  ) 
a  u 


(D 


Given  Ap  <1,  and  Ap  <1,  then  by  equation  D-16 


II  I 

C^(Ap  ,  X^)  <  C^(Ap,X^)  is  equivalent  to 


{-  *p  0^  (V  ri-s„i  -  4'  G,  s„E  I  +  a;;  c^)  < 

Ap  G^  (uj  (l-S^j  -  4  G^  („_^)  S^E  I  B^  +  Ap  G^  (u^) 


or 


itR 


Ap 


^  ic  (to  )  [1-B  +SB]-G  (u))S  Eb1> 
I  a  u  '■  e  w  e^  s  u  w  *  /  ” 


G  (u)  )  [1  -  B  +  S  B  ]  -  G  (w  )  S  EB 
a  u  '■  e  w  e^  s  u  w  e 


(D-17) 


If 


Ap  ■  0  ,  then  by  equation  D-17,  C  (Ap  ,  X^)  <  C  (0»X^)  and  only  if 


0  >  O.  (-u) 


-  B  +  S  8  -  G  (,j  )  S  EB 

C  W  cj  •  u  W  € 


V“u>  '  •  ,  X  po. 


(D-18) 


-16) 
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t  ft 

Now,  for  any  Ap  <  Ap  ,  satisfaction  of  condition  D-18  implies  satisfaction  of 

It  t 

condition  D-17.  Thus,  with  Ap  >  Ap,  if  all  apertures  are  above  the  detector,  then 

X"  >  X  <=Bs>  C  (Ap  ,  X  )  <  C  (0,X  )  <■—■>  C  (Ap,  X  )  <  C  (Ap,X  ) 

e  -  e  g  '  ’  e  -  g  e  g  e  ~  g  e' 

<«>  C  (0,  X")  <  C  (0,  X')  <3=s?>  X"  >  X'  . 

g  ’  e  —  g  ’  e  e  —  e 

I 

For  a  given  area,  then  X^  plotted  to  X^  with  Ap  as  a  parameter  has  the  form  of 
Figure  D-17. 


Figure  D-17 

Correct  Form  for  Aperture  Corrections 


The  crossover  point,  ,  is  given  by  the  equality  of  condition  D-18.  In 
condition  D-18,  for  a  given  eccentricity  ratio,  the  right  side  is  independent  of 
area  and  the  left  side  is  independent  of  X^.  The  crossover  point  may  then  be 
found  as  in  Figvire  D-18. 
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Deterainatlon  of  Points  of  Equality  of  Condition  D>18 


IV.  AN  EXAMPLE  ANALYSIS 


The  application  of  the  point  source  model  is  illustrated  by  a  simple,  but 
fairly  realistic  example.  Consider  the  decontamination  situation  of  Figure  D'19. 

In  this  example  the  detector  is  located  t^  cee  feet  above  an  open  area  (a  parking 
lot  or  park,  say)  which  covers  a  city  block.  Area  I  includes  the  parking  lot  and 
the  adjoining  streets.  Area  II  is  the  flat  roof  of  a  building  one  block  long,  75 
feet  wide,  and  15  feet  high.  Area  III  is  another  parking  lot  separated  from  the 
detector  by  the  building. 

The  detailed  analysis  or  these  three  areas  is  shown  in  Table  D-I.  First,  the 
centered  area  is  analysed ;  the  analysis  is  very  simple  in  this  exan^le  since  the 
eccentricity  meets  the  criterion  and  there  .:.s  no  shielding.  If  this  area  does  not 
meet  the  eccentricity  criterion,  it  should  b^?  iubdivided  and  the  off-center  sub¬ 
divisions  should  be  treated  as  other  off-center  areas. 

Since  Areas  II  and  III  are  quite  elongated  and  the  areas  do  not  meet  the 
criterion  of  Figure  D-11,  they  are  subdivided  into  areas  which  are  nearly  square 
(in  this  case  exactly  square)  and  which  meet  the  area  criterion  in  Figure  D-11. 

The  data  for  Area  II  with  |h-z|  «  12  feet  is  interpolated  between  the  10  and  20 
foot  curves.  Note  that  the  horizontal  barrier  factor  is  off-scale  on  Figure  D-12 
and  thus  Is  carried  as  an  inequality.  Note  also  that  the  aperture  correction  for 
vertical  barriers  is  given  by  «  (1  -  )  X  for  Area  II  where  a  horizontal 

barrier  exists  between  the  vertical  barrier  and  the  source,  and  is  given  by 
Figure  D-15  for  Area  III  where  there  is  no  horizontal  barrier. 

In  the  application  of  the  point-source  model,  a  tabulation  sheet  such  as 
Table  D-I  will  be  found  very  useful. 
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Example  An.' 


.70(-l  =  5.70  X  10 


n-it 


TABU!  D-l 


Example  Analvsig  of  Coaflgurwtlon  In  Figure  D-19 


Horizontal 
Barrier  Xq 

(psf) 

Eo(XQ.lh-z|,A) 

or 

®oi^^o* 

Chart  12  or 

Chart  13 

*  ®0’ 
or 

AC^  X  Bq^ 

Vert.  Barrier 

X  (psf) 

%  Ap 

Aperture 

Correction 

V  ^  100^'^ 

or  Chart  15 

B*  (X  ) 

V  V 

or 

B  (X  ) 

V  v' 

Chart  14 

Contribution 

Cj  *  ®o  * 

AC,  X  Bn,  X  B 

1  01  V 

or  AC,  X  B 

1  V 

Total  Contribution 
of  Sub-Areaa  In 

Each  area  ■  C^^^ 

0 

1.0 

5.70(-l 

0 

0 

0 

1.0 

5.70  (-1 

5.7  (-1 

25 

<  .01 

<  2.0  (-5 

25 

20 

20 

.49 

<  1.0  (-5 

25 

<  .01 

<  2.5  (-5 

25 

20 

20 

.49 

<  1.25*(-5 

25 

<  .01 

<  2.5  (-5 

25 

20 

20 

.49 

<  1.25  (-5 

25 

<  .01 

<  2.0  (-5 

25 

20 

20 

.49 

<  1.0  (-5 

<  4.5  (-5 

0 

1.0 

2.0  (-3 

50 

20 

51 

.275 

5.5  (-4 

0 

1.0 

2.25(-3 

50 

20 

51 

.275 

6.2  (-4 

0 

1.0 

2.0  (-3 

50 

20 

51 

.275 

5.5  (-4 

1.72(-3 

V.  THE  EQUIVALENT  PIANES  METHOD 


The  Equivalent  Planes  Method  is  based  on  the  analyses  of  the  previous  sections. 
The  keys  to  the  method  are  two  limiting  criteria.  The  first  gives  limits  on  the 
length- to-width  ratio  (eccentricity)  of  an  area  centered  above  or  below  the  detector 
such  that  the  dose  contribution  is  within  ten  percent  of  the  approximating  dose 
from  a  centered  square  of  the  same  area,  nie  other  criterion  gives  limits  on  the 
area  of  an  off-center  square  such  that  the  dose  contribution  is  within  ten  percent 
of  the  approximation  given  by  the  product  of  the  area  times  the  dose  contribution 
per  unit  area  from  the  center  of  the  square. 

Shielding  factors  are  presented  for  above  ground  detectors  which  are  separated 
from  source  areas  by  relatively  light  shields.  These  data  should  be  most  applica¬ 
ble  to  decontamination  studies.  The  dose  contributions  calculated  are  suitable 
for  Inclusion  in  formulas  for  decontamination  analysis. 

A  concise  manuel,  giving  a  step-by-step  procedure  for  applying  the  Equivalent 
Planes  Method,  is  included  at  the  end  of  this  volume.  Several  data  tabulation 
sheets  to  be  used  with  the  manuel  are  also  included. 
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VI.  CCWCLUSIONS  AND  RECOMMENDAnONS 


It  is  concluded  that  the  point-source  model  can  be  a  useful  and  realistic 
method  of  analyzing  fallout  radiation  contributions  for  decontamination  studies. 

The  real  keys  to  the  method  are  the  limiting  criteria  given  in  Figures  D-7  and  D-11. 
The  first  criterion,  Figure  D-7,  gives  limits  on  the  length- to- width  ratio 
(eccentricity)  of  an  area  centered  above  or  below  the  detector  such  that  the  dose 
contribution  is  within  ten  percent  of  the  approximating  dose  from  a  centered  square 
of  the  same  area.  The  other  criterion.  Figure  D-11,  gives  limits  on  the  area  of 
an  off-center  square  such  that  the  dose  contribution  is  within  ten  percent  of  the 
approximation  given  by  the  product  of  the  urea  times  the  dose  per  unit  area  at  the 
center  of  the  square.  These  criteria  are  based  on  unshielded  areas.  In  future 
studies  it  would  be  desirable  to  investigate  these  criteria  for  shielded  areas. 

Shielding  factors  employed  in  this  model  are  obtained  from  well-known  sources 
(References  D-1,  D-2,  D-4  and  D-5).  Hie  approach  taken  is  based  on  the  premise  that 
above  ground  detectors  which  are  separated  from  source  areas  by  relatively  light 
shields  are  of  primary  interest  in  decontamination  studies.  If  applications 
Indicate  that  basement  areas  and  heavier  shields  are  of  more  interest,  the  model 
can  readily  be  extended.  Probably  it  will  be  found  necessary  to  generate  data 
for  Xq  values  between  zero  and  25  psf  to  fill  in  Figure  D-13  since  many  "curtain" 
walls  and  light  roofs  have  mass  thicknesses  in  this  range.  The  data  are  not  avail¬ 
able  from  Reference  D-3,  so  they  will  have  to  be  generated  directly  from  Spencer's 
data,  Reference  D-1. 
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Appendix  E 

A  FORTRAN  Program  for  Decontamination  Analysts 

I.  INTRODUCTION 


A.  General 

This  appendix  describes  a  computer  program  written  in  FORTRAli  lo  compute  the 
effectiveness  parain»»ters  used  to  analyze  municipal  decontamination.  The  program 
was  written  in  FORTRAN  64  to  be  used  on  large  scale  computers  such  as  the  CDC  3600. 
The  program  has  been  debugged  and  tested  and  has  been  used  under  OCD  work  unit  3233B, 
For  a  given  level  of  decontamination  the  program  is  capable  of  determining 
both  the  reductions  in  dose-rate  at  specified  detector  locations  and  the  re 
ductions  in  total  dose  for  persons  spending  prescribed  amounts  of  time  at 
specified  detector  locations.  Other  parameters  used  in  the  analysis  of  municipal 
decontamination,  such  as  the  PF's  at  the  detector  locations  and  the  equivalent 
PF's  associated  with  the  activities  (without  decontamination),  are  also  computed. 

B.  Contents 

Section  II  of  this  Appendix  is  a  description  of  the  computer  program  and 
includes: 

1.  a  list  of  the  equations  which  are  computed 

2.  the  data  input  formats 

3.  the  printed  output  formats 

•jif 

The  basic  equation  defining  the  equivalent  protection  factor  is: 

*  f,/P,  +  f./?^  +  .  ,  .  +  f  /P 
112  2  n  n 

where  is  the  fraction  of  time  spent  with  protection  P,. 
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4.  the  detailed  flow  charts 

5.  a  card  listing  of  the  program  with  a  set  of  sample  inputs 

6.  a  sample  output  of  the  program. 
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II.  THE  COMPUTER  PROGRAM 


A.  Introduction 


The  CDC  3600  FORTRAN  program  computes  the  values  of  the  following:  the  pro¬ 
tection  factor  (PF)  at  each  location,  the  equivalent  protection  factor  (EPF)  for 
each  activity  pattern,  the  fraction  (CF)  of  the  total  intensity  prior  to  decontami¬ 
nation  at  a  given  detector  due  to  a  particular  contaminated  plane,  the  intensity 
reduction  factor  (RN)  for  each  decontamination  strategy  at  each  detector  location, 
and  the  activity  dose  reduction  factor  (RN^)  for  each  combination  of  activity  and 
decon..  .injtion  strategy  at  all  detector  locations. 

Calculations  are  made  from  the  following  input  data:  Contributions  (C)  to 
intensity  at  each  detector  location  from  each  plane  of  contamination  the  fraction 
(f)  of  time  that  an  individual  spends  at  each  detector  location  performing  each 
particular  activity,  the  mass  reduction  factor  (f)  for  each  decontamination  strategy, 
and  the  key  (S)  to  associate  the  appropriate  value  of  F  with  each  composite  decon¬ 
tamination  strategy  for  each  plane. 

The  program  is  designed  to  accept  a  maximum  of  twelve  each:  detector  locations, 
planes  of  contamination,  activity  patterns,  composite  decontamination  strategies; 
and  twenty  mass  reduction  factor  values. 

IT.j  equations  used  to  compute  the  above  parameters  are  the  following: 

PF.  -  - ^ -  (E  n 

J  P 


EPF, 

k 


CF 


ij 


P 

}  C.  , 


<E-2) 


(K  3) 


E-3 


(E-4) 


(E-5) 


where: 

j  identifies  the  detector  location;  d,  the  number  of  detector  locations 
j  “  1,  2,  d  (d  mcximum  =  20) 

i  identifies  the  plane  of  contamination;  p,  the  number  of  planes  of  contamination 
1  =  1,  2,  p  (p  maximum  =  12) 

k  identifies  the  activity  pattern;  a,  the  number  of  activity  patterns 
k=l,  2,  ...,a(a  maximum  =  20) 

H  identifies  the  decontamination  strategy;  s,  the  num/ner  of  decontamination 
strategies.  £  =  1,  2,  s  (s  maximum  *  20) 

h  identifies  the  composite  decontamination  strategy;  c,  the  number  of  composite 
decontamination  strategies,  h  =  1,  2,  ...,  c  (c  maximum  »  100) 

C  .  identifies  the  contribution  to  intjusity  at  ,  detector  locatictn  from  ith 
i.l 

plane  of  contamination 

f  identifies  the  fraction  of  time  spen^  at  detector  location  j  perfcrmlng 
the  activity  denoted  r>y  k 

F_  identifies  the  mass  reduction  factor  (fraction  of  fallout  remaining)  for 
/ 

decontamination  strategy  c 

Identiiies  whif-h  strategy  £  Is  used  to  decontaminate  plane  I  for  composite 
strategy  h 

Data  Input  Fonnat 

The  first  input  card,  identified  by  a  2  in  card  column  one,  specifics  the 
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B. 


number  of  values  to  be  entered  for  each  variable  subscript  in  Integer  form  (XX)  as 
shown  below; 


Each  value  for  each  of  the  Input  variables  and  f^^  Is  entered  on  a  separate 

card  and  entered  in  sequence  by  row  (A,  A  A,  A,  ;  A,  A„  A. 

L,i  i,j  j.  ,n  z,  L  ZfZ 

^0  „ »  I  >  5 1  o  >  •  •  *  ^  — )  •  * 

z,n  m,z  in,J  in,n 


The  values  for  input  variables  F  (up  to  twenty)  are  entered  from  a  single 

H 

card,  The  first  value,  which  serves  also  as  the  identification  number,  in  card 
columns  1-4  must  always  be  1.00, 
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The  values  for  input  variable  (up  to  twelve  planes)  are  untered  by  row 
each  row  completed  in  one  card  which  is  numbered  in  sequence: 


Serial  ^3  “h,7  h,8  \,9  ^h,10  ^h,ll  ®h,]2 

'"t 


Data  cards  are  to  be  entered  in  the  following  sequence: 


c 


Printed  Output  Format 


D=XX,  P=XX,  A=XX,  S=XX,  C*XX 


Values  of  C(J,I)-- 


Detector 

Location 

1 

2 

PLANE  NUMBER 

3 

p 

No.  1 

.xxxxxxxx 

.XXXXXXXX 

.XXXXXXXX  .  .  . 

.xxxxxxxx 

No.  2 

.xxxxxxxx 

.xxxxxxxx 

.xxxxxxxx  .  .  . 

.xxxxxxxx 

No.  3 

.XXXX30DCX 

.xxxxxxxx 

.xxxxxxxx  .  .  . 

.xxxxxxxx 

No.  d 

• 

.xxxxxxxx 

• 

.xxxxxxxx 

•  • 

.xxxxxxxx  .  .  . 

.xxxxxxxx 

Values  of  FR(J,K)  — 


Detector 

Location 

1 

ACTIVITY 

2 

PATTERN  NUMBER 

3 

a 

No.  1 

X.XX 

X.XX 

X.XX  .  .  . 

X.XX 

No.  2 

x.xx 

X.XX 

X.XX  .  .  . 

X.XX 

No.  3 

x.xx 

X.XX 

X.XX  .  .  . 

x.xx 

No.  d 

« 

x.xx 

« 

X.XX 

•  • 

X.XX  .  .  . 

X.XX 

Values  of  S(H,I)-- 


X.XX 

X.XX 

X.XX 

Value 

■3  of 

S(H,I)-- 

Composite 

PLANE  NUMBER 

Strategy 

1 

2 

3 

No. 

1 

XX 

XX 

XX 

No. 

2 

XX 

XX 

XX 

No. 

3 

XX 

XX 

XX 

XX  XX  XX 


x.xx 


p 

XX 


XX 


XX 


E? 


No.  c 


XX 


PF(J) 


Detector 

Location 

No.  1  XXXX.XX 

No.  2  XXXX.XX 

No.  d  XXXX.XX 

EPF(K)  = 

Activity 

Pattern 

No.  1  XXXX.XX 

No.  2  XXXX.XX 

f  •  • 

No.  a  XXXX.XX 


CF(J,I)  = 

Detector 

Location 

1 

2 

PUNE 

3 

NUFBER 

4 

p 

No.  1 

X.XX 

X.XX 

X.XX 

X.XX 

.  .  .  X.XX 

No.  2 

X.XX 

X.XX 

X.XX 

X.XX 

.  .  .  X.XX 

No.  3 

X.XX 

X.XX 

X.XX 

X.XX 

.  .  .  X.XX 

• 

« 

* 

• 

•  * 

No.  d 

X.XX 

X.XX 

X.XX 

X.XX 

.  .  .  X.XX 

Check  by 

summing 

CF(I,J)  on  each  J. 

Should  equal  1. 

X.XX 

X.XX 

X.XX 

X.XX 

X.XX 

K-8 


RN(H,J)  = 

Compos  ife. 
Strategy 

1 

2 

DETECTOR  LOCATION  NUMBER 

3  A 

d 

No.  1 

x.xx 

x.xx 

X.XX 

x.xx  .  .  . 

X.XX 

No.  2 

x.xx 

x.xx 

x.xx 

X.XX  .  .  . 

X.XX 

No.  3 

x.xx 

x.xx 

X.XX 

x.xx  .  .  . 

X.XX 

• 

« 

• 

•  • 

• 

No.  c 

x.xx 

x.xx 

x.xx 

x.xx  .  .  . 

X.XX 

RNA(H,K)  = 

Composite 

Strategy 

1 

2 

ACTIVITY 

3 

PATTERN  NUMBER 

A 

a 

No.  1 

x.xx 

x.xx 

X.XX 

X.XX  .  .  . 

X.XX 

No.  2 

x.xx 

x.xx 

x.xx 

X.XX  .  .  . 

x.:oc 

No.  3 

x.xx 

x.xx 

X.XX 

X.XX  .... 

x.xx 

• 

• 

• 

•  • 

• 

No.  c 

x.xx 

x.xx 

X.XX 

X.XX  .  .  . 

x.xx 

PRINT 

ERROR  MESSAGE 


DOUBLE 

SPACE 


SET  £»l 


K-17 


TRIPLE  SPACE 


PRINT 

"CF(J,I)  =  '' 


DOUBLE  SPACE 


"DETECTOR. . . 
PLANE  NUMBER 


"LOCATION 


DOUBLE  SPACE 


CF(J,I)ON  EACH  J 
HOULD  EQUAL  IjjJ 


PRXKT  SCFXU) 


■ 

SET  j 

1 

RN  > 
jh 

0 

r 

SET  i- 

□ 

1 

£*S 

sn* 

RMjh  *  RN  jh 
+  CF  ij  F  £ 


G 


C-2i 


8-2J 


S-24 


E 


Card  Listing  of  Sample  Input  Deck 


'SEduENCE»339  ■"  . . . 

» J03» 12345« JRYAN 
*FTN«L«X 

PROGRAM  JRYAN 

DIMENSION  C(  1C»20)  *FR  (  20 .20 )  .F  ( 20  )  .LSI  1  0 . 1 0  )  .  PF  (  20 )  .EPf-  (20) 

_ 1  . 10> .RN( 100.20) .RNA(  100.20 .SCFl (20)  »NUM( 100) 

PRINT  lObo  '■  . 

iOCC  FORMAT! SOX 16H  GLOSSARY) 

PRINT  1001 

IPCI  FORMAJIfp!:!  I  .J)  0EN0T£S_  CONTRIBUTIONS  TO  INTENSITY  AT  JTH  DETEC'^;> 
1  LOCATION  FROM  ITH  PLANE  OF  CONTAMINATION*) 

_ PRINT  1002  _ _ _ 

1002  FORMAT  (*OFR(i<«  J)DENOTES' THE  FRACTION  OF  TIME  SPENT  AT  DETECTOR  J 
SERFOR^NG  THE  ACTIVITY  DENOTED  BY  K* ) 

PrTnT  1003 

1003  FORMAT! *OF(L)  DENOTES  THE  I^ASS  REDUCTION  FACTOR  (FRACTION  OF  FAllO 
3UT  REMAINING)  FOr’ DECONTAMINATiON  STRATEGY  L*) 

PR IN.I- 1 .9?.'?.  ,  _ _ _ _ 

13>.4  ‘f'o’rMAT  !*6£'(H.  i  I'CENOTES"  WHI’ch~STRATE’gy’'l  IS  USED  TO  DECONTAM I  l  l 
4  PLANE  I  FOR  COMPOSITE  STRATEGY  H* ) 

PRINT  lOOE 

1005  FpRMAT(*OPF!J)  "NOTES  THE  PROTECTION  FACTOR  AT  EACH  L'ETaCTOR  LC.  - 
5TI0N*) 

_ _  PBiNXJ.9j:'6 _ _  _ _ _ 

Toos' FORmXt(*CEPF(k:T  DErX'TCS  THE  EQu'l VAirENT'  PROTECfrCN’ factor  FOR  LAC’-- 

CJJ  y  I TY .  p J . 

PR~INT"i0'07 

1007  FORMAT !  *GCF_  ( Jj_I  )  DENOTES  THE. FRACTION  OF  THE  TOTAL  INTENSITY  RF  U - 
7  TO  d’e'cONTAm'inATION  AT”  A  GIVEN  DETECYoR*) 

_ PRINT  1008 _ _ _ _ _ _ _  .... 

ToOe  FORMAT!*  DUE  TO  A  PARTICULAR  CONTaVI  NAtXd  PLANE**')' 

P91N1  .IPP?  .  _ _ _ _ 

1009  FORMAT! VoRNlH.J)  DENOTES  THE  INTENSITY  REDUCTION  FACTOR  FOR  EAC:-^ 
BDECONTAMINATION  STRATEGY  AT  EACH  DETECTOR  LOCATION*) 

PRINT  1010 

1010  FORMAT (*ORNA(H.<)  DENOTES  THE  ACTIVITY  .DOSE.  REDUCT I  ON  FACTOR  FO- 
9ACH"'C0M8INATi0N  OF  ACTIVITY  AND  OECONTa'm  I  NAT  I  ON*  ) 

PRINT  1011 

1011  FORMAT!*  STRATEGY  At  EACH  uETECTCR  LOCATION** 

99  READ  1 . ICK. ID. IP. 1  A, IS. IC 

1  FORMAT! II .512) 

_ IF(£0F.6C )41 .42  _ _ _ _ 

42  IF! ICK-2»6l .43.61 

43  DO  10  I«l .IP 
00  10  J« I . 1 0 
READ  2. ICK.C! I «J) 

2  fo«maY( li .F9*ei 

_ IF!  IC<-3)62. 10.62 _ _ 

lo’cONTiNUE 

DC  2^  Xal  .  I  A 

00  20  J«r.  lD 

_ REf^  _ 

i  FORMAT!  1  I  .>4, *27 

_ IF(|CK»4)e3_.2aj6_3 . . . . 


CONT  I  Nut 

READ  4« (F(L) ♦Lsl » IS) 

4  F0RMAT(2CF4,a ) 

IF(FC 1 )-l , )64«30»64 
3C  00  40  M=1 . IC  ' 

.  REAP  5«  1C.K»NUM(M  )  «  (LS(M«  1  )  *  I  =  I1.IPJ _  _ 

!>( ICK-5J65«40«65  ~  "  . 

40  CONTINUE  __  _ 

PWINT  21  «  lb  .  1P«  I  a‘«  Is  V'fC  . . . . . . 

F0RMAT(4H1D  =«I3«3X.3MP  =.I3»3X.3MA  s:,13»jx«3HS  *.I3«3X.3HC  = «  I  j 
'PRINf  22  . . .  . .  '  '  . . 

22  FORMAT  ( /»OVALUES  OF  C(  1  «J) — tO _ _  _ _ _  _ 

‘prTnT  £20’ 

FORMAT (9H0DETECT0R.53X.22HP  LANE  NUMBER) 

PRINT  221  .n  .’1  =  1  .IP)  . 

£.:„l  FORMAT  (16H  LOCATION  .13.10111) 

DO  222  J=1.1D 

c2£  PRINT  1 1 . J. (C(  I . J)  . 1^1 . IP) _ _ 

n  f6rMAT(4h6nO. . iE.2X. 1 IFl r.8)  "  ~  . .  ^  - . 

PRINT  23 

23  FORMAT (/»OVALUES  OF  FRCJ.K) — *) 

PRINT  230 

23C  FORMAT  (9HCDETECTCR.42X.'43hA  C  T  1  V  1  T  Y  PA  T  T  E  R  N  N  U  c 

_  IE  R)  _ _ _  _ _  _ 

PR  1  NT  ~23  r.  {  1  .  1  =T'.~l  A ) 

231  FORMATdlH  LOCATION  .2016) 

DO  232" J=1 . ib  ” 

232  PRINT  l2«'Ji  «FR(I<«  J)  .Kal  ilA) 

12  FORMAT(4h6nO.'.  I5.3X«2bF6.2) 

PRJLNT  24  _  _ _ _ _ _  _ „  .  .  _ _ 

24  FORMAT  (  /#OVALuifj '"of  F(L  )--*•/)  "  “  '  •  “  . - . - 

PRINT  13. (F(L) .La  1 . IS) 

13  FORMAT( 12X.2bF6.2) 

PRINT  25  _ 

c5  FORMAT (/*OVALUES  OF  S(H*I)— *) 

_ PRINT  2f)£ _ _ _ _ _ _ 

2-.C  FORMAT(  lOHOCOMFOSlTi.52X.22HP  L  A  N  t"  N  U  M  B  E  R) . . 

PRINT  251  «  {  I  .  I_=  I  .  ’P)., 

2::i  F0RMAT(12H  STRATEGY  '  .1216) 

DO  252  M=  Ij  I C  . . 

232  PRINT  14,NUM"{M  >  ,  <L3{*M.  I  )  ,  I-'  1  .  IH# 

1^^  _ 

DOao'ja'r.'lD  -  --  -  —  . 

SUMCIJsC. 

00  so  l*i  ."I'p”'  ■  ’ 

uo  5UMC  i  J  *  3U  MC  I  J  C  ^  J  ._  _) 

6C  'PF(U).el.7SuMCVj  . 

■  •  FvR  I  NT  1 S _ ; _ _ _ _ 

15  rORMAT(»Jl.PF(  J.i  fc*  ) 

.  PRINT  ISO 

ITC  A0RMAT{9HbbeTEC'TbR./9H  LCCATiON)*  " 

PR  I  NT  1 51  .J  J  *_PFJ  J)  *jml  .1 0 

IT.  i  r  ORmAT  »4f.^Hb>''.'~d'.'F972)  '  ■  . 

_ _  30'  oO  xa  )  «  I A 
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SQ  =  0. 

.  .,.^,70  J=1*IP _ _  _ 

7C  SQ=SQ+FR(K, J)/PF( J) 

0Q__EE’f .y<>^U»/.SQ _  .  _  _ _ 

PRINT  16  . . 

16  FORMAT  (*0EPFtK)  s*) _ _ _  _ 

PRINT  160  “  . 

160  FORMAT  (  9H0 ACT  I  VI.TY  _*/8H  PATTERN  ) 

PRINT  161  * (KtEPFCO  ,Kal ♦ lA) 

161  F0RMAT(,4_HCN0.  »  I5*F9.2) 

DO  no  j’=r«iD  .  . . 

_ 5C1J  =  0 _ _  __  _ 

DO  100  1  =  1 «  IP 

ICO  _SCI^aSCl^+C(  i  »J_)_ _ _ 

'  do'  lio  ’i  =  i\"'ip . .  . . . 

no  CF(  uj)=c(  I « j)/scij  _  _ 

. 'print  17'  ''  . .  "  ■  '  . . . 

-  17  FORMAT(/»JCF(  jn  )  =  *) _ _ _ 

PRINT  220 

PR  1  NT  23 1.(1  ♦  1_=  1_«J  PJ  _ 

DO  r2d”'j='l  .Vd . . 

1 2  0  PR  I  NT_i2 .  J  «  (  CF  (  1 .0  )  ♦  I  =  1  .  l.PJ  . 

PrTnT  9 

_ 9  FORMAT  (»0  CHECK  BY  SUM'- ING  CFCJ.l)  ON  EACH  J . _  SHOULD,  EQUAL  1  /  / 

'  DO  130  Jal.ID  . . 

_ SCFl ( J)aQ _ 

DO  130  1  =  1  «  IP  ~  . . . 

1  3C_JC,F  1J_^)  f^CF  n  jn-CF(  I  .  JJ _ _ 

PRINT  13.  (SCF'h  J)  «  J=1  .  ID)  “  . . . . 

DO  140  Msl  «  IC _ _ 

DO  140  J=1 . ID 

_  RN^J.M)_=0.  _  _ _  _ 

”d6  u'd”  fsT.'iP 

L  =  LS{M.  1  )  _ _ _ _  __ 

14C  RnI J.M)'=RN{  j7m')’+CF(  I  .J)*F(L) 

_ PRINT  18 _ _ _  _ 

18  F0RMAT{/*0RN(H. J)  =*)  '  ~  . . 

PRINT_:30  _ _  __  _ 

l'e6~  FORMAT!  ro"HOc6MPOsVTE.3exT4*5HD  E  >  E  C^'t'O  R  LOCATION  N,. 

1M_B_E  _R)_ _  _ _ _ _ _ 

PR  I  NT'~254  ,'<  J  ♦  j=  r,  1 0 ) 

_ DO  253  Ms  1  .  1C _ _ 

’233  PRINT  12.NUM(M) , (RN( J.M) .J*l » ID) 

.  _PO  5*  1  '  _  -  .  _ _  _ 

"  do"  310  Mi  n  ic  '  '  . 

SUSQaQ.  ,,  _ 

'  suspi'd’i .  "  '  --  - 

DO  300  Ja 1 « ID _ 

SUSQ»SUSQ+FR(K.U)/PF{U)  -  -  - - . 

300  SySPaSUSP4-FR(>«;,U)/PF(  J)«RN(U.M) _ 

3 1  6~ RNA  ( K  Vm  )  ■  SUSP/SUSO  '  . 

PRINT  19 _ _ _ 

19  FORMAT  (/♦CRNAiM'.KI*  •{ 

PRINT  190  _ _ 
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150  FOR.vrATC  lOHOCON'.f-O?  ITE*36X«43HA  C  T  i  V  1  T  Y  p  A  T  T  i,  P 

13  E  _RJ  .  . . 

PR  I  NT  254 »  <  K » K=  i  «  I A  )* . . 

DO  191  M=1«IC 

1  <1  PRINT  12.  NUV!  ( M )  ♦  <  RNA  f  <  ♦  M )  « K=  1  .  lA  ) 

_ PRJ.NT_.2S _  _ _ _  -  . 

58  FORMATTlHi)  . 

GO  TO  99 

61  PRINT "71  . 

71  F0Rj^AT(*0  incorrect  D/^A_CAR_D— ^OT  VALUES  OF  SUBSCRIPTS.*) 
"  GO  to' 41  '  . . . .  . . . 

62  PRINT  72  „  _ _  __ 

'  72>ORMaT(*6'inc6rRE:CT  DATA“"cARD-'-NOf~v'ALUE  of  cfl.’j).*) 

GO  TO  41 

o3  PRINT  73 

73  FORMAT (*0 INCORRECT  DATA  CARD — NOT  VALUE  OF  FR{K.J).*> 

GO  TO  41 

.i4  J^^NT_7_‘» _ _ _ _ 

74  FORMAT  (  *0  fNCORRECf  DATA  1: A WD--Nof‘  VALUES  OF*  F  { Ci  •*  ) 

GO  TC  41 

o5  PRINT  75 

75  FORMAT (*0 INCORRECT  DATA  CARD — NOT  VALUES  OF  S(H«1/,*) 

GO  TO  41 

!.♦  12l$L  _ _ 

._54  FORMAT!  IIH  STRATEGY  .1216) 

4 1  CONT INUE 
END 

SCOPE 

•LOAD 

«RUN«3.2CCC.  1 _ _ _ _ _ 

4.040405060*6*  ■  '  . . 

■5*2S»9 

0.336 

-^.024 

3.702 

3.  nj _ _ _ _ _ _  _ _ 

3.031 

0.209 

0.005 

U£D6 

.  1 BZ _  _  _  _ _ 

>.'00  1 
-..071 
..j77 
•>.005 

.  ■- . v>  2 

4 .  0  6  __  _ 

.40.07""'  . . . 

-..20 

~.C.67  "  '  -  -  - 

"•  V  »  0  7 
4.;,  .17 
4  0 . 0. 9 


l-tf. 


I 


ft 


-*  0  »  M 
•40.02 

•4  .  j  .  w  3 

•^'.67 

•*-  •  13  _ 

•4  ^) .  0  5 
•^u  •  1  5 

-4  Vy  .  67 

-..,..19 

-+  Vy  .  C  8 
•■^.06 
-•0.67 

: .  ooc.02C.cJo.i-yo.  loc.  1 3 

■::0  1C1G5C102 
:  0?03'’3C401 
0301010401 
■..'04  0  1030601 
■,■'0304020  101 
-..^601010506 

*  I 


« 


* 
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F .  Sample  Output 


J  *5473  JHriV  1440? 

1  *4213  IU».  0:347 

1  S7411  TS4.  • 

*R03))4M 

•lORR 

IRBSLEO 

•iOMF 


1  630Sn  lOS,  C0423 
1  »774C  ILIOC.  00733 
1  07344  OSaiFlOC  00033 


1 

*1444 

USOfe^RO.^ 

00202 

1 

»7»13 

RSTAR 

00240 

1 

00141 

ION, 

02210 

1  *i>42  asoENrur  oooto 

1  S744S  STM,  000)* 


•■aMfi*iir.o  aoMMOM 

NON* 

txTiir  FoiMis 


0*074 

ifcOtMf 

1 

70334 

*19*2 

i*odi:t, 

1 

07340 

*73*1 

omsiosl. 

1 

47000 

fO-T. 

*7242 

OtOCMlIN 

1 

*2173 

oaontsT. 

*1*4* 

ISOfeMOOM 

•1 

41703 

OaUfeRSET 

07**1 

OSCHkCO 

1 

07414 

RSClOS^ 

077*4 

4tTj4'l, 

1 

40072 

Rust* 

001*0 

lan. 

1 

001*1 

Rci'Bur. 

07213 

.***C47, 

klECUtlON  ST4MIEU  <r  1417  >02 


1 

41044 

aOBEOTRY 

1  *2347 

THFNB) 

1 

57410 

TSn. 

irTirt 

STh, 

1 

42373 

10*. 

1  *2*70 

10$) 

1 

*23*4 

aonoagL, 

- T1**i* 

TO*;  ' 

i 

41043 

t»tT 

1  0791 3 

•  STAR 

1 

*7077 

MQISkLkO 

1  *011* 

•Luoe, 

1 

*0101 

IROruRN, 

«LLOelN 

1 

073*1 

aooirioc 

1  *»79» 

"IlS;" 

B>30 


ULOSSANY 


C<I,V»  nENJTESi  CUNtiJHuTlONS  TO  INTENSITY  AT  JIW  ntlfcCTOH  tOCAtfO;.  :  WDH  ITH  PLANE  Of  CONTANI^I^ 

f H(Nt JiOENOTES  THb  FMAOTION  OT  T|Nfc  SPfcNT  aT  DETECIOH  J  PEREORHInG  THE  ACTIVITY  DENOTED  BY  K 

MLI  denotes  The  NASS  NEUUCTIOn  EACTUH  IERaCTION  Of  fallout  REHA1NIN6)  FOP  DECUnTANINATION  STRATEOT  l 

STn.II  denotes  hhICH  STHATEQY  L  is  used  to  DECOnTahINATE  plane  I  FOR  COMPOSITE  STRATEttV  M 

PF(J)  denotes  the  PROTECTION  FACTOR  AT  PaCH  UETFCTOR  LOCAIIOh 

kPF(R)  DENOTES  TnE  EOJIVALENT  PROTECTION  FACTOR  FOH  EACH  ACTIVITY  PATTERN 

criJil)  DENOTES  The  FRACTION  OF  The  TOTAL  INTENSITY  PRIOR  TO  DECONTAMINATION  aT  A  OlVEN  DETECTOR 
DOE  TO  A  »AMTICOLAH  contaminated  PLANE 

rnihiJ)  uenotes  The  intensity  reduction  factor  fur  each  ufcontamination  syrateqy  at  Each  detector  location 

HNATHtKI  DENOTES  'HE  ACTIVITY  OOSE  REDUCTION  FACTOR  FOR  EACH  CUMRINaTION  OF  ACTiVITT  AND  DECONTATTIRATTON 
STMaTEQT  at  each  OeTECIOR  location 
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0*4  *■)  S*6  Ca* 

VALUES  or 

ditectoo  plane  n  (j“H  f  TTr 

LOCATION  1834 


NO, 

1 

0, 89900900 

O.iOlOOOOO 

0.92900090  Q.onOOOOO 

NO. 

t 

0.30600000 

0.03100000 

0.20608000  0.37/00000 

**0, 

3 

0.02490000 

U. 20900000 

0.10200000  O.U89000UO 

NO, 

4 

0.08200000 

9.00900080 

0.00100300  8.0028008P 

values 

or 

rR(j,n)- 

• 

DBTECTON 

LOCATION 

1 

2 

3 

4 

5 

NO, 

1 

0.06 

0.0/ 

0,29 

0.13 

0.19 

NO, 

1 

0.0/ 

0.1/ 

0,02 

0.05 

0.08 

3 

0.20 

0.09 

0,03 

0,19 

0.06 

NO, 

4 

0.6/ 

0.6/ 

0.67 

0.67 

0,67 

values 

or 

Ml)** 

1.00 

0.02 

0,03 

0,09 

0.10 

values 

or 

SIM.l)** 

conposite 

ITNaTEQT 

1 

2 

3 

4 

NO, 

1 

1 

9 

1 

2 

NO, 

1 

3 

3 

4 

1 

NO, 

1 

1 

4 

1 

NO, 

4 

1 

3 

6 

J, 

NO, 

9 

4 

3 

1 

1 

NO, 

6 

1 

1 

9 

5 

A  c  T  I  V  I  T  T  p  A^r  r  iT'N  “wu  wTrfir- 


plane  N  U  W  B  E  1 
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Mtj)  • 

-OlTICTOr 

lOCtTlOW 


NO,  1 

1,00 

NO,  t 

1.00 

NO,  S 

2.00 

NO,  A 

100,00 

EoriKi  • 

ACTfVITT 

OATTiON 

NO,  1 

4,22 

NO,  t 

0.43 

NO,  1 

3.11 

NO,  « 

3,02 

NO,  S 

3,26 

crij.iiA 

3ETECT0* 

location 

1 

NO,  1 

0,30 

NO,  1 

0,39 

NO,  I 

0.0» 

NO,  4 

0,20 

CmECA  OV 

SUNNINS 

1.00 

NNiMtJ)  ■ 

CONROSITI 

STRaTIEV 

1 

NO,  t 

0.04 

NO,  t 

0.13 

NO,  ) 

0.02 

NO,  4 

0.44 

NO,  • 

0.A3 

NO.  4 

O.A* 

2 

3 

4 

0.10 

0,03 

0.07 

0.03 

#.ai 

0.3R 

0.42 

0,34 

O.IT 

0.09 

0.10 

0.20 

cr!j,n  ON 

EACH 

1.00 

1,00 

1.00 

2 

3 

4 

0.40 

0,44 

0.30 

0.41 

0.22 

0.23 

O.Sl 

C,42 

0.91 

0.T9 

0,21 

0.43 

0.42 

«.?» 

0.33 

0.44 

0.N2 

0.73 

^L*NE  NUN8Eli 


SmOuLO  Equal  i, 


BETICTO**  location  MUNRIR 
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COMOSITG 

STRATEOr 

1 

2 

3 

4 

9 

RO, 

i 

0.99 

0.64 

0,00 

0.67 

0.73 

RO, 

1 

0.09 

0.31 

0,19 

0.21 

0.21 

RO. 

] 

0.60 

0.72 

0,99 

0.6i 

0.62 

KJ* 

4 

0.40 

0,62 

0,46 

0.46 

0.92 

ROi 

9 

0.91 

0.60 

1.62 

0.60 

0.41 

NO( 

* 

0,90 

0.40 

0,47 

0.49 

0.40 

activity  p  a  t  t  je  n  h  nuh>ep 
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The  Nature  and  Scope  of  ComMind  and  Control  System  Eletnente 
Required  for  Conductlag  Effective  Decontanilnetton 

Analyele 


Note:  The  materiel  in  this  Appendix  was  originally  submitted 
to  USNPJDL  as  Research  Memorandtm  RH>0U>214-9'^. 
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Appendix  F 


The  Nature  and  Scope  of  Comaand  and  Control  System  Elements  Required 
for  Conducting  Effective  Decontamination  In  Municipalities 

I.  INTRODUCTION 


A.  Purpose 

The  purpose  of  this  study  Is  ro  determine  the  nature  and  scope  of  the  comnand 
and  control  system  elements  which  are  required  to  effect  practical  decontamination. 
Emphasis  Is  placed  on  decontamination  within  municipalities.  Such  a  study  Is 
necessary  In  order  to  answer  a  number  of  questions  basic  to  decontamination  analy* 
sis.  The  five  questions  on  which  this  study  focuses  are: 

1.  What  are  the  preattack  and  postrttack  data  required  to  effect  decontami¬ 
nation  operations? 

2.  What  are  the  essential  components  of  the  Information  system  needed  to  ef¬ 
fect  decontamination  operations? 

3.  How  should  trained  and  untrained  personnel  and  decontamination  equipment 
be  pre positioned,  organized,  and  controlled? 

4.  How  can  a  decontamination  system  In  a  municipality  be  evaluated? 

3.  How  can  a  decontamination  system  In  a  municipality  be  most  effectively 
modeled  to  provide  a  ready  vehicle  for  system  analysis? 

Of  coarse,  all  of  these  questions  are  asked  for  various  levels  of  decontamina¬ 
tion  capabilities,  requirements,  and  attack  envlrotoents.  As  implied  by  the  ques¬ 
tions  above,  a  subsidiary  purpose  of  this  study  is  to  develop  a  procedure  for  ana¬ 
lyzing  a  decontamination  system  in  a  municipality. 

B.  Svsisaary 

This  appendix  first  describes  elements  and  purposes  of  coosand  and  control 
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systems  in  general.  The  various  components  of  a  decontamination  system  in  a  municipality 
are  then  identified  and  einbedded  in  the  gereral  command  and  control  system  framework.  It 
is  shown  that  a  connand  and  control  system  for  decontamination  operations  must  pro¬ 
vide  both  for  decisions  on  W.  .ether  or  not  to  undertake  a  mission  and  for  manpower 
and  decontamination  resource  commitment  and  allocation  decisions.  These  decision 
functions  require  an  elaborate  information  subsystem  consisting  of  organized  data 
files  containing  prestored  (preattack)  data  and  postattack  assessments  (including 
system  feedback). 

The  detailed  characteristics  of  the  individual  components  in  each  of  the  essen¬ 
tial  subsystems  of  a  decontamination  command  and  control  system  are  studied.  The 
interrelationships  among  the  individual  subsystem  components  are  identified  and  dis¬ 
played. 

In  order  to  determine  the  conesand  and  control  system  elements  required  to  ac¬ 
complish  practical  decontamination  missions  in  a  municipality,  the  environment  and 
the  system  goals  are  reviewed  and  analyzed.  Basic  system  evaluation  criteria  are 
also  discussed  and  the  essential  decontamination  system  evaluation  criteria  are 
identified. 

Recomnendacions  and  guides  leading  to  the  design  of  a  basic  comiand  and  control 


system  for  municipal  decontamination  are  indicated. 


II.  GENERAL  CHARACTERISTICS  OF  COIMAND  AND  CONTROL  SYSTEMS 


A.  Introduction 

Most  persons  who  are  concerned  with  business,  government,  or  military  operations 
are  also  acquainted  with  the  basic  characteristics  of  coimnand  and  control  systems. 
However,  a  review  of  the  elements  and  purposes  of  command  and  control  systems  is 
useful  in  providing  a  framework  into  which  the  components  of  a  decontamination  sys¬ 
tem  can  be  embedded. 

Further,  since  the  technology  of  command  and  control  systems  is  advancing  rap¬ 
idly,  che  terms  which  define  the  basic  components  of  such  systems  are  changing  day 
by  day.  Thus,  it  is  necessary  to  define  precisely  the  terms  and  expressions  used  in 
'  this  appendix  to  define  the  system  components  and  system  goals. 

The  succeeding  parts  of  this  explanatory  section  are: 

,  Section  II-B  The  Purpose  of  Command  and  Control  Systems 

Section  II-C  The  Basic  Components  of  Command  and  Control  Systems 
Section  II-D  Some  Examples  of  Command  and  Control  Systems 

B.  The  Purpose  of  Command  and  Control  Systems 

The  purpose  of  any  command  and  control  system  is  to:  (1)  determine  for  a  given 
operation  or  set  of  operations,  what  must  be  done,  who  (or  what)  does  it,  and  when, 
where,  and  how  they  do  it;  and  (2)  direct  the  activities  which  accomplish  it. 

Obviously  from  the  above,  it  is  seen  that  a  coosoand  and  control  system  must  have 
the  capacity  both  to  initiate  and  to  monitor  activities.  As  operations  or  sets  of 
operations  become  more  complex,  systematized  cosnand  and  control  mechanisms  become 
more  necessary. 

C.  The  Basic  Coeronenta  of  Coqnand  and  Control  Systems 

Most  command  and  control  systems  are  composed  of  a  set  of  basic  elements  called 
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subsystems.  The  expressions  used  to  identify  these  subsystems  as  well  as  their  de¬ 


fining  characteristics  are  presented  here  along  the  lines  indicated  in  Reference  F-1. 
The  following  subsystems  are  common  to  all  command  and  control  systems: 

1 .  Sensor  subsystem; 

2.  Effector  subsystem; 

3.  External  conmunication  subsystem; 

4.  Internal  cixnmunication  subsystem; 

5.  Information  subsystem; 

6.  Decision  subsystem; 

7.  On-line  control  subsystem;  and 

8.  Off-line  control  subsystem. 

The  coninand  and  control  center  is  composed  only  of  the  Internal  communication 
subsystem,  the  Information  subsystem,  the  decision  subsystem,  and  the  on-line  control 
subsystem.  Figure  F-1  shows  the  functional  relationships  among  the  basic  subsystems 
of  a  command  and  control  system. 

1.  Sensor  Subsystem 

The  basi''  purpose  of  the  sensor  subsystem  is  to  acquire  information  inputs. 
This  informp-ticu  input  m^i’•.erial  is  classified  into  three  categories: 

a.  observed  information  inputs; 

b.  derived  information  inputs;  and 

c.  cotmand  information  inputs. 

Unfortunately,  sensor  subsystems  also  acqu  .re  material  which  is  not  Information 
input  at  all.  This  non- information  input  or  "noise"  represents  Irrelevant  or 
erroneous  material. 

Observed  information  inputs  can  be  such  things  as  visual  observations, 
verbal  communications,  radar  scope  readings,  radiological  dosimeter  readings, 
or  on-alt»  damage  computations.  Derived  information  inputs  are  data  which  re¬ 
sult  from  deductive  processes.  It  may  be  the  product  of  either  experimental 


Figure  F-1  Functional  Oiagran  of  Coonand  and  Control  Syatem 


or  theoretical  research,  for  example.  Thie  Is  usually  prestored  data  and  thus 
the  external  sensors  which  acquire  such  Information  Inputs  may  be  Inactive  dur¬ 
ing  most  of  the  time  that  the  system  Is  operating.  Command  Infoiinatlon  Inputs 
are  In  the  form  of  directives  which  enter  the  system  from  external  sources. 
Decisions  made  outside  of  the  command  and  control  system  but  affecting  the 
operation  of  the  system  are  a  form  of  command  Information  Input.  Figure  F-2 
shows  the  structure  of  Information  Inputs  which  arc  acquired  by  the  sensor  sub¬ 
system. 

2.  Effector  Subsystem 

The  purpose  of  the  effector  subsystem  Is  to  carry  out:  the  directives  of 
the  decision  subsystem.  This  may  consist  of  simply  disseminating  Information 
or  at  the  other  extreme  carrying  out  the  over-all  system  missions.  The 
effectors  are  thus  controlled  by  the  output  conmnds  and  decisions  of  the 
command  and  control  center.  The  effectors'  task  Is  to  apply  the  energy  com¬ 
mitted  for  some  external  action.  An  example  of  an  effector  would  be  an  offen¬ 
sive  aircraft  directed  to  attack  a  target  as  prescribed  In  the  command  and 
control  center. 

3.  External  Comaunlcatlons  Subsystem 

The  external  cocnunlcatlons  subsystm  Is  the  network  of  paths  which 
carries  material  (information  Inputs,  Information,  etc.)  to  and  from  the  func¬ 
tional  blocks  tied  to  thfs  subsystem  In  Figure  F-1.  These  paths  can  consist 
of  radio  communications,  telephone  links,  verbal  Infonaatlon  transmission,  etc. 
Often  an  external  sensor  or  effector  can  also  be  a  part  of  the  external  cosisu- 
nicatlons  subsystem  as  well  as  the  sensor  or  effector  subsystem.  The  organi¬ 
sation  and  ocher  systosm  variables  Involved  in  this  subsystem  depend  largely 
on  the  characteristics  ot  the  particular  cotmaand  and  control  system  of  which 
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Internal  Communlcatlona  Subsyatem 


The  internal  comiDunicationa  subayatem  la  alnllar  to  the  exterral  comnuni- 
cationa  aubayaten  except  that  it  ia  entirely  contained  in  the  conmand  and  con¬ 
trol  center.  Figure  F>1  ahowa  ita  functional  relatlonahlp  to  the  ayatem  aa  a 
whole. 

5.  Information  Subavatem 

The  information  aubayatem  ia  the  part  of  the  command  and  control  center 
which  handlea,  atorea,  tranaforma,  or  diaplaya  the  Information  which  la  trans¬ 
mitted  to  the  conmand  and  control  center.  It  conaiata  of  all  of  the  people, 
machines,  reports,  and  files  within  the  conmand  and  control  center  which  are 
employed  to  perform  this  function. 

6.  Decision  Subsystem 

The  decision  subsystem  (illustrated  in  Figure  F-3)  might  also  be  called  the 
conmand  subsystem.  It  is  the  component  of  the  conmand  and  control  system  which 
tranaforma  information  into  decisions.  These  decisions  are  usually  translated 
into  conmands  to  perform  an  activity,  or  conmit  resources  to  perform  prescribed 
actions  during  a  prescribed  time  interval.  Figure  F-3  illustrates  how  the  de¬ 
cision  subsystem  tranaforma  information  into  decisions  and  coomanda. 

7.  On-line  Control  Subsystem 

Ideally,  the  on-line  control  subsystem  (illustrated  in  Figure  F-4)  controls 
all  information  transmission  and  monitors  all  actions  performed  by  the  conmand 
and  control  system.  In  practice,  the  on-line  control  mechanism  only  attempts 
to  perform  these  tasks  in  moat  real  situations.  Many  activltisa  are  performed 
by  the  effectors  by  their  own  volition  because  there  is  not  enough  time  for  che 
decision  subsystem  to  outline  a  course  of  action.  For  example,  a  Polaris  sub¬ 
marine  under  attack  prestsaably  would  either  attack  its  adversary  or  take  eva¬ 
sive  action  immediately  without  communicating  through  a  higher  echelon  decision 
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subsystem. 

A  very  important  aspect  of  such  on-line  control  is  the  regulation  of  sys¬ 
tem  accuracy  and/or  response  time.  This  means  that  the  on-line  control  sub¬ 
system  attempts  to  control  delays  and  errors  in  the  transmission  and  trans¬ 
formation  of  information  in  the  comnand  and  control  system.  This  role  of  the 
on-line  control  subsystem  makes  it  the  most  important  element  in  the  system. 
Unfortunately,  it  is  also  the  most  difficult  to  identify  precisely  in  most 
physical  systems.  In  the  above  Polaris  example,  the  effector  subsystem  per¬ 
formed  a  function  which,  for  at  least  one  situation,  ioentlfled  the  effector 
as  a  part  of  the  on-line  control  syste  i.  Usually  those  machines,  people,  etc.,, 
which  perform  on-line  monitoring  and  control  functions  from  within  the  command 
i^nd  control  center  are  considered  to  comprise  the  on-line  control  subsystem. 
Figure  F-4  shows  schematically  how  the  on-line  control  subsystem  performs  in 
the  system  environment, 

8.  Off-line  Control  Subsystem 

The  off-line  control  subsystem  monitors  and  controls  processes  which  affect 
the  system  but  are  not  perfonned  on-line.  That  is,  they  are  either  performed 
before  the  system  is  operational  or  are  performed  outside  of  the  control  of  the 
on-line  control  subsystem.  One  important  function  of  the  off-line  control 
subsystem  is  to  collect  and  process  statistics  and  research  for  future  augmen¬ 
tation  of  the  static  variables  of  the  system.  On  the  basis  of  the  observed 
statistics  and  research,  the  off-line  control  subsystem  makes  recamenduclonii 
to  the  on-line  control  subsystem  for  possible  future  control  of  the  dynamic 
variables.  Basic  I'esearch  performed  long  In  advance  of  an  operational  system 
must  be  controlled,  it  is  the  purpose  of  the  off-line  control  subsystem  tc 
control  this  research.  Thus,  most  often,  acme  segments  of  che  off-line  con¬ 
trol  subsystem  ara  the  first  operational  parts  of  a  cgwnsand  and  control  sys¬ 
tem. 
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Figure  F-4  Functional  Characteristic!!  of  On-line  Control  Subsystca 
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D.  Some  Examples  of  Conmand  and  Control  Systems 

Some  of  the  Air  Force's  "L"  Systems  are  excellent  examples  of  military  command 
and  control  systems.  A  few  are  sisnaarlzed  below  (See  Reference  F>2  for  these  and 
additional  examples): 

412L  -  Air  Weapons  Control  System 

An  overseas  theater  tactical  air  weapons  control  and  warning  system. 

416L  •  SAGE  Air  Defense  System 

A  semi-automatic  area  air  weapons  control  and  warning  system  including 
the  Back  Up  Interceptor  Control  (BUIC)  for  detecting,  identifying, 
cracking,  and  providing  interceptor  weapon  direction  against  air- 
breathing  threats  to  the  United  States  and  Canada. 

474L  -  Ballistic  Missile  Early  Warning  System 

A  system  to  provide  early  warning  of  a  mass  ballistic  missile  attack  on 
the  North  American  continent  from  the  North. 

4811.  -  Postattack  Coranand  and  Control  System 

A  system  to  enable  the  Conmander- In- Chief,  Strategic  Air  Conmand,  to 
control  his  forces  In  the  event  that  a  nuclear  attack  destroys  or 
seriously  degrades  his  normal  facilities. 

Many  of  these  as  well  as  other  of  the  *’L"  Systems  are  subsystems  of  higher 
echelon  command  and  control  systems.  474L  (BMEWS),  described  above,  is  also  a 
sensor  subsystem  of  higher  echelon  systems. 
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III.  A  COMMAND  AND  CONTROL  SYSl’EM  FOR  KUrilCIPAL  DECONTAMINATION 


A.  General 

The  previous  section  defined  the  characteristics  of  subsystem  interrelationships 
conmon  to  all  coiranand  and  control  systems.  It  is  the  purpose  of  this  section  to 
identify  the  basic  system  elements  required  by  municipal  decontamination  and  to 
embed  these  elements  into  the  general  connand  and  control  systems  framework  of  the 
previous  section.  In  order  to  accomplish  this,  the  following  steps  are  taken  in 
this  section: 

1.  Analyze  the  enviroianent,  within  which  municipal  decontamination  is  expected 
to  operate. 

2.  Identify  the  goals  of  municipal  decontamination. 

3.  Prescribe  the  decontamination  system  performance  requirements. 

4.  Identify  and  relate  the  elements  of  a  decontamination  system. 

5.  Relate  the  elements  Identified  In  4  above  to  the  general  cooEiand  and 
control  system  framework  defined  in  Section  II. 

B.  Environmental  Aspects  of  Municipal  Decontamination 

Municipal  decontaraination  by  itself  is  only  a  part  of  the  broadsr  system  of 
postattack  radiological  defense.  The  time  during  which  a  radiological  defense  sys¬ 
tem  operates  can  be  logically  divided  into  th^ee  time  phases:  the  emergency  phase, 
the  operational  recovery  phase,  and  the  final  recovery  phase  (Reference  F-3).  The 
underlying  reason  for  dividing  radiological  defense  into  three  time  phases  is  the 
change  due  to  fallout  decay. 

Figure  -5  shows  that  the  emergency  phase  begins  with  a  preattack  warning--an 
event  highly  probable  during  tines  of  International  crisis-- if  it  occurs,  or  attack 
warning  otherwise.  It  may  last  for  several  days,  weeks,  or  month*  depending  on  the 
warning  obtained  icoakimd  with  the  population  discipline  attained  for  such  an  event) 
and  the  attack  psrameter*.  The  decision-makers  during  this  phase  utilize  preattack 
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planning,  before  and  after  warning,  the  extent  and  effectiveness  of  which  is  de¬ 
termined  by  the  amount  of  preattack  warning  given,  'fhis  is  a  function  of  intelli¬ 
gence  indicators. 

The  objective  of  radiological  countermeasures  during  this  phase  is  survival,  for 
it  is  durixig  this  period  that  the  fallout  arrives,  accumulates,  attains  some  maximum 
intensity,  and  then  begins  to  decay.  It  has  been  shown  that  the  optimal  or  primary 
radiological  countermeasure  during  this  period  is  enshelterment.  The  central  re¬ 
quirement  of  the  system  during  this  period  is  to  provide  a  sufficient  number  of 
adequate  shelters  so  located  as  to  minimize  casualties.  The  potential  effectiveness 
of  the  system  during  this  period  is  determined  by  the  degree  of  protection  afforded, 
space  availability  and  accessibility,  and  warning  time. 

The  operational  recovery  phase  follows  the  emergency  phase  and  encompasses  the 
operations  prerequisite  to; 

!•  Sustaining  life  in  a  hostile  environment;  and 

The  recovery  of  essential  postattack  activities  and/or  facilities. 

From  Figure  F-5  it  is  seen  that  the  transition  to  the  operational  recovery  phase 
takes  place  when  limited  egress  from  shelters  is  feasible. 

Tlie  final  recovery  phase  begins  when  radiation  intensity  decays  to  an  insignif¬ 
icant  level  and  primal^  considerations  can  be  focused  on  functional  restoration  of 
the  target  area  as  nearly  as  possible  to  its  preattack  condition. 

Municipal  decontamination  will  be  activated  primarily  during  the  operational 
recovery  phase.  The  following  environmental  conditions  will  then  prevail: 

1.  Fallout  deposition,  generally,  will  be  complete. 

2.  Fallout  radiation  will  still  be  a  hazard  in  some  areas.  Host  of  the 
civilian  population  will  be  sheltered  during  much  of  the  time  that  de¬ 
contamination  is  being  performed. 

3.  There  may  be  damage  In  the  area  from  direct  weapons  effects  (blasts, 
thermal,  EMP,  etc.). 
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4.  The  entire  Civil  Defence  System  will  be  very  active,  i.e.,  communications 
lines,  logistics  supplies  and  equipment,  etc.  will  be  very  busy. 

This  last  item  is  perhaps  the  most  important  condition  affecting  the  environment 
of  municipal  decontamination  systems.  It  is  a  factor  over  which  the  decision-maker 
has  some  control.  As  is  pointed  out  in  Reference  F-3,  the  decision-maker  has  little 
or  no  control  over  most  environmental  factors. 

Reference  F-3  also  points  out  that  the  environment  is  an  important  factor  lim¬ 
iting  the  range  of  decisions  which  can  be  made  in  the  postattack  period.  Because 
of  the  wide  variation  of  damage,  yielding  a  wide  variety  of  physical  environments 
(in  which  municipal  decontamination  might  operate),  the  decision  subsystem  of  any 
proposed  decontamination  system  must  be  extremely  flexible  and  capable  of  making 
very  complex  decisions.  Further,  it  should  be  capable  of  making  decisions  based  on 
varying  amounts  of  reliable  (and  unreliable)  information.  This  will  subsequently 
be  discussed  in  some  detail. 

C.  Goals  of  Municipal  Decontamination 

The  primary  goals  of  decontamination  within  a  municipality  (in  a  very  broad 
sense)  are: 

1.  Continued  survival  within  the  municipality, 

2.  Accelerated  municipal  recovery,  and 

3.  Assistance  In  the  recovery  of  a  region  (not  necessarily  containing  the 
municipality  '.’here  the  decontamination  systm  is  operating). 

Since  the  three  goals  listed  above  are  also  primary  goals  of  the  overall  radio¬ 
logical  defense  system,  any  decontamination  system  must  be  closely  coordinated  v’ith 
other  radiological  defense  systems.  This  Implies  the  existence  of  an  overall  radio¬ 
logical  defense  comaand  and  control  system.  Although  it  is  not  the  purpose  of  this 
report  to  elaborate  on  the  overall  system,  It  should  be  recognized  that  much  of  the 
coeinand  information  input  to  the  decontamination  counand  and  control  system  will  be 
output  from  such  a  higher  echelon  decision  subsystem. 
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D.  Municipal  Decontamination  System  Performance  Criteria 

In  light  of  the  existence  of  a  higher  echelon  command  and  control  system,  the 
question  might  be  asked:  "Why  a  command  and  control  system  for  decontamination?" 
There  are  two  basic  reasons  for  proposing  a  separate  system  to  command  and  control 
decontamination  opera'^ions  in  a  municipality: 

1.  The  specific  decisions  involved  with  effcctiig  decontamination  operations 
in  a  municipality  are  in  themselves  complex  enough  to  warrant  a  "decontami¬ 
nation-level"  decision  subsystem. 

2.  The  information  requirements  for  effecting  useful  decontamination  decisions 
are  in  themselves  very  involved  and  require  a  separate  "decontamination- 
level"  information  subsystem.  This  is  not  to  say  that  a  decontamination 
information  subsystem  would  not  borrow  components  from,  or  overlap  with, 
parts  of  other  radiological  defense  information  subsystems,  e.g.,  RADEF. 

I.ater  in  this  appendix,  the  information  input  and  decision  structure  requirements 
for  effecting  decontamination  will  be  identified  and  will  substantiate  the  above 
reasoning  for  a  separate  decontamination  command  and  control  system. 

The  basic  criteria  upon  which  to  judge  a  decontamination  system  are: 

1.  To  what  extent  does  the  system  Increase  the  effectiveness  of  decontamination 
in  achieving  the  goals  and  requirements  prescribed  by  higher  echelon  decision 
subsystems? 

2.  How  effectively  does  the  system  reduce  on-line  decontamination  costs  (rn 
manpower,  crew  doses,  expenditures  of  fuel,  water,  etc.)  for  specified 
levels  of  performance? 

3.  How  effectively  does  the  system  reduce  off-line  decontamination  costs  (re¬ 
search  funds  required,  training  of  personnel,  procureiwnt  of  equipment,  etc.) 
for  specified  levels  of  performance. 

Tltis  last  criterion  is  separated  from  the  second.  In  that  it  can  be  treated  as 
a  static  goal  of  .he  system  (i.c.,  one  which  is  evaluated  by  off-line  simulation 
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and/or  analysis  of  on-line  system  performance).  The  first  two  are  dynamic  goals 
(i.e.,  those  evaluated  by  on-line  performance  and  vrhich  can  be  controlled  on-line). 

E.  Elements  of  a  Decontamination  System 

The  purpose  of  this  section  is  to  identify  the  functional  decontamination  de¬ 
cision  elements  as  well  as  to  determine  the  data  base  and  information  flow  required 
by  these  elements.  These  data  and  information  flow  may  be  internal  to  the  decon- 
t.smination  system  (between  subsystems)  or  external  to  the  system  (to  and  from  higher 
echelon  or  para'/ lei  systems).  The  higher  echelon  command  and  control  system  with 
regard  to  decontamination  is  the  Postattack  Civil  Defense  Emergency  Operating  Sys¬ 
tem.  Parallel  systems  would  include  systems  to  control  such  functions  as  rescue, 
law  and  order,  engineering,  welfare,  firefighting,  medical,  etc.  Certain  other 
subsystems  of  the  Emergency  Operating  System  would  provide  sensor  information  inputs 
for  the  decontamination  system.  These  include  RADEF,  damage  assessments,  NUDETS, 
etc.  Figure  F-6  illustrates  schematically  the  information  flow  between  the  various 
command  and  control  systems  related  to  decontamination. 

The  subsystems  of  the  decontamination  system  must  provide  the  basis  for  command¬ 
ing  and  controlling  all  of  the  activities  and  decisions  which  comprise  decontamina¬ 
tion.  Thus,  a  first  step  towards  identifying  the  elements  of  a  decontamination 
command  and  control  system  is  to  list  all  of  the  important  activities  and  decisions 
which  are  required  to  perform  effective  decontamination.  In  order  to  determine 
whether  or  not  to  decontaminate  a  given  facility  or  activity  and  thence  to  schedule 
the  operation,  Che  following  steps  must  be  taken: 

1.  Determine  Chat  the  facility  or  activity  is  essential  to  sustaining  life  or 
accelerating  recovery.  Tills  may  Involve  no  more  than  being  cold  by  the 
higher  echelon  decision  subsystem  that  a  given  facility  or  activity  is 
essenC iai . 

2.  Determine  that  the  given  facility  or  activity  is  presently  denied  bv  fall¬ 
out  or  will  be  when  you  want  it. 
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3.  Determine  that  dacontamlnation  can  reduce  the  fallout  denial  time  of  this 
facility  or  activity  to  or  below  the  minimum  desired. 

4.  Determine  that  equipment  and  supplies  (water,  etc.)  are  available  to  decon¬ 
taminate  around  this  facility  or  activity. 

5.  Determine  that  sufficient  manpower  (taking  dosa  histories  into  account)  is 
available  to  decont.iminate  this  facility  or  activity. 

6.  Determine  that  the  facility  or  activity  can  be  decontaminated  to  a  specified 
level  of  effectiveness  without  overexposing  the  decontamination  crews. 

7.  Determine  that  direct  weapons  effects  or  congestion  will  not  impede  the  de> 
contamination  equipment  in  reaching  the  facility  or  activity. 

8.  Decide  to  decontaminate  the  given  facility  or  activity. 

9.  Commit  resources  and  schedule  the  decontamination  of  the  given  facility  or 
activity. 

10.  Initiate  and  conduct  the  decontamination  operation. 

11.  Monitor  and  control  the  decontamination  operation. 

12.  Determine  whether  additional  resources  are  required  and  available  to  complete 
the  decontamination  of  the  facility  or  activity. 

13.  Conmit  additional  resources  to  decontaminate  the  facility  or  activity  if 
needed . 

14.  Determine  that  the  ct'*contamlnacion  operation  is  ccmpleted. 

Each  of  these  activities  triggers  a  response  within  the  system.  Thus,  none  of 
these  activities  can  be  omitted  in  the  analysis  without  reflecting  a  gap  in  the  re- 
qulreu  data  base  and  organUat lonal  structure.  Each  activity  demands  a  supporting 
data  base,  some  resources,  and  some  organisational  (coosand  and  control)  structure. 
This  organisational  structure  should  provide  for  both  Information  channels  and  de¬ 
cision  itroctuTcs. 

Figure  F-7  show*  the  supporting  structure  and  data  bate  for  conducting  the  de¬ 
contamination  of  a  given  facility  or  activity.  Each  of  the  important  activities 
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involved  with  such  an  operation  is  delineated  and  tlie  data  base  fov  each  activity 
specified.  Clearly  the  required  information  channels  comprise  a  complex  network 
and  imply  the  necessity  for  a  well  organized  information  subsystem. 

Close  examination  of  Figure  F-7  reveals  the  necessity  for  other  subsystems 
iiimllar  to  those  described  in  Section  II  of  this  report.  For  instance,  an  on- lint 
control  system  will  be  required  to  monitor  and  control  the  decontamination  operations 
and  also  to  direct  and  insure  that  information  is  forwarded  to  the  proper  functional 
elements  of  the  decontamination  system.  The  parallelism  between  required  decontami¬ 
nation  system  w.Iements  and  classical  coimand  and  control  subsystems  is  considered 
next. 

F.  Decontamination  System  elements 

The  purpose  here  is  to  identify  the  command  and  control  system  elements  re¬ 
quired  to  conduct  effective  municipal  decontamination.  The  interrelationships  be¬ 
tween  the  functional  elements  are  discussed  and  displayed. 

In  Section  II  of  this  report,  each  of  the  necessary  subsystems  of  any  command 
and  control  system  are  identified  and  defined.  The  basic  supporting  data  base  as 
well  as  the  activities  and  decisions  necessary  to  conducting  postr.ttaci'  decontamina¬ 
tion  are  illustrated  in  Figure  F-7.  It  is  nov'  necessary  to  relate  these  necessary 
decontamination  system  activities  to  the  command  and  control  framework  described  in 
Section  II.  Table  F-I  lists  some  of  the  basic  dif. contamination  functions  performed 
by  each  of  the  cooanand  and  control  subsystems.  Figure  F-8  Illustrates  the  functional 
relationships  among  these  decontarolnatton  activities.  It  is  seen  that  some  of  these 
functions  are  also  parallel  systems  to  decontamination.  Radiological  monitoring, 
for  Instance,  would  be  required  for  almost  all  of  the  parallel  systems  shown  In 
Figure  F-b. 

What  i#  possible  to  note  from  Tsble  F-1  and  Figures  F  7  and  F-8.  however,  is  that 
ayai.'TO  r<K|uired  to  conduct  effective  dr  ontaminatlon  does  have  all  of  the  character¬ 
istics  of  a  full  -•itcale  coffansr-.d  and  control  cystea.  All  of  the  usual  sub-systems  of 
ctwsmand  and  control  *-yste«s  are  needed  In  some  fashion  to  perform  municipal  decon¬ 
tamination. 
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TABLE  F-I 


Decontamination  Subsystem  Functions 


Sensor  Subsystem 

1.  Radiological  monitoring. 

2.  Field  Surveys,  etc. 

Effector  Subsystem 

1.  Move  decontamination  equipment  to  sites. 

2.  Perform  decontamination. 

External  Communications  Subsystem 

1.  Transmit  data  and  Information  to  and  from  radiological  monitors,  field 
Inspection  teams,  etc.,  and  decontamination  crews  to  the  Emergency 
Operating  Center  (EOC). 

Internal  Comnunlcatlons  Subsystem 

1.  Transmit  data  and  Information  within  the  EOC. 

Information  Subsystem 

1.  Score  and  catalog  decontamination  planning  guides  and  other  related 
(prestored)  materials. 

2.  Display  environmental  information  (RADEF,  etc.)  within  the  EOC. 

3.  Process  data  Inputs  and  information  coming  in  from  the  external  sensors 
into  a  usable  form  for  persons  with  Information  needs. 

4.  Provide  Information  to  decision  subsystem  on  request. 

Decision  Subsystem 

1.  Decide  when,  what,  and  how  to  decontaminate. 

2.  Decide  when  to  conmlc  .'-esources  and  which  resources  to  comnlt. 

3.  Decide  whether  additional  Information  is  required  to  make  a  decontami¬ 
nation  decision. 

On-line  Control  Subsystem 

1.  Monitor  and  control  decontamination. 

2.  Coomit  resources  to  decontamination  tasks. 

3.  Control  InfirmaClon  flow,  assign  tasks  to  radiological  monitors  and  field 
Inspection  teams.  This  Included  regulating  response  times  and  information 
accuracy. 

Off-line  Control  Subsystem 

1.  Monitor  and  control  decontamination  research. 

2.  Preposition  and  inventory  decontamination  personnel,  equipi«nc,  and 
supplies. 

3.  Organise  decontmninatlon  subsystems  within  the  EOC. 


« 
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Figure  F-8  Illustration  of  Functional  Relationships  Among 

Decontamination  System  Activities 


IV.  AN  INFORMATION  SYSTEM  (SUBSYSTEM)  FOR  DECONTAMINATION 


A.  Introduction 

Figure  F-7  in  the  last  section  depicts  a  very  complex  information- flow  network 
required  in  the  decision  whether  or  not  to  conduct  a  particular  decontamination 
operation.  Certainly  some  control  over  the  collecting,  scoring,  and  disseminating 
of  all  of  these  data  is  necessary.  It  is  Che  purpose  of  this  section  to  indicate 
the  basis  upon  which  an  information  subsystem  for  decontamination  could  be  designed. 
The  steps  taken  towards  this  end  are  along  the  lines  suggested  in  Chapter  3  of 
Reference  F-3: 

1.  The  information  input  requirements  are  discussed. 

2.  The  information  format  including  the  criteria  upon  which  to  judge  the  use¬ 
fulness  of  Che  information  are  discussed. 

3.  The  means  for  obtaining,  storing  and  disseminating  Che  required  information 
are  identified. 

The  last  two  steps  required  in  des:.gnlng  an  information  system  (identifying 
alternative  systems  and  choosing  a  best  one)  are  not  taken  up  in  this  report. 

B.  Information  Input  Requirements  for  Conducting  Decontamination 

Most  of  the  information  inputs  shown  in  Figure  F-7  are  not  required  for  every 
decontamination  decision.  For  instance  if  no  blast  damage  has  occurred  in  the 
city  where  decontamination  is  being  performed,  specific  blast-damage  information 
inputs  would  be  unnecessary.  It  a  higher  echelon  decisions  subsystem  commands  that  a 
facility  be  decontaminated,  the  information  inputs  required  to  determine  Chat  the 
job  is  necessary  are  no  longer  relevant  to  the  decisi'>n.  Thus,  there  are  different 
infonuitlon  input  requirements  depending  on  the  environment  and  the  kind  of  decon¬ 
tamination  decision  being  made.  The  environment  Includes  not  only  the  attack  damage 
but  rhi:  overall  decontamination  resource  situation  and  the  population  dispersal. 

(Tliis  laj^t  item  is  a  function  of  the  attack  warning  time  of  day  of  attack,  pre- 
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conditioning  and  training,  etc.)  Clearly  the  problem  of  identifying  all  of  the  pre¬ 
attack  and  po^tattack  data  requirements  for  effecting  municipal  decontamination  for 
all  situations  would  be  a  mammoth  undertaking.  VHiether  or  not  such  an  effort  would 
be  fruitful  is  in  itself  a  difficult  question.  This  report  will  simply  list  some  of 
the  basic  requirements  which  appear  to  be  a  necessary  input  to  most  decontamination 
decisions.  They  are: 

1.  Fallout  Information  -  This  includes  both  amounts  of  fallout  (mass)  and  the 
dose-rates  where  the  fallout  is  deposited. 

2.  The  location  of  essential  facilities  and  the  urgency  of  their  recovery  - 
This  information  may  be  in  the  form  of  a  command  from  a  higher  echelon 
decision  subsystem  or  prestored  Information  in  the  EOC. 

3.  Decontamination  resource  inventories  and  schedules  -  This  information 
should  be  kept  up  to  date  in  the  EOC. 

4.  Planning  guides  for  conducting  decontamination  -  This  information  is  pre¬ 
stored  in  the  EOC. 

These  basic  information  requirements  are  presented  in  a  general  manner  only; 
any  additional  refinement  would  require  a  complete  situation  analysis  including  a 
description  of  the  environment  and  the  kind  of  decontamination  decision  being  made. 

C .  Information  Utility  Criteria 

As  stated  in  Reference  F-3,  the  basic  criteria  upon  which  to  Judge  information 

are: 

1.  Timeliness; 

2.  Pertinence; 

3.  Reliability;  and 

4.  Comprehensibility 
1 .  Timeliness 

The  value  of  radiation  information  particularly  is  very  dependent  on  when 
it  is  received.  This  is  because  of  the  decay  rate  of  the  fallout  intensity. 
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Radiation  reports  during  the  build-up  phase  and  shortly  thereafter  are  not 
essential  for  future  decontamination.  This  information,  in  fact,  is  best  kept 
out  of  the  decontamination  information  subsystem.  It  simply  clutters  the  system. 
Radiation  information,  to  be  useful,  must  be  reasonably  up  to  date  and  reflect 
the  amount  of  fallout  or  intensity.  As  the  decay  progresses  it  becomes  less  and 
less  important  that  the  information  be  processed  quickly.  Of  course  late  re¬ 
ports  might  delay  a  decision  to  decontaminate  an  essential  facility. 

Prestored  data  should  be  organized  so  that  the  decision-maker  can  obtain 
information  from  his  files  as  he  needs  it. 

2.  Pertinence 

The  Information  should  be  organized  and  disseminated  in  a  manner  to  minimize 
the  amount  of  u’mecessary  information  being  pro 'ided  to  the  decontamination  de¬ 
cision-maker.  Prestored  information  (decontamination  planning  guides,  etc.) 
should  be  concise  and  to  the  point.  As  is  pointed  out  in  Reference  F-3,  "Extra¬ 
neous  material  not  only  wastes  the  time  of  those  who  must  use  it,  it  may  also 
lead  to  incorrect  decisions  by  obscuring  the  most  important  considerations." 

3.  Reliability 

Information  should  be  accurate.  Reference  F-3  points  out  chat  accuracy  is 
influenced  by  Che  quality  of  the  information  inputs  as  well  as  the  ability  of 
the  system  to  transmit  and  transform  data  without  making  errors. 

Unfortunately,  this  criterion  is  particularly  difficult  to  meet  for  decon¬ 
tamination  decisions.  Radiological  monitoring  devices  are  quite  variable  in 
their  accuracy  of  intensity  measurements.  Indeed,  winds  or  redistributed  fnllout 
cause  information  inputa  to  reflect  incorrectly  the  status  of  the  radiological 
environment.  Because  of  this  uncontrollable  rellabtUcy,  decontamination  de¬ 
cisions  will  have  to  be  made  despite  aoet*  inaccuracy  in  infonsatlon  inputs. 
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4.  Comprehensibility 


Information  inputs  from  external  sensors  must  be  collected  (time,  location, 
etc.)  and  presented  in  a  comprehensible  manner.  Prestored  data  must  be  clear 
and  written  or  presented  in  a  manner  readable  by  the  user,  e.g.,  decontamination 
planning  guides  should  be  readable  to  the  local  decision-maker  in  the  £0C. 

D.  Means  for  Obtaining.  Storing,  and  Disseminating  the  Information 

Radiation  information  can  be  obtained  from  rt  liation  monitors.  This  informa¬ 
tion  can  be  measured  with  instruments  or  visually  observed.  It  may  be  forwarded  to 
the  EOC  via  telephone,  radio,  or  hand-carried.  Such  information  can  be  posted  and, 
kept  up  to  date  on  display  boards,  map  books,  or  typed  and  hand-written  reports. 

Field  survey  data  can  be  handled  in  the  same  manner 

Prestored  data  are  obtained  by  preattack  uata  acquisition  and  research.  The 
kinds  of  prestored  data  obtained  are  determined  by  careful  preattack  assessment  of 
postattack  data  needs.  They  can  be  stored  as  printed  reports,  maps,  etc.,  or  if 
retrie  al  machinery  is  available,  they  could  be  stored  on  punched  cards,  magnetic 
tapes,  drums,  discs,  etc. 

C.  Sumnary 

It  is  not  the  purpose  in  this  appendix  to  prescribe  preclse'i  '  an  information 
subsystem  for  decontamination.  Whether  or  not  a  separate  information  subsystem  for 
decontamination  as  opposed  co,  say,  firefighting  is  needed  is  not  answered  in  this 
apiendix  The  answer  to  this  question  dependit  largely  o.  the  organizational  structure 
and  systeti  design  of  whatever  higher  d'helon  systems  are  developed  (e.g..  Emergency 
Operating  Systems). 

This  appendix  simply  indicates  a  basis  for  dosigning  a  decontamination  informa¬ 
tion  system. 
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V.  CONCLUSIONS  AND  RECOMMENDATIUNS 


A.  Conclusions 

Some  command  and  control  structure  will  be  necessary  to  conduct  effective  munici¬ 
pal  decontamination.  VRiether  or  not  a  separate  decontamination  command  and  control 
system  is  necessary  is  not  determined  in  this  report--nor  will  it  be  determined  with¬ 
out  a  far  more  extensive  research  effort.  The  decisions  and  activities  related  to 
decontamination  do,  however,  require  an  elaborate  information  subsystem  consisting 
of  organized  data  flies  containing  prestored  (NFSS  data,  decontamination  planning 
guides,  etc.)  data  and  postattack  assessments  (including  system  feedback).  A  good 
argument  for  a  decontamination  information  subsystem  could  be  presented  on  the  basis 
of  the  complex  network  of  Information  flow  shown  in  Figure  F-7. 

B.  Recommendations 

It  is  recommended  that  decontamination  be  considered  in  the  development  and  de¬ 
sign  of  any  Postattack  Emergency  Operating  System  to  be  operdtional  at  the  municipal 
or  community  level.  Ic  is  further  recommended  that  the  systems  requirements  for  con¬ 
ducting  decontamination  be  analyzed  more  thoroug  :ly  and  the  necessary  system  functions 
and  minimal  organizational  structure  for  conducting  effective  municipal  decontamina¬ 
tion  be  more  precisely  identified. 
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